■  •••  i 

- 

US  Army  Corps 
of  Engineers 


mu# 


MISCELLANEOUS  PAPER  GL-87-25 


AN  EXAMINATION  OF  SLOPE 
STABILITY  COMPUTATION  PROCEDURES 
FOR  SUDDEN  DRAWDOWN 


Stephen  G.  Wright 

Geotechnical  Engineering  Department 

-\V 

•  "J 

University  of  Texas  at  Austin 

■*.  -O 

Austin,  Texas  78712 

and 

J.  M.  Duncan 


Department  of  Civil  Engineering 
Virginia  Polytechnic  Institute 
Blacksburg,  Virginia  24061 


Jr  f - 


3  \  '  i  z 

\  .•  A*/ 


W0^0  6  1987 fj  1 


September  1987 
Final  Report 


Approvt'ff  f(if  Public  Roh\iS<‘  f ) i S T r ;hu!n  in  (JrFn  N 


Prep.m,)  f(),  DEPARTMENT  OF  THE  ARMY 
US  Army  Corps  of  Engineers 
Washington,  DC  20314-1000 

Umjei  Computer  Applications  in 
Geotechnical  Engineering  (CAGE)  Proiect 
Contract  Nos.  DACW39-86-M-0982.  DACW39  86-M-1  73 ' 


m. .r.itor,-. j  hy  Geotechnical  Laboratory 
US  Army  Engineer  Waterways  Experiment  Station 
PO  Box  631,  Vicksburg.  Mississippi  39180-0631 


87  i-'  ■  '  -> 


Destroy  this  report  when  no  longer  needed.  Do  not  return 
it  to  the  originator. 


The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position  unless  so  designated 
by  other  authorized  documents. 


The  contents  of  this  repoi'  are  not  to  be  used  for 
advertising,  publication,  or  promotional  purposes. 
Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of 
such  commercial  products. 


Unclassified 

SECURITY  CLASSIFICATION  OF  THIS  PAGE 


la  REPORT  SECURITY  CLASSIFICATION 

Unclassified 

2a  SECURITY  CLASSIFICATION  AUTHORITY 
2b  DECLASSIFICATION /DOWNGRADING  SCHEDULE 
4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 


REPORT  DOCUMENTATION  PAGE 

Ton  I  lb  RESTRICTIVE  MARKINGS 


form  Appro* ect 
OMfl  No  0W4  0188 
£*p  Date  iun  )0  198b 


3  DISTRIBUTION 'AVAILABILITY  OF  REPORT 

Approved  for  public  release; 
distribution  unlimited 

5  MONITORING  ORGANIZATION  REPORT  NuMBfRIS) 

Miscellaneous  Paper  GL -87-25 


6a  NAME  OF  PERFORMING  ORGANIZATION  6b  OFFICE  SYMBOL  T  7a  NAME  OF  MONITORING  ORGAN'ZA'ION 

See  reverse  O' appi, cable)  USAEWFS 


inapplicable)  (JSAFWFS 

Geotechnical  Laboratory 

7b  ADDRESS  (City.  State  and /ip  Code) 

PO  Box  631 

Vicksburg,  MS  39180-0631 

8b  OFFICE  SYMBOL  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NuMBE  R 

(If  appi, cable)  DACW39-86-M-0982 

_ DACW39-86-M-1731 _ 

10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM  I  PRC 

ELEMENT  NO  I  NO 


6<  ADDRESS  (City,  State,  and  ZIP  Code) 

Austin,  TX  78712  and 
Blacksburg,  VA  24061 

8a  NAME  OF  FUNDING  .  SPONSORING 
ORGANIZATION 

US  Army  Corps  of  Engineers 

8c  ADDRESS  (City.  State,  and  ZIP  Code) 

Washington,  DC  20314-1000 


11  title  (Include  Security  Claudication) 

An  Examination  of  Slope  Stability  Computation  Procedures  for  Sudden  Drawdown 

12  PERSONAtTauThOR(S}  _  _________ 

Wright,  Stephen  G. ,  and  Duncan,  J.  M. 

1  3a  TYPE  of  REPORT  Mb  T iME  COVE  RE D_  DATE  qu  REPORT  (Year  Nronth  Oay)~  1 S  f 

final  Report  from  to  September  19 87 

16  supplementary  notation 

Available  from  National  Technical  Information  Service,  SR25  Port  Royal  Road, 
Springfield,  VA  21*  161. 


PROJECT 

TASK 

NO 

NO 

WORK  jNiT 
ACCESSION  NO 


14  DATE  OF  REPORT  { Year  toonth  Day)  IS  PAGF  COUNT 

September  19 87  177 


1  7 

COSATl 

CODE  S 

10  SuBifC  T  r  E  R  VS  Continue  on  reverse 

it  neressary  and  identify  by  block  number) 

FIELD 

GROUP 

SUB  (jROuP 

t  f  t » ■  <.  t  i  v . >  s  t  res  s  and  lysis 

S i ope  s tab i  1  i  f  v 

i‘orr  pef'SSure  ujeffii  ii-nfs 

Sudden  drawdown 

Shear  Sti'E.'TiUth  enveloiies 

’wu-staoe  j t i a  1  V  a  1 

19  ABSTRACT  i Continue  on  reverse  >f  necessary  and  identify  by  blcxk  number) 

“Sudden"  or  "rapid"  drawdown  is  considered  to  oi<u»  wUcii  •  vn  i  ••  n»  nth**  ad  ■  ,i . 
water  levels  die  lowered  at.  such  a  rate  that  little  nr  rm  drainage  n< <  <n  \  n;  ,im  *  * , » •  t  *  • 

slope  at  the  time  the  water  1  ♦  •  v  *  •  1  is  being  lowered.  1  :ipe  t  ah  i  1  ’  t  v  analyse*.  a»‘e  •  *  *  . ' 

performed  to  c.aliulat**  the  tac  tor  <>♦  safety  for  earth  -J.  ,..*•<  tJb  ie<  t  ,•  1  t ,  *  t  *■  . . ip 

Several  different  prm  eOures  e»ist  and  air  ureentl,  used  by  des  ;  :*!*■*•••.  t.n  iir  .  ,*> 
slope  stability  ti.11  the  I  use  nt  sudden  drawdown.  T  he*.»  vlt’n.js  ,  *  .•••In*-,  j**-  Im  •  ' 

t  undai'ien  t  a  1  1  y  d  ?  f  f  e  r  epi  f  app  r/  i,i'  hr  * ;  and  ire  known  to  ;  •  r  <  !  ,i  »  • !  i f  1 » • »  *  a,  t  »  •  ■  , 1  1  ■  1  t  *  ’  * 


some  instan  i 

obtained  usi 

1  he  var 
I  S  <1ef  ’  fil'd  a 

.)■  e.hi'S-  .  » 


I  tie  Of  111. 
'.{  the-,e  ; ir'( <* 
l *  u ’ .  ;>  i  i  ii  edm 

d  i  ac  tie  p  i ) 
a  f ■ '  *v t ;  1 1  *  ii' 
,u  1  re  t  i  :  ♦  ; 


1  i  ir  'S  ,  the  e.  ,  j"  p  t  1 1 » f  i  *  *i'  p  i  n  v 1  d  ,  aid  t  ►  1 1  ■  *  \  • ,  ♦  t  1  <  a  ! 
1u»  es  air*  1  he  '  uti  I  •  ♦  *• ;  •.  .  a.  *•*  . 

.  '  1  1  f  t  » ■  r  p  r  Ml.  !  p  a  i  V  in  Ml.  ■  .  t  * .  II .  • »  •  I  A  *  i  ■  f  T  *  •  * 
I'd  l  ’it  i )  t  wm  ; *  ■ 1  •  *  a  •  1  t  i  •  J-  • »  l  •  •  .  e ,  *  *  •  *  ;  *  t  . 

! -a ne  lh  ,  .  ■  .  T  t  ...  •  ,t  '  »  .  .  '  M  ‘ . 

1  i  ’  vi  1 1 » ■  f  ,  i  s  ■  *  •  !  ’  .  i f  *  ■  !  •  T  •*..’;  *  ,  ! 


DO  FORM  147 J,  hj  va- 


mcuwty  CWAMine*T10W  ow  rm>  >»■«  _ 


6a.  NAME  OF  PERFORMING  ORGANIZATION. 

Geotechnical  Engineering  Center,  Civil  Engineering  Department,  University  of  Texas 
at  Austin  and  Department  of  Civil  Engineering,  Virginia  Polytechnic  Institute 

19.  ABSTRACT  (Continued). 

and  the  differences  among  the  various  effective  stress  procedures  are  related  directly  to 
differences  in  the  procedures  used  to  estimate  pore  water  pressures.  The  second  category 
of  procedures  for  sudden  drawdown  analysis  is  comprised  of  two-stage  procedures:  one  set 
(stage)  of  calculations  is  performed  for  conditions  existing  immediately  prior  to  draw¬ 
down  and  is  used  to  estimate  effective  stresses  and  corresponding  undrained  strengths;  the 
second  set  (stage)  of  computations  is  performed  for  conditions  immediately  after  drawdown 
and  employs  undrained  shear  strengths  estimated  from  the  first  set  of  computations. 

This  study  provides  a  basis  for  a  number  of  important  conclusions  regarding  the 
accuracy  and  the  degree  of  conservatism  of  the  procedures  for  analyzing  rapid  drawdown 
slope  stability.  From  these  results,-  the  best  procedure  should  have  the  following 
characteristics: 

a.  The  procedure  should  employ  a  measure  of  shear  strength  and  pore  pressure 
that  reflects  the  actual  properties  of  the  soils  in  the  slopes. 

b.  The  procedures  shoul d ^employ  a  soundly  based  technique  for  relating  shear 
strength  to  effective  consolidation  pressures. 

c.  The  procedure  should  not  include  strength  components  due  to  negative  pore 
water  pressures  to  avoid  overestimating  factors  of  safety  for  conditions 
where  partial  drainage  occurs  during  drawdown. 

A  method  with  these  characteristics  can  be  conceived  as  a  combination  of  methods 
currently  in  use.  It  could  be  described  as  being  a  modification  of  the  method  currently 
used  by  the  US  Army  Corps  of  Engineers,  the  modification  being  an  adoption  of  Lowe  and 
Karafiath's  procedure  for  representing  the  undrained  strength  of  the  soil  rather  than 
using  the  R-envelope  as  is  now  done. 
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PREFACE 


Although  considerable  advances  have  been  made  during  the  past  several 
years  in  the  understanding  of  shear  strength  and  computation  procedures  for 
estimating  slope  stability,  the  procedures  for  estimating  slope  stability 
during  sudden  drawdown  have  remained  relatively  poorly  understood.  The  lack 
of  understanding  of  procedures  for  slope  stability  computations  for  sudden 
drawdown  became  especially  evident  to  one  of  the  writers  (Wright)  during 
recent  efforts  to  develop  a  computer  program  for  slope  stability  for  the  US 
Army  Corps  of  Engineers  (CE).  Rather  than  settle  for  adoption  of  a  procedure 
in  the  computer  program  somewhat  arbitrarily,  without  understanding  the 
potential  differences  among  procedures,  the  US  Army  Engineer  Waterways  Experi¬ 
ment  Station  (WES)  contracted  (Contract  Nos.  DACW39-86-M-0982  and  DACW39-86-M- 
1731)  with  the  writers  to  prepare  a  state-of-the-art  paper  on  the  subject  of 
stability  computations  for  sudden  drawdown.  The  results  of  that  effort  are 
presented  in  this  paper. 

An  initial  draft  of  this  paper  was  first  prepared  by  Dr.  Stephen  G.  Wright 
and  reviewed  by  Dr.  J.  M.  Duncan.  A  separate  draft  set  of  comments,  conclu¬ 
sions,  and  recommendations  was  subsequently  prepared  by  Dr.  Duncan.  The  two 
drafts  were  then  combined  in  the  final  form  presented  herein.  Various  CE 
personnel  and  consultants  also  reviewed  the  original  draft  of  this  paper  and 
offered  many  useful  suggestions  and  comments.  Their  contributions  are  grate¬ 
fully  acknowledged  and,  to  the  extent  feasible,  their  suggestions  have  been 
incorporated  into  this  paper. 

The  writers  express  their  appreciation  to  the  members  of  the  Computer 
Applications  in  Geotechnical  Engineering  (CAGE)  Committee  of  the  CE  for  their 
many  comments  and  useful  suggestions.  Special  appreciation  is  extended  to 
Messrs.  Dale  Munger,  Office,  Chief  of  Engineers  (OCE);  and  Earl  Edris,  Engi¬ 
neering  Group,  Soil  Mechanics  Division,  Geotechnical  Laboratory  (GL),  WES; 
Professors  Tom  Wolff,  Michigan  State  University;  and  Milton  Harr,  Purdue 
University,  for  their  many  detailed  comments  on  the  first  draft  of  this  report. 
The  interest  and  support  of  the  CE  provided  both  the  catalyst  and  means  that 
made  this  study  possible. 

The  report  was  funded  by  the  CAGE  project  sponsored  by  OCT. 

Mr.  William  E.  Strohm,  Jr.,  Engineering  Geology  and  Rock  Mechanics  Division, 

GL,  was  Principal  Investigator  of  the  CAGE  project.  At  WES,  the  work  was 


under  the  management  of  the  Soil  Mechanics  Division  (SMD),  GL. 

The  Contract  Monitor  was  Mr.  Earl  Edris.  Mr.  Clifford  L.  McAnear  was 
Chief,  SMD;  and  Dr.  William  F.  Marcuson  III  was  Chief,  GL,  during  the  prepara 
tion  of  this  report. 

COL  Dwayne  G.  Lee,  CE,  is  Commander  and  Director  of  WES.  Dr.  Robert  W. 
Whalin  is  Technical  Director. 
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-  Skempton's  pore  pressure  coefficient 

-  Skempton's  pore  pressure  coefficient  (=  AB) 
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-  Skempton's  pore  pressure  coefficient 
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-  cohesion  intercept  for  effective  stress  envelope  (from  either  R 
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circle  of  stress  representing  stresses  on  the  failure  plane 

-  void  ratio 

-  acceleration  due  to  gravity 
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procedure  for  estimating  effects  of  anisotropic  consolidation  on 
undrained  shear  strengths  using  results  of  tests  with  isotropic 
consol idation 


SECTION  1 


INTRODUCTION 

(Wright) 

"Sudden"  or  "rapid"  drawdown  is  considered  to  occur  when  reservoir 
or  other  adjacent  water  levels  are  lowered  at  such  a  rate  that  little  or 
no  drainage  occurs  in  an  earth  slope  while  the  water  level  is  being 
lowered.  Slope  stability  analyses  are  routinely  performed  to  calculate 
the  factor  of  safety  for  earth  slopes  subjected  to  this  condition. 
Frequently  the  design  of  such  earth  slopes  is  governed  by  the  sudden 
drawdown  condition  rather  than  by  other  conditions,  such  as  those 
immediately  after  construction  and  with  steady  state  seepage. 

Several  different  procedures  exist  and  are  currently  used  by 
designers  to  compute  slope  stability  for  the  case  of  sudden  drawdown. 
These  various  procedures  are  based  on  fundamentally  different  approaches 
and  are  known  to  produce  different  results  in  at  least  some  instances. 
The  procedures,  the  assumptions  employed  and  the  numerical  results 
obtained  using  these  procedures  are  the  subject  of  this  paper. 

BACKGROUND  AND  GENERAL  ASSUMPTIONS 

Most  stability  computations  for  sudden  drawdown  are  based  on 
several  common  assumptions.  First,  most  procedures  are  based  on  the 
assumption  that  the  drawdown  is  instantaneous.  As  the  soil  drains,  the 
shear  strength  will  change  and,  eventually,  will  approach  values 
corresponding  to  the  fully  "drained"  or  steady  seepage  condition.  By 
considering  the  case  of  instantaneous  drawdown  as  well  as  the  steady 
seepage  condition  it  is  assumed  that  the  worst  condition  will  be 
detected  and  can  be  used  to  judge  the  design.  Intermediate  conditions 


between  the  undrained  and  fully  drained  condition  are  usually  ignored. 
Only  the  condition  of  instantaneous  drawdown  with  no  drainage  will  be 
considered  in  this  paper.  Stability  computations  for  the  steady  seepage 
condition  are  routinely  performed  and  do  not  need  to  be  covered  here. 
Although  partial  drainage  during  drawdown  is  not  covered  in  this  paper, 
the  techniques  usually  used  are  based  on  solutions  to  governing 
equations  for  transient  flow,  which  are  discussed  in  Section  2  of  this 
paper . 

Procedures  for  sudden  drawdown  stability  computations  are  also  all 
based  on  the  assumption  that  the  soil  is  saturated.  This  assumption  is 
consistent  with  the  usual  assumption  that  sudden  drawdown  will  follow  a 
period  of  relatively  steady  reservoir  or  river  levels  and  the  soil  will 
have  ample  opportunity  to  become  saturated.  Although  there  are  probably 
many  cases  where  the  soil  does  not  become  saturated,  especially  when 
drawdown  follows  the  first  filling  of  a  reservoir  or  an  unprecedented 
flood,  the  possibility  still  exists  that  the  soil  will  at  some  time 
become  saturated.  Shear  strengths  will  be  at  their  minimum  values  when 
such  saturation  occurs,  and,  accordingly  the  assumption  of  saturated 
conditions  appears  to  be  appropriate  for  design.  However,  the  fact  that 
the  soil  may  not  be  saturated  at  the  time  of  some  failures  may  need  to 
be  recognized  when  actual  failures  an-  e«ammed 

Factors  of  safety  calculated  by  the  various  procedures  for  sudden 
drawdown  will  also  be  influenced  by  the  cart'u  ar  mt  equ’',b*ium 
slope  stability  analysis  procedure  emp’oyed  'or  e«ampie.  'he  mriji  , 
Method  of  Slices  will  typically  underest ’mate  the  '  •  >*  sa'et,  wh  -  ' 
the  Corps  of  Engineers'  Modified  Swedsn  procedure  .  .e-est  ^a'e  ’ 
factor  of  safety  in  comparison  to  "e  .a  ues  D'amed  .s  •',)  e  :  . • es 


based  on  more  rational  assumptions.  However,  most  of  the  differences 
among  factors  of  safety  caused  by  the  procedures  used  to  compute  the 
factor  of  safety  can  be  eliminated  by  use  only  of  procedures  wMch 
satisfy  complete  static  equilibrium.  Such  complete  equilibrium 
procedures  are  readily  available  and  have  been  implemented  in  many 
computer  programs  for  slope  stability  analysis,  including  the  program 
UTEXAS2  which  is  currently  being  implemented  by  the  Corps  of  Engineers. 
Accordingly,  the  subject  of  assumptions  and  differences  among  various 
limit  equilibrium  procedures  is  of  little  concern  and  is  ignored  for  the 
balance  o*  this  paper.  It  will  be  assumed  that  any  of  the  various 
procedures  discussed  for  stability  analysis  for  sudden  drawdown  would  be 
used  with  a  limit  equilibrium  procedure  which  satisfies  complete  static 
equ i I i b  r i urn . 

ORGANIZATION  OF  REPORT 

Procedures  for  stability  calculations  for  sudden  drawdown  can  De 
separated  into  two  general  groups:  (1)  Effective  stress  procedures,  and 
(?)  Two-stage  procedures.  Effective  stress  procedures  are  based  on 
defining  shear  strengths  m  terms  of  effective  stresses  and  require  that 
pore  water  pressures  'mmediately  after  drawdown  be  determined.  Effective 
st-ess  procedures  are  described  and  discussed  in  Section  2.  Two-stage 
ana  /s' s  procedures  require  mat  two  separate  sets  of  stabi'ity 
a  >  j  ;  at  ions  be  performed  the  *  mt  set  is  performed  for  ,:ond't'0’'s 
•  ,sT  before  drawdown  to  compute  the  effective  stresses,  wh’ch  a>e 


iSed  '  esf'-mate  j'’1'  a’"e<l  snea'  smengms  'c  subse  t  a 

'if  j'd  'j-'-ed  snea*-  s*  •-e-gtr's  a»  e  *ne«  .sed  •  •  e 
•  jmp  it  at  S  *  :  ■  m{;  ,  r  e  ’  ne  •  1  »  qt  .  *  s  i  *  ef  ,  i  f 

'  is  owe’'  eg  ' wo  - sf  aue  i>  a  .  ■-  *-  :  •  ■  e  !  i  •  •• 


3.  One  of  the  two-stage  analysis  procedures  (Lowe  and  Karafiath's  , 
1960a)  as  well  as  an  alternative  effective  stress  analysis  procedure 
which  is  developed  and  presented  in  this  paper  take  into  account  the 
fact  that  the  soil  is  consolidated  anisotropical ly  under  unequal 
principal  stresses  (o^,  o^)  before  drawdown.  Procedures  for  estimating 
the  effects  of  anisotropic  consolidation  from  the  results  of  laboratory 
tests  employing  isotropic  consol idation  are  discussed  in  Section  4.  The 
alternative  effective  stress  procedure  considered  in  this  paper  is 
presented  in  Section  5.  Results  of  several-  series  of  comparative 
calculations  for  the  factor  of  safety  are  presented  in  Section  6. 
Finally,  a  summary  and  recommendations  are  presented  in  Section  7. 


SECTION  2 


EFFECTIVE  STRESS  (ONE-STEP)  ANALYSIS  PROCEDURES 

(Wright) 

Effective  stress  analysis  procedures  are  based  on  expressing  the 
shear  strength  in  terms  of  effective  stresses.  Accordingly,  analyses  are 
based  on  effective  stresses  and  pore  water  pressures  must  be  estimated 
for  the  condition  immediately  after  drawdown.  The  only  difference 
between  an  effective  stress  analysis  for  conditions  after  drawdown  and 
an  effectives  stress  analysis  for  long-term,  steady  seepage  conditions 
is  in  the  pore  water  pressures  and  surface  loads  used  in  the 
computations.  In  the  case  of  sudden  drawdown  analyses  the  pore  water 
pressures  should  represent  the  pressures  expected  immediately  after 
drawdown;  for  long-term,  steady-state  seepage  analyses  the  pore  water 
pressures  represent  those  associated  with  steady  seepage. 

SHEAR  STRENGTH  PARAMETERS 

Effective  stress  shear  strength  parameters  may  be  determined  from 

either  consol idated-undrained  tests  with  pore  water  pressure 

measurements  or  from  consol i dated-drai ned  tests.  In  the  case  of 

consol idated-undrained  tests  two  candidate  criteria  exist  for  defining 

"failure"  and  determining  the  effective  stress  shear  strength 

parameters:  Failure  may  be  defined  to  be  either  at  the  point  of  maximum 

principal  stress  difference,  (a.  -  a.,)  ,  or  at  the  point  where  the 

l  j  n?d  x 

effective  stress  path  becomes  tangent  to  the  failure  envelope  (Fig. 


2.1)1.  jn  case  rapid  drawdown  the  loading  is  considered  to  be 
undrained.  Consequently,  it  is  fundamentally  more  correct  to  define 
failure  and  the  corresponding  effective  stress  shear  strength  parameters 
at  the  point  of  peak  principal  stress  difference;  the  peak  principal 
stress  difference  represents  the  point  of  peak  load  (shear  resistance) 
for  undrained  loading.  However,  differences  between  the  effective  stress 
shear  strength  parameters  at  the  point  of  peak  principal  stress 
difference  and  stress  path  tangency  are  likely  to  be  relatively  small 
for  compacted  soils  and  soils  with  low  sensitivities.  Such  soils  are  the 
most  likely  types  of  fine-grained  soils  encountered  in  earth  dams  and 
riverbank  slopes  where  drawdown  is  significant. 

Consol i dated-drai ned  shear  test  procedures  may  also  be  used  to 

determine  the  effective  stress  shear  strength  parameters.  Fundamentally, 

the  effective  stress  shear  strength  parameters  determined  in 

consol idated-drained  shear  tests  should  be  different  from  those 

applicable  to  undrained  loading  associated  with  rapid  drawdown  for  two 

reasons:  First,  the  effective  stress  strength  parameters  determined  in 

consol i dated-drai ned  shear  tests  correspond  to  points  of  stress  path 

tangency  on  the  effective  stress  failure  envelope,  (o./oJ  ;  the 

i.  o  m3  x . 

effective  stress  strength  parameters  applicable  at  failure  for  undrained 

^If  no  cohesion  (c)  intercept  exists,  the  point  where  the  effective 
stress  path  becomes  tangent  to  the  failure  envelope  coincides  with  the 
point  of  maximum  effective  principal  stress  ratio,  (c./d ..)  .  Thus, 

"maximum  effective  principal  stress  ratio"  and  "stress  p3tWa£angency" 
are  synonymous  failure  criteria.  If  a  cohesion  intercept  exists,  the  two 
failure  criteria  are  not  identical,  and  it  is  fundamentally  more  correct 
to  use  the  point  of  stress  path  tangency,  although  the  point  of 
effective  principal  stress  ratio  is  often  taken  to  be  the  point  of 
failure  for  convenience. 


loading  are  those  corresponding  to  the  peak  principal  stress  difference. 


(o,  -  o,)  ,v  .  Secondly,  in  consol i dated-drai ned  shear  tests  a  portion 
of  the  load  applied  to  cause  failure  goes  into  producing  volume  change 
and  the  strength  parameters  measured  reflect  this  additional  energy 
required  to  produce  the  volume  changes  (Bishop  and  Eldin,  1953;  Bishop, 
1954a;,  Rowe,  1962;  Rowe,  Barden  and  Lee,  1964).  For  undrained  loading 
as  in  the  cases  of  sudden  drawdown  and  consol i dated-undrai ned  shear 
tests  in  the  laboratory,  no  energy  is  used  to  produce  volume  changes 
and,  accordingly,  any  measured  shear  strength  parameters  do  not  reflect 
such  added  energy.  Consequently,  the  effective  stress  shear  strength 
parameters  determined  from  consol idated-drained  and  consol i dated- 
undrained  shear  tests  will  not  be  identical. 

Effective  stress  shear  strength  parameters  determined  for  a 
particular  soil  will  clearly  differ  depending  on  the  test  procedures 
(drained  vs.  undrained)  employed  and  the  manner  in  which  -.ilure  is 
selected.  In  addition  to  the  two  failure  criteria  discussed  earlier,  a 
limiting  strain,  e.  g.  15  percent,  is  sometimes  selected  as  the  point 
defining  failure  in  a  laboratory  shear  test.  Thus,  a  third  potential 
"failure  criterion"  is  introduced  and  further  differences  among  shear 
strength  parameters  may  exist.  Additional  differences  may  also  be  caused 
by  the  type  of  laboratory  shear  device  employed  (triaxial  shear,  direct 
shear,  direct  simple  shear  and  plane  strain).  In  most  cases  all  of  the 
factors  discussed  pertaining  to  laboratory  measurements  of  shear 
strength  will  probably  be  overshadowed  by  scatter  in  data  and  the 
variability  of  most  soils  in  the  field.  However,  all  of  the  factors  will 
contribute  to  at  least  some  extent  to  the  computed  results  and  may 
influence  conclusions  drawn  in  any  particular  instance. 
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ESTIMATION  OF  PORE  WATER  PRESSURES 

All  of  the  effective  stress  analysis  procedures  require  that  pore 

water  pressures  at  the  end  of  drawdown  be  estimated.  At  least  four 
fundamentally  different  approaches  have  been  used  to  estimate  the  pore 
water  pressures.  The  first  procedure  was  first  suggested  by  Bishop 
(1954b)  and  later  employed  by  Morgenstern  (1963)  using  Skempton's  (1954) 
pore  pressure  coefficients.  The  procedure  involves  estimating  the  pore 
water  pressures  by  first  estimating  what  the  pore  pressures  would  be 
before  drawdown  and  then  estimating  the  changes  in  pore  water  pressure 
which  would  occur  during  drawdown.  The  changes  in  pore  water  pressure 
are  estimated  by  estimating  the  changes  in  total  stress  and  relating  the 
changes  in  total  stress  to  changes  in  pore  water  pressure  using 
Skempton's  pore  pressure  coefficients.  The  second  approach  to  estimating 
pore  water  pressures  involves  solving  either  analytically  or  numerically 
the  basic  governing  differential  equation  for  the  transient  (non-steady) 
flow  of  water  associated  with  the  drawdown  of  the  reservoir.  Procedures 
based  on  the  second  approach  generally  employ  solutions  developed  for 
either  soil  consolidation  or  groundwater  flow.  The  third  approach  is 
actually  a  simplification  of  the  second  approach,  consisting  of 
representing  the  flow  as  steady-state  flow  for  selected  instants  of  time 
following  drawdown.  The  fourth  approach  is  the  least  well-defined 
approach  and  essentially  encompasses  those  approaches  in  which  the  pore 
water  pressures  are  estimated  based  on  past  experience  and  judgement. 

BISHOP-MORGENSTERN  PROCEDURE 

Bishop's  (1954b)  and  Morgenstern ' s  (1963)  procedure  for  estimating 
pore  water  pressures  after  drawdown  is  based  on  Skempton's  pore  pressure 
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coefficient,  B,  and  estimated  changes  in  total  stress  during  drawdown. 
The  pore  pressure  coefficient  B  is  defined  as, 


where  Au  is  the  change  in  pore  water  pressure  and  Ao^  is  the  major 
principal  stress  change.  Thus,  the  change  in  pore  water  pressure  during 
drawdown  can  be  expressed  as, 

au  =  Bao^  2.2 

The  smaller  the  value  of  B,  the  smaller  the  changes  in  pore  water 
pressure  during  drawdown.  Because  the  changes  are  expected  to  represent 
decreases  in  pressure  the  smaller  the  decrease  in  pressure  during 
drawdown,  the  higher  will  be  the  final  values  of  pore  water  pressure. 
Higher  values  of  pressure  correspond  directly  to  lower  factors  of 
safety.  Thus,  the  factor  of  safety  will  decrease  as  B  becomes  smaller, 
all  other  factors  remaining  the  same.  Bishop  and  Morgenstern  each 
considered  probable  values  for  B  and  suggested  that  a  value  of  unity 
represented  a  lower-bound  and,  thus,  conservative  value.  Consequently, 
and  based  on  Eq.  2.2,  the  change  in  pore  water  pressure  becomes  equal  to 
the  major  principal  stress  change,  Ao^  (B  =  1.0  in  Eq.  2.2). 

Bishop  and  Morgenstern  also  assumed  that  the  major  principal  stress 
change,  Aop  was  equal  to  the  change  in  total  vertical  stress  on  the 
slope  directly  above  any  particular  point  of  interest.  Thus, 

Ao,  =  Ah  y  2.3 

1  w'w 

where  Ah^  is  the  change  in  height  of  water  above  the  surface  of  the 
slope  directly  above  the  point  of  interest  (Fig.  2.2)  and  yw  is  the  unit 
weight  of  water.  Combining  Eqs.  2.2  and  2.3  and  setting  B  equal  to  unity 
then  gives, 

AU  =  Ah  Y  2.4 

w'w 
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Bishop  and  Morgenstern  also  assumed  that  the  head  loss  in  the  slope 
prior  to  drawdown  is  small  and  can  be  ignored.  The  consequence  of  this 
assumption  combined  with  the  previous  assumptions  and  Eq.  2.4  is  that 
the  pore  pressures  after  drawdown  can  be  expressed  by  a  simple 
piezometric  surface  as  illustrated  in  Fig.  2.3.  The  piezometric  surface 
is  coincident  with  the  surface  of  the  slope  above  the  level  of  the 
drawdown  and  is  coincident  with  the  reservoir  surface  at  and  below  the 
level  of  the  drawdown.  Pore  water  pressures  are  simply  the  product  of 
the  unit  weight  of  water  and  the  depth  below  the  final  piezometric 
surface. 

Bishop's  and  Morgenstern 1 s  procedure  involves  the  following  four 
assumptions: 

(1)  The  head  loss  due  to  seepage  in  the  slope  before  drawdown  is 
negligible. 

(2)  Skempton's  (1954)  equation  and  pore  pressure  coefficients  can 
be  used  to  compute  changes  in  pore  water  pressures. 

(3)  The  pore  pressure  coefficient  B  is  equal  to  unity. 

(4)  The  major  principal  stress  change,  Ao^,  is  equal  to  the  change 
in  vertical  stress  on  the  surface  of  the  slope  directly  above  any 
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point  of  interest. 

The  assumption  that  the  head  loss  in  the  slope  is  negligible  appears  to 
be  reasonable  inasmuch  as  slides  due  to  sudden  drawdown  are  typically 
shallow.  In  addition  the  assumption  of  negligible  head  loss  is  made 
primarily  for  simplicity;  the  concepts  employed  in  Bishop's  and 
Morgenstern ' s  procedure  could  be  easily  extended  to  include  head  losses 
in  the  slope  before  drawdown  if  they  were  important.  Bishop's  and 
Morgenstern ' s  procedure  is  also  based  on  the  fundamental  assumption  that 


Skempton's  (1954)  equation  for  pore  water  pressures  can  De  app’ied  to 
predicting  pore  water  pressures  in  earth  slopes  due  to  sudden  d^awdowr 
For  the  present  this  will  be  considered  valid  and  tne  use  of  Sxempton  s 
equation  is  discussed  in  further  detail  in  Section  5  o*  tn’s  paper. 
However,  the  remaining  two  assumptions  listed  above  requ^'e  some  fjrtner 
discussion  at  this  point. 

The  assumption  that  the  pore  pressure  coefficient,  B,  is  equal  to 
unity  is  best  examined  by  considering  Skempton's  original  expression  for 
the  changes  in  pore  water  pressure  due  to  undrained  loading.  Skempton 
expressed  the  change  in  pore  pressure  during  undrained  loading  as, 

Au  =  BAo^  +  AB(Ao^  -  Ao^  )  2.5 

where,  Au  is  the  change  in  pore  water  pressure,  Ao^  and  a a^,  are  the 
major  and  minor  principal  stress  changes,  respectively,  and  A  and  B  are 
pore  pressure  coefficients.  The  coefficient  B  represents  tendencies  for 
pore  water  pressures  to  change  due  to  a  change  in  all-around  (isotropic) 
pressure;  while  the  coefficient  A  represents  tendencies  for  pore  water 
pressures  to  change  due  to  a  change  in  shear  stress  (principal  stress 
difference).  By  dividing  both  sides  of  Eq.  2.5  through  by  the  major 
principal  stress  change  (Ao^)  and  rearranging  terms  the  following 
equation  can  be  written: 

^  =  B[1  -  (1  -  A) ( 1  -  — )]  2.6 

AOj  AOj 

The  quantity  on  the  left-hand  side  of  this  equation  is  equal  to  the  pore 
pressure  coefficient,  B.  In  the  case  of  most  saturated  soils  the  value 
of  B  will  be  essentially  unity  and,  thus. 


Accordingly,  the  pore  pressure  coefficient  B  can  be  related  to  the  pore 


pressure  coefficient  A  and  the  principal  stress  changes,  Ao^  and  At 

this  point  the  "principal  stress  changes,"  Ao^  and  which  have  been 

used  in  Eqs.  2.1  through  2.7  should  be  distinguished  from  "changes  in 
principal  stress."  The  major  principal  stress  change,  Ao,,  in  these 
equations  represents  the  largest  algebraic  (most  positive)  change  in 
stress.  In  the  case  of  drawdown,  the  changes  in  normal  stress  are  likely 
to  all  be  negative  in  value.  Thus,  the  least  negative  value  is  the  major 
principal  stress  change,  while  the  most  negative  value  is  the  minor 
principal  stress  change  (Fig.  2.4).  Accordingly,  the  absolute  value  of 
the  major  principal  stress  change,  Ao ^ ,  is  likely  to  be  smaller  than  the 
absolute  value  of  the  minor  principal  stress  change,  A o^.  Referring  to 
Eq.  2.7  it  is  then  likely  that  both  Ao^  and  Ao^  will  be  negative  and  the 
absolute  value  of  Ao^  will  be  larger  than  the  absolute  value  of  Ao^. 
Thus,  the  ratio,  Ao^/Ao^,  will  be  greater  than  unity.  Consequently,  the 
value  B  given  by  Eq.  2.7  will  always  be  greater  than  unity  for  all 
values  of  A  less  than  or  equal  to  unity.  Most  compacted  soils  as  well  as 
natural  soils  in  slopes  likely  co  be  subjected  to  sudden  drawdown  would 
be  expected  to  produce  values  of  A  which  are  less  than  unity. 
Accordingly,  Bishop’s  and  Morgenstern 1 s  assumption  that  a  conservative 
lower-bound  value  for  B  is  unity  appears  to  be  valid.  The  degree  of 
conservatism  associated  with  this  assumption  may  be  seen  py  examining 
values  of  B  computed  for  various  values  of  A  and  ratio  of  principal 
stress  change  (Ao^/Ao^)  as  shown  in  Fig.  2.5.  Referring  to  this  figure 
it  can  be  seen  that  the  values  of  B  are  at  least  unity  and  in  many  cases 
may  be  several  times  greater  than  unity.  For  example,  if  the  ratio  of 
principal  stress  changes  is  2.0  and  A  is  0.2,  a  value  which  is  certainly 


reasonable  for  a  compacted  clay  at  low  confining  pressures,  the  value  of 
B  is  1.80.  This  value  (1.80)  is  approaching  twice  the  value  assumed  d 
Morgenstern.  Consequently,  the  change  in  pore  water  pressure  during 
drawdown  would  be  nearly  twice  the  va'ue  predicted  based  on  Bishop's  ana 
Morgenstern ‘ s  assumed  values. 

The  final  assumption  to  oe  considered  in  Bishop's  and  Morgenstern ' s 

procedure  is  that  the  major  principal  stress  change  is  equal  to  the 

change  in  vertical  stress  directly  above  any  point  of  interest.  At  the 

face  of  the  slope  tne  change  in  stress  normal  (perpendicular)  to  the 

face  of  the  slope  will  be  equal  to  the  change  in  water  pressure  (or 

vertical  stress)  on  the  face  of  the  slope.  The  change  in  normal  stress 

on  the  face  of  the  slope  will  be  negative  but  greater  in  magnitude  than 

the  change  in  normal  stress  on  a  plane  perpendicular  to  the  slope  face. 

Thus,  at  the  face  of  the  slope  the  change  in  vertical  stress.  *h  v  will 

w  w 

actually  be  equal  to  the  minor  principal  stress  change.  while  the 

major  principal  stress  change.  An,,  will  be  smaller  in  absolute  value 
than  the  change  in  vertical  stress  (Fig.  2.6).  The  consequence  of  this 
will  oe  that  the  actual  change  in  pore  pressure  computed  from  Eq.  2.2 
would  be  smaller  than  what  is  computed  by  Morgenstern 1 s  assumption, 
provided  that  all  other  factors  remain  the  same,  including  the 
assumption  that  B  is  unity. 

Based  on  the  above  discussion  it  appears  that  Bishop's  and 
Morgenstern ' s  procedure  may  underestimate  the  actual  va'ue  of  B  while  at 
tne  same  time  the  major  principal  stress  change,  'o,.  ' s  overes t ; mated . 
Consequent 1 y .  it  is  entirely  possible  that  two  errors  may  compensate 
each  other  in  tne  procedure 


TRANSIENT  FLOW  SOLUTIONS 


During  the  past  15  years  the  writer  has  gained  the  impression  that 
in  a  number  of  instances  pore  water  pressures  due  to  drawdown  have  been 
estimated  from  theoretical  solutions  based  on  the  equations  of  transient 
flow.  Although  the  details  of  the  solutions  are  not  available,  the 
solutions  appear  to  have  been  based  generally  on  numerical  solutions  of 
the  equations  for  consolidation  or  groundwater  flow.  In  this  section  the 
governing  equations  for  transient  flow  are  presented  and  the  assumptions 
employed  in  the  classical  versions  of  these  equations  for  soil 
consolidation  and  groundwater  flow  are  examined.  Requirements  for 
correct  modeling  in  the  case  of  sudden  drawdown  are  also  presented. 

The  general  governing  partial  differential  equation  which  must  be 
satisfied  for  two-dimensional  transient  (time-dependent)  flow  in  porous 
media  is. 
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where,  h  is  total  head,  n  is  porosity,  p  is  the  mass  density  of  the 
fluid  (y^/g) ,  and  is  the  degree  of  saturation  (e.  g.  Freeze  and 
Cherry,  1979).  The  partial  derivatives  on  the  left  are  taken  with 
respect  to  the  spatial  coordinates  (x  and  y);  the  partial  derivatives  on 
the  right-hand  side  are  taken  with  respect  to  time  (t).  In  the  case  of 


saturated  flow  (S^  =  100%  -  aS^/at  =  0)  and  an  incompressible  fluid 


(ap/at  =  0),  Eq.  2.8  reduces  to 


pk 


2.9 


In  order  to  solve  Eq.  2.9  either  the  head  (h)  must  be  related  to  the 
porosity  (n),  or  the  porosity  must  be  related  to  the  head.  The  most 
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widely  employed  procedure  for  relating  porosity  to  head  follows  the  one 
employed  in  the  classical  consol idation  theory  first  suggested  by 
Terzaghi.  This  approach  can  be  derived  by  first  expressing  the  partial 
derivative  of  porosity  with  respect  to  time  as, 


dn  _  1  3e  9o  3u  3h  „ 

at  ”  1  +  e  ao  au  ah  at 
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where,  3e/3o  represents  the  change  in  void  ratio  with  respect  to 
effective  stress,  3o/3u  represents  the  change  in  effective  stress  with 
respect  to  pore  pressure,  3u/3h  represents  the  change  in  pore  pressure 
with  respect  to  total  head,  and  3h/3t  represents  the  change  in  total 
head  with  respect  to  time.  The  change  in  void  ratio  with  respect  to 
change  in  effective  stress  can  be  expressed  by  the  coefficient  of 
compressibility,  ay,  where. 
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2.11 


The  change  in  effective  stress  is  assumed  to  be  equal  in  magnitude  and 
opposite  in  sign  to  the  change  in  pore  water  pressure.  Thus, 


The  change  in  pore  pressure  with  respect  to  head  can  is  expressed  as, 


3u 
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2.13 


Substituting  Eqs.  2.11,  2.12  and  2.13  into  2.10  gives, 

3n  1  3h 
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o 

Finally,  substituting  Eq.  2.14  into  Eq.  2.5  and  rearranging 

3^h  +  3^h  _  av  \  ah 
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terms  gives. 


2.15 


In  the  case  of  one-dimensional  flow  Eq.  2.15  reduces  to. 
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which  is  Terzaghi's  classical  equation  for  one-dimensional  equation 
expressed  in  terms  of  total  head.  Equation  2.16  can  also  be  expressed  in 
terms  of  excess  pore  pressure,  u,  as, 
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An  alternate  form  of  Eq.  2.15  is  often  used  in  modeling  groundwater 
problems.  The  equation  for  two-dimensional  transient  flow  is  written  as, 
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"specific  storage." 


2.18 

The  specific  storage  is  expressed 


Ss  =  Yw(a  +  np)  2.19 
where,  a  and  3  are  the  compressibilities  of  the  soil  skeleton  and  water, 
respect i vc 1 y .  That  is, 
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and , 
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Equations  2.15  and  2.18  are  fundamentally  identical  except  that  Eq.  2.18 
also  takes  into  account  the  compressibility  of  water.  Solutions  to  the 
two  equations  should  be  essentially  identical. 

Discussion  of  Assumptions 


purposes.  Equations  2.15  and  2.18  are  also  based  on  a  number  of  other 
important  assumptions: 

1.  The  soil  is  saturated  and  remains  saturated  during  flow. 

2.  Darcy's  Law  is  valid  (Assumed  in  deriving  Eq.  2.8). 

3.  The  void  ratio  is  related  to  a  single  effective  normal  stress 
component  (a),  regardless  of  the  state  of  stress,  including  shear 
stresses . 

4.  The  change  in  effective  stress  is  equal  in  magnitude  to  the 
change  in  pore  water  pressure,  regardless  of  changes  in  total 
stress . 

The  assumption  that  the  fluid  is  incompressible  and  Darcy's  Law  is  valid 
seems  acceptable.  However,  the  remaining  three  assumptions  listed  may 
not  be  valid  and  are  discussed  in  further  detail  below. 

Soil  is  saturated.  In  the  case  of  drawdown  the  soil  may  become 
partly  saturated  above  the  phreatic  surface  as  drainage  occurs. 
Accordingly,  the  third  term,  npSS^/at,  in  the  governing  equation  (Eq. 
2.8)  for  transient  flow  must  be  retained.  This  requires  that  a 
relationship  be  established  between  the  degree  of  saturation  of  the  soil 
and  the  total  head  in  the  soil.  Although  such  relationships  can  be 
determined,  they  require  that  additional  soil  parameters  be  evaluated. 
Such  relationships  between  degree  of  saturation  and  total  head  do  not 
appear  to  have  been  included  in  past  uses  of  the  transient  flow 
equations  for  predicting  pore  water  pressures  during  sudden  drawdown. 

Single  effective  stress  component.  The  assumption  that  the  change 
in  volume  (void  ratio)  is  related  to  a  single  effective  normal  stress 
component  ignores  the  fact  that  shear  stresses  produce  volume  change 
and,  accordingly,  for  undrained  or  partially  drained  loading  shear 


stresses  will  affect  the  pore  water  pressures.  Correct  representation  of 
the  relationship  between  effective  stress  and  volume  change  requires 
that  a  complete  constitutive  (stress-strain)  model  be  used  and  that  all 
components  of  stress  be  considered.  For  two-dimensional  plane-strain 
conditions  three  stress  components  must  be  considered:  c^,  o  ,  and  t^. 
Constitutive  models  based  on  Hooke's  law  are  not  valid  for  the  present 
problem  because  they  do  not  properly  model  the  volume  change  behavior  of 
soil  during  shear  and,  accordingly,  they  will  not  predict  proper  changes 
in  pore  water  pressures.  During  the  past  decade  considerable  progress 
has  been  made  in  developing  constitutive  models  based  on  plasticity 
theory,  which  can  account  for  shear  induced  volume  changes.  Such  models 
show  considerable  promise;  however,  they  do  not  appear  to  have  been 
coupled  with  the  equations  for  transient  flow  and  applied  to  the  case  of 
reservoir  drawdown. 

Constant  total  stress  (A o  =  -  Au).  The  assumption  that  changes  in 
effective  stress  are  equal  to  the  changes  in  pore  water  pressure  is 
invalid  for  cases  where  there  is  a  significant  change  in  total  stress. 
For  example,  at  the  immediate  face  of  a  slope  which  is  subjected  to 
rapid  drawdown  the  total  stress  and  pore  water  pressure  change  by  the 
same  amount,  while  the  effective  stress  remains  unchanged.  Any  solution 
based  on  the  assumption  that  the  change  in  effective  stress  is  equal  to 
the  change  in  pore  pressure  is  clearly  not  valid  in  this  case.  Changes 
in  total  stress  and  their  effects  on  the  effective  stress  are  easily 
accounted  for  in  the  case  of  one-dimensional  problems:  For  example  in 
Terzaghi's  classical  problem  of  one-dimensional  consolidation  the  change 
in  total  stress  can  be  assumed  to  be  constant  over  any  vertical 
distance.  In  two  dimensions  the  changes  in  total  stress  will  vary  from 


point  to  point  in  the  slope  and  can  only  be  computed  from  a  complete 
stress  analysis.  By  coupling  the  stresses  with  the  equations  of 
transient  flow  as  discussed  in  the  previous  paragraph,  changes  in  total 
stress  should  automatically  be  taken  into  account;  however,  as  noted 
above  such  coupled  solutions  have  not  generally  been  used. 

SIMPLIFIED  PROCEDURES  BASED  ON  WATER  FLOW 

A  number  of  procedures  for  estimating  pore  pressures  have  been 
developed  based  on  various  simplified  hydraulic  flow  solutions.  These 
procedures  and  solutions  can  be  generally  grouped  into  two  categories: 
The  first  category  of  solutions  is  based  on  the  Dupuit  assumptions  and  a 
simplified  differential  equation  for  flow.  The  second  set  of  procedures 
is  based  on  a  solution  of  the  governing  partial  differential  equation 
for  two-dimensional  steady-state  flow.  These  procedures  and  the 
assumptions  employed  are  reviewed  below. 

Procedures  Based  on  Dupuit  Assumptions. 

Procedures  based  on  the  Dupuit  assumptions  employ  some  form  of  the 
following  partial  differential  equation; 

kh^  +  k(|£)2=  n  U  2.22 

ax2  3x  3t 

This  equation  is  based  on  the  Dupuit  assumptions  that  (1)  equi potenti al 
lines  are  vertical,  and  (2)  the  hydraulic  gradient  is  equal  to  the  slope 
of  the  phreatic  surface  (line  of  seepage)  and  does  not  vary  in  the 
vertical  direction.  In  addition  the  equation  is  based  on  a  number  of 
other  assumptions: 

1.  Darcy's  Law  is  valid. 

2.  The  soil  is  saturated  below  the  phreatic  surface  at  any  time. 


3.  The  flow  is  horizontal  in  the  x  direction. 


4.  The  soil  or  slope  in  which  flow  is  occurring  is  homogeneous  and 
underlaid  by  an  impervious  base.  Although  not  apparent  from  Eq. 
2.22,  "h"  represents  not  only  head,  but  also  the  elevation  of  the 
phreatic  surface  above  the  assumed  impervious  base.  Thus,  Eq.  2.22 
is  not  a  general  governing  differential  equation  for 
one-dimensional  flow  where  x  may  represent  a  direction  other  than 
horizontal . 

5.  The  soil  is  assumed  to  be  incompressible,  i.  e.  there  is  no 
volume  change  even  though  the  pore  pressures  (heads)  may  change 
with  time. 

6.  The  drainage  of  the  soil  in  the  zone  where  the  free  surface  is 
changing  with  time  is  represented  by  the  effective  porosity,  n 
ASTM  (1986)  defines  effective  porosity  as  “The  ratio  of:  (1)  the 
volume  of  the  voids  of  a  soil  or  rock  mass  that  can  be  drained  by 
gravity,  to  (2)  the  total  volume  of  the  mass."  In  the  case  of  Eq. 
2.22  as  the  phreatic  surface  drops  from  an  elevation,  h^ ,  to  a  new 
elevation,  h^,  the  volume  of  water  that  drains  from  the  soil  is 
equal  to  the  product  of  ng  and  the  total  volume  of  soil  between  the 
elevations  h^  and  h^. 

Equation  2.22  is  usually  simplified  by  ignoring  the  second  term  on 
the  left-hand  side  involving  products  of  derivatives,  which  gives, 
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Equation  2.23  is  also  sometimes  written  as, 
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Browzin  (1961)  developed  a  simple  analytical  solution  for  Eq.  2.23  while 
Desai  and  Sherman  (1971)  presented  a  numerical  solution.  Several  other 
published  solutions  (Brahma  and  Harr,  1963;  Newlin  and  Rossier,  1967) 
have  been  based  on  further  simplification  of  Eq.  2.23  achieved  by 
replacing  the  quantity  "h"  in  the  first  term  by  a  constant  quantity,  fi, 
representing  the  average  height  of  the  phreatic  surface  above  the 
impervious  base: 

2 
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The  various  solutions  to  Eqs.  2.23  or  2.25  have  been  directed  toward 
determining  the  drop  in  the  phreatic  surface  due  to  partial  drainage 
during  drawdown,  rather  than  toward  determining  the  pore  water  pressures 
at  the  end  of  a  sudden,  instantaneous  drawdown.  However,  regardless  of 
the  degree  of  drainage,  the  pore  water  pressures  are  defined  by  the 
location  of  the  phreatic  surface  and  the  assumption  of  vertical 
equipotential  lines.  For  an  instantaneous  drawdown  the  pore  pressures 
should  be  identical  to  those  determined  using  Bishop's  and  Morgenstern ' s 
procedures,  although  the  two  approaches  ( Bi shop-Morgenstern  and 
Dupuit-based  flow  solutions)  represent  totally  different  approaches 
fundamentally.  By  either  procedure  the  pore  pressures  are  simply 
calculated  by  taking  the  distance  from  the  surface  of  the  slope  (or 
reservoir  level  if  below  the  final  drawdown  level)  and  multiplying  the 
distance  by  the  unit  weight  of  water. 

Procedures  Based  on  Steady-State  Flow . 

A  number  of  procedures  have  employed  the  solutions  for  steady-state 
flow  for  the  purpose  of  calculating  pore  water  pressures  and  flow 
velocities  during  drawdown.  In  the  case  of  isotropic  materials  the 
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following  LaPlace  equation  is  the  governing  equation  which  must  be 
solved: 


2  2 
9x  9y 


2.26 


Cedergren  (1948)  appears  to  have  been  one  of  the  first  to  apply  this 
equation  to  the  case  of  drawdown.  He  constructed  flow  nets  to  solve  the 
equation  and  suggested  an  approach  for  accounting  for  partial  drainage 
and  the  drop  in  the  phreatic  surface  during  or  after  drawdown.  He  first 
constructed  a  flow  net  for  the  original  reservoir  level  and 
steady-seepage  (pre-drawdown)  condition.  Next,  using  the  original 
location  of  the  phreatic  surface  but  with  a  new  reservoir  level 
representing  drawdown  at  the  slope  face  he  constructed  a  new  low  net. 
For  the  new  flow  net  the  portion  of  the  upstream  slope  between  the 
original  and  new  reservoir  levels  as  well  as  the  phreatic  surface  were 
now  considered  zero  pressure,  drainage  (discharge)  boundaries.  Based  on 
the  velocities  of  flow  in  a  direction  perpendicular  to  the  phreatic 
surface  (obtained  from  the  new  flow  net),  Cedergren  estimated  the  amount 
of  water  which  would  have  flowed  away  from  the  original  phreatic  surface 
for  a  selected  increment  of  time.  This  allowed  a  new  position  for  the 
phreatic  surface  to  be  estimated  and  the  process  was  repeated.  This 
process  was  repeated  in  selected  increments  of  time  until  the  desired 


final  solution  in  time  was  obtained.  If  appropriate,  additional  changes 
in  the  reservoir  level  could  be  considered  at  any  stage  in  the  process. 
Desai  (1972,  1977)  used  the  same  approach  employed  by  Cedergren,  but 
employed  the  finite  element,  rather  than  graphical  (flow  net)  procedures 
to  solve  the  LaPlace  equation  ( Eq .  2.26). 

Although  procedures  based  on  the  governing  equation  for 


steady-state  flow  appear  to  have  been  used  in  a  number  of  instances  to 


obtain  pore  pressures  associated  with  reservoir  drawdown,  there  appears 
to  be  no  fundamental  basis  for  the  use  of  such  procedures.  As  shown 
previously  the  governing  equation  for  transient  flow  is, 
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Equation  2.8  is  clearly  not  the  same  as  Eq.  2.26  and  only  by  ignoring 
significant  terms  on  the  right-hand  side  of  Eq.  2.8  do  the  two  equations 
become  i denti cal  . 


EMPIRICAL  PROCEDURES  BASED-  ON  EXPERIENCE 
Pore  pressures  due  to  sudden  drawdown  appear  to  be  predicted 
sometimes  using  various  empirical  procedures  based  on  the  results  from 
either  model  tests  or  field  measurements .  For  example,  Nagy  (1967) 
presented  a  series  of  charts  showing  the  critical  rate  of  drawdown 
(required  to  cause  instability)  for  various  slope  angles,  effective 
stress  friction  angles  and  mean  soil  grain  diameter.  Nagy  considered  the 
mean  grain  diameter  to  suitably  characterize  the  soil  permeability  (and 
presumably  soil  volume  change  characteristics)  and  developed  his  chart 
based  on  the  results  of  a  series  of  model  experiments.  Although 
empirical  procedures  based  on  experimental  measurements  and  field 
experience  can  be  used  to  predict  pore  water  pressures  for  use  in  sudden 
drawdown  analyses,  the  lack  of  sufficient  data  and  the  complexity  of  the 
problem  appears  to  have  limited  the  use  of  this  approach  to  at  most  a 
few  isolated  cases. 


DISCUSSION 

None  of  the  procedures  examined  for  estimating  pore  pressures  due 
to  sudden  drawdown  can  be  justified  on  theoretical  bases;  only  the 
procedure  developed  by  Bishop  (1954b)  and  Morgenstern  (1963)  even 


attempts  to  account  for  changes  in  pore  water  pressure  which  occur  due 
to  induced  shear  stress  during  drawdown.  However,  even  Bishop's  ana 
Morgenstern 1 s  procedure  does  not  account  for  differences  among  various 
soils.  Instead  Bishop  and  Morgenstern  adopted  what  they  considered  to  be 
conservative  assumptions. 

The  only  way  in  which  any  of  the  effective  stress  procedures  can  be 
checked  is  by  comparison  of  results  obtained  using  the  procedures  with 
measured  results  for  actual  dams.  Such  comparisons  may  take  several 
forms.  Pore  pressures  may  be  measured  and  compared  directly  with  values 
which  are  predicted  by  the  various  theoretical  and  empirical  approaches. 
This  is  difficult  to  accomplish  and  due  to  point-to-point  variations  in 
pore  pressures  it  is  difficult  to  judge  the  overall  validity  of  various 
approaches.  Another  approach  is  to  examine  case  histories  where  slope 
failures  occurred  due  to  drawdown  and  to  compute  the  factor  of  safety  by 
various  procedures.  Wong,  Duncan  and  Seed  (1983)  computed  the  factors 
of  safety  for  two  dams  which  failed  by  sliding  due  to  sudden  drawdown: 
Pilarcitos  Dam  and  the  Walter  Bouldin  Dam.  They  found  that  the  factor  of 
safety  computed  by  Bishop's  and  Morgenstern ' s  procedures  were  1.15  and 
0.98,  respectively.  These  values  (1.15,  0.98)  are  sufficiently  close  to 
unity  for  these  dams  which  failed  to  suggest  that  the  pore  pressures 
predicted  were  reasonably  correct.  Wong,  Duncan  and  Seed  only  considered 
Bishop's  and  Morgenstern ' s  procedure  for  predicting  pore  water  pressures 
and,  thus,  no  assessment  was  made  of  other  effective  stress  procedures. 
Analyses  like  those  performed  by  Wong,  Duncan  and  Seed  are  very  useful; 
however,  the  number  of  wel 1 -documented  case  histories  of  actual  slope 
failures  due  to  sudden  drawdown  is  small  and  such  analyses  are  seldom 
performed . 


Much  of  the  existing  field  data  pertaining  to  embankment 
performance  during  sudden  drawdown  consists  of  measured  pore  pressures 
in  stable  embankments.  Winkley  (1982)  examined  measured  pore  pressure 
data  for  eight  embankment  dams  and  dikes  reported  in  the  literature  and 
summarized  in  Table  2.1.  He  computed  factors  of  safety  for  these  eight 
dams  using  measured  pore  pressure  information  and  pore  pressures  based 
on  Bishop's  and  Morgenstern ' s  procedure.  Because  sufficient  strength 
data  were  generally  not  available  for  the  dams.  Winx'iey  assumed  seve-'d 
sets  of  strength  values  and  computed  the  factor  of  safety  for  each  set; 
there  was  no  way  of  estaolishing  which  set  of  values  was  actually 
correct,  because  the  slopes  were  stable  and  the  factor  of  safety  was  not 
known.  The  factors  of  safety  are  summarized  in  Table  2.2  for  each  slope 
and  set  of  shear  strength  values  used.  Factors  of  safety  computed 
employing  pore  pressures  estimated  by  Bishop's  and  Morgenstern ' s 
procedure  are  plotted  versus  the  corresponding  values  computed  using  tne 
measured  pore  pressure  data  in  Fig.  2.7.  A  total  of  30  pairs  of  values 
are  plotted;  only  6  of  the  factors  of  safety  which  were  computed  using 
estimated  pore  water  pressures  exceed  those  based  on  the  measured  pore 
pressures.  The  mean  ratio  of  factors  of  safety  based  on 
estimated-to-measured  values  is  0.98,  with  a  standard  deviation  of  0.09. 
The  results  shown  in  Fig.  2.7  confirm  Bishop's  and  Morgenstern ' s 
hypothesis  that  their  procedure  is  conservative  and  are  in  substantia1 
agreement  with  the  results  of  Wong,  Duncan  and  Seed  for  embankments 
which  failed. 

Few  comparative  studies  appear  to  have  been  conducted  emp'cy'r'o 
pore  pressures  estimated  by  procedures  other  than  Bishop's  and 
Morgenstern '  s .  Some  insight  into  the  probable  relationship  r.etwee'--  'r'P 
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EMBANKMENT  DAMS  STUDIED  BY  WINKLEY  (1982) 


Amount 

Duration 

Height 

Upstream 

of 

of 

Name 

of  Dam 

Slope 

Drawdown 

Drawdown 

Ref  . 

(ft.  ) 

(ft.  ) 

( days ) 

Si r  Adam  Peak  Di ke 

52 

2.5:1 

25 

1 

Brazett 

(1961) 

Lonvarka  Dam 

537 

1.9:1 

53.5 

- 

Nonvei 1 1 er 

77.8 

- 

(1957) 

Gler.  Shi  ra  lower  Dam 

53 

3.0:1 

16 

2 

Patton  & 

24 

3 

Sempl e 

30 

4 

( 1961) 

Alcova  Dam 

239 

3.0:1 

126.6 

49 

Glover  et 

al  (1948) 

Cast i 1 etto  Dam 

246 

3.0:1 

194 

180 

Meyer- 

Peter 

(1964) 

Sardis  Dam 

97 

3.5:1 

15 

81 

Corps  of 
Engrs . 

(1967) 

Grenada  Dam 

80 

4.5:1 

16 

90 

Corps  of 
Engrs . 

(  1970) 

Goscheneral p 

345 

4.0:1 

283 

70 

Meyer- 

TABLE  2.2 


COMPARISON  OF  FACTORS  OF  SAFETY  COMPUTED  USING  MEASURED 
PORE  WATER  PRESSURES  AND  PORE  WATER  PRESSURES  ESTIMATED  BY 
BISHOP'S  AND  MORGENSTERN' S  PROCEDURE  (From  Winkley,  1982) 


(Assumed  Factor  of  Factor  of 

Strength  Safety  Safety  F 

Values)  (Meas.  (Bishop-  h  Bi sho 

Dam  A  Data)  Morgenstern)  F.. 

C0  /  »  /  Meas 


Sir  Adam  Peak  Dike 

11.8 

1.06 

1.36 

Sir  Adam  Peak  Dike 

5.9 

1.39 

1.45 

Sir  Adam  Peak  Dike 

3.9 

1.79 

1.75 

Sir  Adam  Peak  Dike 

3 

2.09 

2.04 

Lokvarka  -  Drawdown 

No.  1 

35.8 

1.29 

1.25 

Lokvarka  -  Drawdown 

No.  1 

17.9 

1.32 

1.29 

Lokvarka  -  Drawdown 

No.  1 

11.9 

1.36 

1.32 

Lokvarka  -  Drawdown 

No.  1 

8.9 

1.39 

1.36 

Lokvarka  -  Drawdown 

No.  2 

35.8 

1.28 

1.25 

Lokvarka  -  Drawdown 

No.  2 

17.9 

1.31 

1.26 

Lokvarka  -  Drawdown 

No.  2 

11.9 

1.35 

1.29 

Lokvarka  -  Drawdown 

No.  2 

8.9 

1.39 

1.32 

Glen  Shira  -  Drawdown  No. 

1 

0 

1.46 

1.65 

Glen  Shira  -  Drawdown  No. 

1 

12 

1.48 

1.52 

Glen  Shira  -  Drawdown  No. 

2 

0 

1.45 

1.66 

Glen  Shira  -  Drawdown  No. 

2 

12 

1.49 

1.38 

Glen  Shira  -  Drawdown  No. 

3 

0 

1.73 

1.77 

Glen  Shira  -  Drawdown  No. 

3 

12 

1.43 

1.37 

A1 cova 

54.4 

1.16 

0.95 

A1 cova 

27.2 

1.25 

1.02 

A1 cova 

18.1 

1.33 

1.09 

A1 cova 

13.6 

1.31 

1.15 

Casti letto 

56 

2.11 

2.03 

Casti 1 etto 

14 

2.15 

2.07 

Sardis  Dam 

20 

2.19 

2.19 

Sardis  Dam 

5 

2.19 

2.19 

Grenada 

18 

1.75 

1.65 

Grenada 

4.5 

2.26 

2.16 

Goscheneral p 

78.5 

2.62 

2.61 

Goscheneral p 

19.6 

2.63 

2.62 

1.28 

1.04 

0.98 

0.98 

0.97 

0.98 

0.97 

0.98 

0.98 

0.96 

0.96 

0.95 

1.13 
1.03 

1.14 
0.93 
1.02 
0.96 
0.82 
0.82 
0.82 
0.88 
0.96 
0.96 
1.00 
1.00 
0.94 
0.96 
1.00 
1.00 


pore  pressures  obtained  using  Bishop's  and  Morgenstern ' s  approach  and 
the  procedures  based  on  various  solutions  of  equations  of  hydraulic  flow 
appears  in  the  early  work  of  Glover,  Gibbs  and  Daehn  (1948).  They  first 
observed : 

" . . . .  the  removal  of  the  reservoir  load  caused  an  immediate  change 
in  both  the  contact  pressures  [effective  stress]  and  the  pore 
pressures.  The  latter  assume  a  distribution  resembling  that  due  to 
a  gravity  flow  system  but  with  a  pressure  gradient  in  the  vertical 
direction  considerably  less.  This  distribution  will  rapidly  revert 
co  the  gravity  flow  system  as  soon  as  enough  percolation  can  take 
place  to  satisfy  the  tendency  of  the  grain  structure  to  change 
volume.  The  pore  pressure  distribution  referred  to  herein  as  a 
gravity  flow  system  is  one  which  can  be  maintained  by  gravity 
forces  only  as  distinguished  from  a  system  wherein  the  pore 
pressures  are  influenced  by  an  unsatisfied  tendency  of  the  grain 
structure  to  expand." 

Although  they  failed  to  recognize  that  the  soil  may  exhibit  a  tendency 
to  compress  due  to  shear,  rather  than  expand,  Glover  et  al_,  recognized 
that  pore  pressures  will  be  influenced  by  the  tendency  of  the  soil  to 
change  volume  and  will  differ  from  those  associated  with  gravity  flow. 
However,  based  on  their  assumption  that  the  soil  would  tend  to  expand 
they  continued: 

"The  gravity  flow  system  will  exhibit  an  approximately  hydrostatic 

pressure  increase  along  a  vertical  line .  It  is  concluded  that 

no  pore  pressure  distribution  less  favorable  to  stability  than  the 
gravity  flow  system  will  be  encountered,  following  a  drawdown,  no 
matter  how  rapidly  it  may  be  accomplished." 

Thus,  while  recognizing  that  pore  pressures  will  be  influenced  by 

tendencies  for  the  soil  to  change  volume,  they  considered  the  assumption 

of  a  gravity  flow  condition,  i.  e.  that  Eq.  2.26  can  be  considered  the 

governing  equation,  should  be  conservative.  Morgenstern  (1963)  also 

noted  that  Terzaghi  and  Peck  suggested  as  a  first  approximation  that 

horizontal  flow  and  vertical  equipotential  lines  could  be  assumed  for 

conditions  after  sudden  drawdown  and  that  the  assumption  led  to 

essentially  identical  pore  pressures  to  those  obtained  using  Bishop's 
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and  Morgenstern ' s  approach.  Also  as  noted  earlier  the  transient  flow 
procedures  based  on  the  Dupuit  assumptions  lead  to  pore  pressures 
identical  to  those  obtained  by  Bishop's  and  Morgenstern 1 s  procedure. 
Consequently,  it  appears  that  pore  pressures  predicted  using  at  least  a 
number  of  the  procedures  based  on  equations  of  hydraulic  flow  are 
identical  to  those  based  on  Bishop's  and  Morgenstern ' s  approach.  It  then 
also  follows  that  the  validity  of  the  procedures  may  be  comparable  to 
the  validity  of  Bishop's  and  Morgenstern ' s  procedure. 

SUMMARY 

Effective  stress  analysis  procedures  all  require  that  the  pore 
pressures  be  estimated  for  conditions  at  the  end  of  sudden  drawdown;  the 
various  effective  stress  "procedures"  differ  only  in  the  manner  in  which 
the  pore  water  pressures  are  estimated.  Several  procedures  for 
estimating  pore  water  pressures  have  been  developed  and  discussed  in  the 
previous  sections.  Based  on  the  discussion  several  important  conclusions 
can  be  drawn: 

(1)  The  various  procedures  for  estimating  pore  water  pressures  are 
based  on  totally  different  approaches  and  assumptions  and  none  can 
be  verified  strictly  on  theoretical  grounds  due  to  the  simplifying 
nature  of  the  assumptions  employed. 

(2)  Bishop's  and  Morgenstern ' s  procedure  has  produced  good 
agreement  with  measured  field  data  including  cases  of  slope 
failure. 

(3)  Bishop's  and  Morgenstern ' s  procedure  is  the  easiest  procedure 
to  use  to  estimate  pore  water  pressures. 

(4)  Several  of  the  approximate  procedures  based  on  equations  of 
hydraulic  flow  yield  pore  pressures  which  are  identical  to  those 


obtained  by  Bishop's  and  Morgenstern ' s  procedures  and,  thus,  such 
procedures  based  on  hydraulic  flow  should  show  equally  good 
agreement  with  field  performance  data. 

(5)  None  of  the  existing  procedures  for  predicting  pore  water 
pressures  account  for  the  differences  among  various  soil  types  and 
their  tendency  to  chanae  volume  and  generate  pore  water  pressures 
during  shear.  Only  Bishop's  and  Morgenstern ' s  procedure  even 
considers  pore  pressures  generated  during  snear. 

(6)  Numerical  analyses  which  couple  the  equilibrium  equations  for 
the  stresses  with  the  governing  equation  of  hydraulic  flow  offer 
the  most  promising  alternative  for  obtaining  a  fundamentally, 
rigorous  solution  to  the  problem  of  predicting  pore  water  pressures 
due  to  sudden  drawdown.  However,  such  an  approach  is  extremely 
complex,  which  may  account  for  the  fact  that  the  approach  has  not 
been  used  to  predict  pore  pressures  due  to  sudden  drawdown. 


SECTION  3 


TWO-STAGE  ANALYSIS  PROCEDURES 
(Wright) 


INTRODUCTION 

Two  separate  sets  of  stability  calculations  are  performed  in  the 
two-stage  analysis  procedures  (Fig.  3.1).  The  first  set  of  computations 
is  performed  for  conditions  just  prior  to  drawdown.  The  purpose  of  the 
first  set  of  computations  is  to  estimate  the  stresses  under  which  the 
soil  is  consolidated  before  drawdown  and  the  corresponding  values  of 
undrained  shear  strength.  The  second  set  of  computations  is  performed 
for  conditions  immediately  after  drawdown  using  the  loads  after  drawdown 
and  the  undrained  shear  strengths  established  from  the  first  set  of 
computations.  The  two  sets  of  computations  are  performed  for  a  selected, 
individual  shear  (sliding)  surface;  various  trial  shear  surfaces  are 
selected  until  the  one  producing  the  minimum  factor  of  safety  (after 
drawdown)  is  found.  The  most  critical  shear  surface  found  in  this  manner 
will  usually  be  different  from  the  most  critical  shear  surface  for 
conditions  before  drawdown  (steady  seepage,  etc.). 

Stage  1  Stability  Computations 

The  first  set  of  stability  computations  in  the  two-stage  procedures 
is  performed  using  shear  strength  parameters  and  loads  representing 
conditions  in  the  slope  just  prior  to  drawdown.  The  shear  strength 
parameters  are  expressed  in  terms  of  effective  stresses  and  would 
normally  be  btained  from  either  consol i dated-drai ned  (CD,  S)  tests  or 
consol i dated-undrai ned  (CD,  R)  tests  with  pore  pressure  measurements. 

The  pore  water  pressures  used  in  the  first  set  of  stability  computations 
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(a)  Stage  1  -  Before  Drawdown 
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(b)  Stage  2  -  After  Drawdown 
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Figure  3.1  -  Schematic  Illustration  of  the  Two  Stages  Considered  in 
Analyses  by  Two-Stage  Procedures  for  Sudden  Drawdown 


lt>  *>. 


are  determined  either  from  a  steady  state  seepage  analysis  or  calculated 
directly  from  the  reservoir  level  before  drawdown,  neglecting  head 
losses  in  the  slope.  The  purpose  of  the  first  set  of  stability 
calculations  is  to  compute  the  stresses  on  the  shear  surface  prior  to 
drawdown.  The  factor  of  safety  calculated  in  the  first  set  of 
computations  is  of  no  interest;  however,  the  factor  of  safety  must 
ordinarily  be  calculated  as  part  of  the  limit  equilibrium  solution  for 
the  stresses.  Once  the  stresses  have  been  calculated  they  are  used  to 
estimate  the  undrained  shear  strength,  s  ,  at  the  shear  surface.  The 
undrained  shear  strength  is  estimated  for  each  slice  and,  ordinarily, 
will  be  different  for  each  slice  because  the  stresses  vary  from  slice  to 
slice.  Several  procedures  have  been  used  to  estimate  the  undrained  shear 
strength  along  the  shear  surface;  the  differences  in  these  procedures 
constitute  the  primary  difference  among  the  various  two-stage  procedures 
of  analysis  and  are  described  in  further  detail  in  the  various 
subsections  below. 

Staae  2  Stability  Calculations 


The  second  set  of  stability  computations  is  performed  for 
conditions  immediately  after  drawdown.  For  all  but  freely  draining 
materials  the  shear  strengths  used  are  those  estimated  based  on  the 
first  set  of  computations  and  are  represented  as  =  0"  strengths  (c  = 
s  ,  0  =  0).  Total  stresses  are  used;  pore  pressures  are  not  considered 
explicitly.  In  the  case  of  freely-draining  materials  effective  stresses 
are  used;  the  shear  strength  parameters  are  identical  to  those  used  in 
the  first  set  of  computations  and  pore  pressures  representing  those 
after  drawdown  are  used.  Regardless  of  whether  materials  are  freely 


draining  or  not,  any  surface  loads  due  to  water  which  remains  after 


drawdown  are  applied  for  the  second  set  of  computations.  The  factor  of 
safety  calculated  in  the  second  set  of  computations  is  the  factor  of 
safety  after  drawdown. 

CORPS  OF  ENGINEERS'  PROCEDURE 

The  Corps  of  Engineers'  procedure  is  based  on  the  use  of  a 
composite  (bilinear)  shear  strength  envelope  which  represents  the 
minimum  of  the  R  and  S  envelopes,  from  consol idated-drained  and 
consol i dated-undrai ned  shear  tests,  respectively,  as  shown  in  Fig.  3.2. 

The  Corps  of  Engineers'  Stability  Manual  (Engineer  Manual  1110-2-1902) 
indicates  that  the  same  composite  (bilinear)  shear  strength  envelope  is 
used  for  both  the  first  and  the  second  states  of  the  computations, 
although  it  is  not  clear  that  procedures  in  the  Manual  are  always  followed 
in  practice.  In  general  geotechnical  engineering  practice,  including 
possibly  the  Corps  of  Engineers'  own  practice,  it  appears  that  at  least 
three  variations  of  the  Corps  of  Engineers'  procedure  might  be  adopted  for 
the  shear  strength  used  in  the  first  stage  computations: 

1.  The  S  ("drained,"  effective  stress)  envelope  could  be  adopted 
exclusively  to  define  the  shear  strength  in  the  first  stage 
computations.  This  would  be  consistent  with  the  general  concept 
that  the  first  stage  computations  in  two-stage  analysis  procedures 
represent  conditions  with  steady-seepage  before  drawdown  and  the 
conventional  geotechnical  practice  of  employing  the  S  envelope  for 
such  conditions. 

2.  A  bilinear  envelope  representing  the  S  envelope  or  a  line  of 
"(R  +  S)/2,"  whichever  is  less,  as  shown  in  Fig.  3.3.  Such  an 
envelope  as  the  one  illustrated  in  Fig.  3.3  is  recommended  by  the 
Corps  of  Engineers  for  the  steady-seepage  stability  condition. 


figure  3.2  -  Bilinear  Shear  Strength  Envelope  Used  in  Cor| 
Sudden  Drawdown  Stability  Calculations 


Fiqure  3.3 


Thus,  use  of  the  envelope  would  again  be  consistent  with  the 
concept  that  the  first  stage  stability  computations  are  performed 
to  represent  conditions  with  steady  seepage  before  drawdown,  and 
would  be  consistent  with  the  Corps  of  Engineers'  recommendations 
for  the  steady-seepage  condition. 

3.  The  composite  (bilinear)  envelope  used  for  the  second  stage 
computations  and  shown  in  Fig.  3.2  could  be  used  for  the  first 
stage  computations  as  discussed  earlier  above  and,  apparently, 
intended  by  the  procedure  given  in  Corps  of  Engineers 
(1970). 

Inasmuch  as  the  first-stage  computations  are  only  used  to  compute  the 
effective  normal  stress  on  the  shear  surface  it  is  probably  not 
important  as  to  which  of  the  three  alternatives  is  selected  to  describe 
the  failure  envelope.  However,  some  differences  among  computed  factors 
of  safety  can  be  expected  depending  on  the  envelope  selected. 

Several  alternatives  may  be  used  in  the  first-stage  computations  to 
estimate  the  pore  pressures  in  the  slope  before  drawdown.  The  most 
fundamentally  correct  approach  is  to  estimate  pore  pressures  based  on 
either  graphical  (flow  net)  or  numerical  solutions  for  the  steady  state 
seepage  through  the  slope  prior  to  drawdown.  A  simpler  approach  is  to 
ignore  head  losses  in  the  slope  in  which  case  the  water  surface  adjacent 
to  the  slope  can  be  extended  horizontally  into  the  slope  and  treated  as 
a  horizontal  piezometric  line  to  define  pore  pressures  at  all  points 
within  the  slope.  Examples  presented  in  the  Corps  of  Engineers' 

Stability  Manual  indicate  that  head  losses  are  ignored,  i.  e.  a 
horizontal  water  surface  is  assumed  to  exist  before  drawdown.  In  this 
case  effective  stresses  can  be  accounted  for  in  either  of  two  ways:  (1) 


submerged  unit  weights  can  be  used  below  the  water  surface  and  total 


unit  weights  can  be  used  above  the  water  surface,  or  {2)  total  un’t 
weights  can  be  used  for  all  materials  and  pore  pressures  can  oe 
calculated  using  the  upstream  reservoir  surface  as  a  piezometric  surface 
for  the  entire  slope.  The  approach  employing  submerged  unit  weignts  is 
only  valid  when  there  are  no  seepage  forces  (no  hydraulic  gradient).  The 
writer  recommends  use  of  total  unit  weights  and  pore  water  pressures  in 
all  cases  because  seepage  forces  are  then  ignored  and  the  potential  for 
errors  due  to  their  inadvertent  omission  are  eliminated. 

Use  of  a  horizontal  water  surface,  rather  than  accounting  for  head 
losses  in  the  slope  due  to  steady  seepage  will  lead  to  higher  values  of 
estimated  pore  pressure  and,  thus,  lower  effective  stresses. 
Consequently,  the  undrained  shear  strengths  which  are  estimated  using 
the  effective  stresses  will  be  lower  when  a  horizontal  water  surface  is 
assumed  and,  thus,  will  be  in  error  on  the  safe  side.  The  negligiole 
head  losses  in  the  upstream  slope  of  relatively  homogeneous  slopes, 
especially  at  the  shallow  depths  where  most  drawdown  failures  may  occur, 
probably  justify  ignoring  such  head  losses.  By  ignoring  such  head  losses 
the  analyses  are  often  substantially  simplified  and  any  errors  are  in 
the  direction  of  increased  safety. 

The  first-stage  computations  are  used  to  calculate  the  effective 
normal  stresses  along  the  shear  surface.  Inasmuch  as  some  procedure  of 
slices  is  ordinarily  used,  the  effective  normal  stresses  will  De 
calculated  at  tne  base  of  each  individual  slice;  the  effective  normal 
stress  is  computed  by  dividing  the  total  norma  1  force  on  the  base  of  the 
slice  by  the  area  of  the  base  of  tne  s’ ice  and  subtracting  tne  DO>'e 


water  pressure*.  Cnee  the  effective  normal  stress  is  calculated  the 
bilinear  shear  strength  envelope  (Fig.  3.2)  is  used  to  determine  an 
appropriate  value  of  shear  strength.  The  shear  strengtn  determined  in 
this  manner  is  then  considered  to  represent  an  undrained  shear  strength 
for  use  in  the  second  stage  stability  computations. 

Second  stage  computations  are  performed  once  appropriate  undrained 
shear  strengths  have  been  estimated  from  the  first  stage  computations. 

In  general ,  the  second  stage  computations  are  performed  in  an  identical 
fashion  for  all  two-stage  analysis  procedures  once  undrained  shear 
strengths  have  been  estimated.  More  specifically,  the  Corps  of 
Engineers'  Stability  Manual  suggests  two  limit  equilibrium  analysis 
procedures  for  computing  the  factor  of  safety  for  the  second  stage:  (1) 
the  Modified  Swedish  procedure,  and  (2)  the  Ordinary  Method  of  Slices. 
Several  observations  may  be  made  with  regard  to  the  recommendation  of 
these  two  procedures: 

1.  Although  the  Ordinary  Method  of  Slices  is  not  generally 
considered  to  be  an  accurate  procedure,  in  the  special  case  of 
sudden  drawdown  (saturated  soil,  undrained  loading)  the  friction 
angle  (<*>)  is  zero  for  the  second  stage  computations.  In  this  case 
( 0  =  0)  the  factor  of  safety  calculated  using  the  Ordinary  Method 
of  Slices  will  be  identical  to  the  value  calculated  by  any 

*The  reader  is  referred  to  Wright  (1986b)  for  the  specific 
equations  used  to  compute  the  normal  force  on  the  base  of  slices  m  the 
various  limit  equilibrium  procedures  employed  in  UTEXAS2  (Wright, 

1986a).  In  the  case  where  submerged  unit  weights  are  used  for  soil  below 
the  water  table,  the  normal  force  which  is  calculated  on  the  base  of  the 
s1 i ce  will  be  the  effective  normal  force  arcl  the  P°re  water  pressures 
should  not  be  subtracted. 


procedure  wn  i  ch  sat  isMes  complete  static  equi  librium.  Thus,  the 
va’iue  f o r  tne  factor  o*  safety  should  be  as  accurate  as  any  that 
can  oe  obtained  Dy  '-m't  equilibrium  procedures. 

.  accurate  ‘c  tne  second-stage  computations,  the  factor  of 

safety  computed  by  tne  Ordinary  Method  of  Slices  will  not  be  the 
same  as  tne  one  computed  by  the  Corps  of  Engineers'  Modified 
Swedisn  or0cedure.  because  the  Modified  Swedish  procedure  does  not 
satisfy  moment  equ i '  i pr 1  urn . 

3.  Tne  Ordinary  Metnod  of  Slices  is  restricted  to  use  with  circular 
snear  surfaces.  Tnus.  the  Modified  Swedish  procedure  would 
apparently  be  used  for  noncircular  shear  surfaces  according  to  the 
Stability  Manual  (excluding  those  special  cases  where  the  Corps' 
Wedge  Method  is  used). 

4.  In  cases  where  some  portions  of  a  slope  may  be  freely  draining 
and,  thus,  shear  strengths  might  be  represented  as  "drained"  and 
expressed  in  terms  of  effective  stresses  for  the  second  stage 
computations,  the  Stability  Manual  is  unclear  regarding  use  of  the 
Ordinary  Method  of  Slices.  All  illustrations  of  applications  using 
the  Ordinary  Method  of  Slices  in  the  Stability  Manual  apply  only  to 
soils  which  are  not  drained.  The  Ordinary  Method  of  Slices  can  lead 
to  large  errors  when  effective  stresses  are  used. 

In  general,  it  is  not  clear  in  the  Corps  of  Engineers'  procedure  if  al  1 
soils  are  to  be  considered  undrained,  or  if  portions  of  the  soil  which 
are  freely  draining  are  to  be  considered  as  "drained."  Consequently,  it 
is  not  clear  regarding  the  manner  in  which  strengths  might  be 
represented  in  a  cross-section  of  a  zoned  earth  dam  for  porfions  which 
are  freely  draining.  For  example,  it  is  not  clear  if  tne  bilinear 


envelope  used  for  conventional  steady-seepage  analyses  (Fig.  3.3)  is  to 
be  used  for  freely-draining  materials  or  if  the  conventional  S  envelope 
by  itself  is  to  be  used  in  these  cases. 

LOWE  AND  KARAFIATH'S  PROCEDURE 

Lowe  and  Karafiath  (1960a)  developed  their  procedure  to  account  for 
the  fact  that  the  soil  is  consolidated  anisotropical ly  prior  to 
drawdown.  Although  Lowe  and  Karafiath  did  not  elaborate  on  the 
first-stage  stability  computations,  the  first-stage  computations  are 
generally  performed  as  conventional  steady-state  seepage  stability 
computations.  Effective  stress  shear  strength  parameters  are  determined 
employing  either  consol i dated-dra i ned  (CD,  S)  shear  test  procedures  or 
consol i dated-undrai ned  (CD,  R)  test  procedures  with  pore  pressure 
measurements;  pore  pressures  are  determined  by  either  an  appropriate 
steady-state  seepage  analysis  or  by  neglecting  head  losses  and  using  the 
reservoir  surface  as  a  horizontal  piezometric  line. 

The  principal  unique  feature  of  Lowe  and  Karafiath's  procedure 
involves  the  use  of  a  family  of  shear  strength  envelopes  to  estimate  the 
undrained  shear  strength  for  the  second-stage  stability  computations 
(Fig.  3.4).  Each  shear  strength  envelop  corresponds  to  a  given  effective 
principal  stress  ratio  at  consolidation,  Kc  (=  °ic/°3c);  the  family  of 
envelopes  is  developed  from  laboratory  tests  employing  anisotropic 
consolidation  and  several,  selected  effective  principal  stress  ratios. 
Each  envelope  expresses  the  variation  in  undrained  shear  strength  (t^) 
with  effective  normal  stress  on  the  failure  plane  at  the  time  of 
consolidation  (o^  )  for  a  given  consolidation  stress  ratio  (Kc).  The 
undrained  shear  strength  is  the  shear  stress  on  the  failure  plane  at 
failure  and  can  be  computed  from 
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Figure  3.4  -  Family  of  Shear  Strength  Envelopes  for  Various  Anis 
K  (=  v.  /T.  ),  Used  in  Lowe  and  Karafiath's  Proced 


where,  and  represent  the  major  and  minor  principal  stresses, 
respecti vely,  at  failure  and  <p  is  the  friction  angle  is  terms  of 
effective  stresses.  The  effective  normal  stress  on  the  failure  plane  at 
consolidation  can  be  computed  from: 


(olc  +  °3c}  (olc  '  °3c} 


sin  $ 


where,  o^c  and  o^c  represent  the  major  and  minor  principal  effective 

stresses  at  consolidation,  respectively.  Values  of  both  and  d^c  can 

be  calculated  from  the  data  from  any  triaxial  shear  test  and  plotted  on 

a  diagram  like  the  one  shown  in  Fig.  3.4.  The  lowest  value  of  effective 

principal  stress  ratio  for  which  failure  envelopes  can  be  determined  is 

unity,  which  corresponds  to  an  isotropic  state  of  stress  at 

consol idation  and  is  the  usual  stress  state  employed  in 

consol idated-undrained  triaxial  shear  tests.  The  upper  limit  on  the 

effective  principal  stress  ratio  to  which  a  soil  can  be  consolidated  is 

equal  to  the  effective  principal  stress  ratio  at  failure  in  a 

consol i dated-drai ned  shear  test.  Accordingly,  the  failure  envelope 

corresponding  to  the  maximum  effective  principal  stress  ratio  at 

consolidation  is  identical  to  the  conventional  (effective  stress) 

failure  envelope  from  consol  idated-drained  shear  tests  (Note:  xr  =  Trr 

f  c  f  f 

and  5f  =  in  the  case  of  consol i dated-drai ned  shear  tests). 

Fundamentally,  Lowe  and  Karafiath's  procedure  requires  that  several 


series  of  consol i dated-undrai ned  shear  tests  be  performed  using  each  of 
several  different  effective  principal  stress  ratios  for  consol i dat i on . 
However,  several  simplified  procedures  have  been  proposed  for  estimating 


the  effects  of  anisotropic  consol idation  using  the  results  of 
isotropically  consol i dated-undrai ned  triaxial  tests  with  pore  water 
pressure  measurements.  The  various  simplified  procedures  are  discussed 
in  detail  in  Section  4  of  this  paper. 

The  undrained  shear  strength,  to  be  used  in  the  second  stage 

of  the  computations  in  Lowe  and  Karafiath's  procedure  is  determined  from 
the  family  of  envelopes  like  those  shown  in  Fig.  3.4.  In  order  to 
determine  the  value  of  both  the  effective  normal  stress  on  the 
failure  plane  at  consolidation,  ,  and  the  effective  principal  stress 
ratio  at  consolidation,  °}C/°3c’  are  determined.  The  effective  normal 
stress  on  the  failure  plane  is  determined  in  the  same  manner  as  in  the 
Corps  of  Engineers'  procedure.  That  is,  the  total  normal  force  on  the 
base  of  each  slice,  computed  from  the  first-stage  stability 
computations,  is  divided  by  the  area  of  the  base  of  the  slice  and  the 
pore  water  pressure  is  then  subtracted.  The  effective  principal  stress 
ratio,  K  ,  is  not  known  and,  thus,  cannot  be  determined  as  easily  as  the 
effective  normal  stress  because  only  the  stresses  on  the  shear  plane  are 
calculated  in  limit  equilibrium  analyses.  To  calculate  the  effective 
principal  stress  ratio  Lowe  and  Karafiath  assumed  that  the  directions  of 
the  principal  stresses  are  the  same  before  drawdown  as  they  would  be  at 
failure.  Thus,  the  effective  principal  stresses  at  consolidation  act  in 
the  same  direction  with  respect  to  the  shear  surface  as  they  would  act 
if  failure  would  occur  (Fig.  3.5).  Knowing  the  effective  normal  stress 
(d^c)  and  shear  stress  ( t ^  )  on  the  shear  surface  at  consolidation,  i.e. 
before  drawdown ,  then  enables  the  principal  stresses  to  be  calculated 
(Fig.  3.6).  The  principal  stresses  are  given  by  the  following: 

°lc  =  °fc  +  Tfc(tan  T  + 


3.3 


°3c  ° 

The  shear  stress,  ,  is  the  shear  stress  on  the  shear  plane  at 

consolidation  and  is  calculated  from  the  first-stage  stability 

computations  by  dividing  the  shear  force  on  the  base  of  each  slice  by 

2 

the  respective  area  .  Once  the  principal  stresses  have  been  calculated 
the  effective  principal  stress  ratio  (K  )  is  computed  and  with  the 
effective  normal  stress  on  the  failure  plane  at  consolidation  (d^  )  also 
known  the  undrained  shear  strength  is  determined  from  the  family  of 
shear  strength  envelopes.  The  undrained  shear  strength,  t„,  becomes  the 
"cohesion"  component  of  shear  strength  for  the  second-stage  stability 
computations  and  <p  is  set  equal  to  zero  (c  =  q  -  0).  Once  the  shear 

strength  is  determined,  the  second-stage  stability  computations  are 
fundamentally  the  same  as  those  by  any  other  two-stage  analysis 
procedure,  including  the  Corps  of  Engineers  procedure. 

Some  difficulties  may  be  encountered  with  Lowe  and  Karafiath's 
procedure  for  soils  where  an  effective  cohesion  intercept  (c)  exists.  In 
such  cases,  the  maximum  effective  principal  stress  ratio  associated  with 
the  consol i dated-drai ned  (effective  stress)  envelope  is  not  a  constant, 
but  varies  with  the  effective  stress.  As  the  effective  stresses  decrease 
and  become  small,  the  effective  principal  stress  ratio  at  failure  first 
approaches  a  value  of  infinity  as  d^  approaches  zero  (Fig.  3.7a)  and, 
then,  the  effective  principal  stress  ratio  becomes  negative;  the 


The  reader  is  again  referred  to  Wright  (1986b)  for  the  specific 
equations  used  to  compute  the  shear  force  on  the  base  of  slices  in  the 
various  limit  equilibrium  procedures  employed  in  UTEXAS2. 


effective  principal  stress  ratio  is  clearly  negative  when  reaches 
zero  (Fig.  3.7b).  Any  attempt  to  interpolate  shear  strength  values 
between  strength  envelopes  corresponding  to  the  maximum  effective 
principal  stress  ratio  (K^),  and  some  lower  effective  principal  stress 
ratio  becomes  complicated  when  the  value  of  becomes  very  large  or 
negative.  In  addition  to  the  problems  associated  with  defining  for 
the  failure  envelope,  either  one  or  both  of  the  effective  principal 
stresses  (at  consol idation)  which  are  calculated  from  Eqs.  3.3  and  3.4 
using  the  stresses  computed  from  the  first-stage  stability  computations 
may  also  become  negative.  In  order  to  avoid  these  difficulties  several 
alternatives  exist  and  one  or  more  may  be  necessary  to  overcome  these 
problems.  First,  low  stresses  and  potentially  negative  values  of  the 
correspond! ng  principal  stresses  calculated  from  the  stability  analyses 
may  be  at  least  partially  eliminated  by  introducing  a  tension  crack  into 
the  stability  computations.  Secondly,  with  regard  to  computing  the 
effective  principal  stress  ratio  associated  with  the  effective  stress 
failure  envelope,  Wong,  Duncan  and  Seed  (1983),  neglected  the  effective 
cohesion  intercept  and  assumed  a  constant  effective  principal  stress 
ratio  at  failure  given  by, 

°lf  2 

K,  =  —  =  tan  ( 45  +  p/2)  3.5 

f  °3f 

In  addition  the  effective  cohesion  intercept  (c)  may  be  neglected 
entirely  in  the  stability  computations,  including  the  first-stage 
computations.  This  would  both  eliminate  the  possibility  of  '-egatne 
values  being  computed  for  the  stresses  in  tne  first-stage  stability 
computations  and  the  effective  principal  stress  ratio  associated  w’tn 
the  effective  stress  failure  envelooe  wou’d  be  a  constan-  m t’”;re 


quantity.  None  of  the  alternatives  in  Lowe  and  Karafiath's  procedure  is 
ideal  nor  have  any  been  investigated  in  sufficient  detail  at  the  present 
time  to  warrant  recommending  one  in  preference  to  another.  However,  the 
effective  cohesion  intercept,  c,  often  has  a  very  important  effect  on 
the  computed  factor  of  safety  for  small  slopes  and  shallow  slides,  and 
the  manner  in  which  the  effective  cohesion  is  treated  may  have  an 
important  effect  on  the  computed  factors  of  safety  for  sudden  drawdown. 

DISCUSSION 

Both  two-stage  procedures  are  conceptually  similar.  Both 
involve  calculating  the  effective  stresses  before  drawdown  from  a  set  of 
stability  calculations  for  the  steady  seepage  condition  and,  then,  using 
the  stresses  to  estimate  undrained  shear  strengths  at  the  time  of 
drawdown.  A  second  set  of  stability  calculations  is  then  performed  using 
the  estimated  undrained  shear  strengths.  The  only  major  difference 
between  the  two  approaches  is  in  the  procedures  used  to  estimate  the 
undrained  shear  strengths.  In  a  broad  sense,  the  Corps  of  Engineers' 
procedure  represents  a  conservative,  lower-bound  approach  to  the  one 
suggested  by  Lowe  and  Karafiath.  However,  the  use  of  the  R-type  failure 
envelope  in  the  Corps  of  Engineers'  procedure  warrants  some  further 
di scussion . 

R  envelope  versus  line.  The  R  envelope  is  obtained  Py  first 

constructing  a  series  of  "Mohr's"  circles,  where  the  minor  principal 
stress  represents  tne  minor  principal  ef  f ecti ve  stress  at  the  time  of 
consolidation  ( .1  ) ,  and  the  diameter  of  the  circle  represents  the 

principa1  stress  difference  at  the  time  of  failure  as  shown 

in  Fig.  3.8.  The  "Mohr's"  circles  and  cor  respond i nq  failure  envelopes 
constructed  using  such  circles  are  referred  to  in  most  textbooks  and 


much  of  the  literature  as  "total  stress"  circles  and  envelopes, 
respectively.  However,  this  nomenclature  is  not  correct  inasmuch  as  the 
minor  principal  stress  is  the  effective  stress  at  the  time  of 
consol idation,  while  the  diameter  of  the  circle  represents  values  at 
f ai 1 ure .  The  stresses  being  used  are  not  actually  total  stresses  nor  do 
they  occur  simultaneously;  the  stresses  are  not  associated  with  a  given 
state  of  stress.  Circles  like  the  ones  plotted  for  the  R  envelope  are 
only  valid  Mohr's  circles  when  the  stresses  occur  at  the  same  instant. 
The  circles  would  be  a  valid  representati on  of  the  total  stresses  at 
failure  provided  that  no  back  pressure  is  used  in  the  laboratory  test  or 
no  pore  pressure  exists  in  the  field.  This  is  seldom  the  case  and,  thus, 
the  circles  in  general  are  invalid.  For  the  remainder  of  this  paper 
circles  like  the  one  shown  in  Fig.  3.8  will  be  referred  to  as  "circles 
of  stress,"  rather  than  Mohr's  circles. 

Although  the  circles  of  stress  on  which  the  R  envelope  is  based 
have  no  fundamental  basis,  the  envelope  which  is  drawn  tangent  to  the 
circles  may  possess  close  similarity  to  the  line  representing  the 
relationship  between  and  o^c  for  the  case  of  isotropic 
consolidation.  The  similarity  can  be  illustrated  by  referring  to  Fig. 
3.9.  The  circle  of  stress  shown  is  the  one  normally  constructed  for 
obtaining  the  R  envelope.  Line  A- A '  is  drawn  tangent  to  the  circle  to 
represent  an  R  envelope  (Only  one  circle  is  shown  for  clarity,  although 
several  circles  are  required  to  draw  the  failure  envelope).  The  stresses 
on  the  failure  plane  at  consolidation,  °fc  (  =  in  the  case  of 
isotropic  consolidation)  and  the  shear  stress  on  the  failure  plane  at 
failure.  are  also  indicated  in  this  figure.  Finally  the  line  B - B ' 

is  shown  drawn  through  points  with  rhe  coordinates  ^ ■  if  can  be 


seen  that  the  lines  A-A'  and  B - B '  are  similar.  In  order  to  examine  the 


degree  to  which  the  two  lines  are  similar  to  one  another  a  senes  oh 
calculations  were  performed  in  which  the  slope  and  intercept  values  for 
the  R  envelope  were  compared  with  the  corresponding  values  for  the  line 
representing  the  relationship  between  •  ^  and  The  intercept  (iQ) 

and  slope  angle  (b  )  fon  the  line  representing  the  relationsnip  between 
and  dfc  a-'e  related  to  the  intercept  and  slope  of  the  R  envelope,  cR 
and  respectively,  and  the  value  of  the  effective  stress  friction 
angle.  <j:  The  intercept,  ,  is  given  by. 

(intercept)  =  c^N^  cos  <t  3.6 

and  the  slope  angle.  is  given  by.. 

b  (slope  angle)  =  tan  *[4(N  -  i)  cos  $1  3.7 

T  c  <t 

where . 


N  -  tan4"  ( 45  +  >  3.8 

■t 

Values  of  the  intercept  and  slope  of  the  line  are  shown  for  a 

relatively  wide  range  ’n  values  of  cR,  jR  and  in  Table  3.1.  The 
percentage  difference  between  the  intercept  of  the  .f  line.  and 

the  "cohesion",  cR,  for  tne  R-en-„  e  '■  ope .  and  the  cor  respond 1  ng 
differences  in  the  s’ooe  ang'-es.  ••  and  *"or  the  two  envelopes  are 
also  shown.  The  differences  showr  are  ’''dependent  of  the  cohesion  value 
(cR)  and  become  1  arge-  -  •  'eas'ng'v  pns'f  ve  as  tne  difference  between 
the  friction  angles,  ;.  3-0  ;.  nr 'eases.  D  ’  *  ‘ererces  of  a'r'jst  80 
percent  '  n  magn  '  t  .<  te  ~  a  ,  jf  *re  R  er.e  "ope  w  ’  ’  '  •  e  ■  d  ’Owe’ 

strengths  than  the  .  •••»  •  o*-  ‘re  -  or  >  by  c*  '  asec  ..sf  r  at  eo  . 


However .  there  are  at 
where  tne  R  er.e  ;■»* 
e — 1 1  r  e  I  v  evident  mv 


>  .,T!e  a  s  e  s  ;  t  -  1  degrees  .  :  A-.  ot-  ;• 


re  s  •'  e  °  q  *  n  s  A  .•  *  •  s  rc' 

*  iny-np.  rc  p.  p  p  p  n  '  ..  ses  t  r  e  A 


TABLE  3.1 


COMPARISON  OF  THE  SLOPE  AND  INTERCEPT  FOR  THE 


vnr  nr\i  Jvn  ui  iiii_  JLur  i_  nnu  1111  i  iwi\  i  i  iu  1  ■f  C  lh  v 

WITH  THE  STRENGTH  PARAMETERS  FROM  THE  R  ENVELOPE  (£„  and 


ENVELOPES 


CR 

*R 

4> 

To 

e 

T 

CR  ‘  To 

T 

o 

0R  ' 

6 

(psf ) 

(degs) 

(degs) 

(psf) 

(degs) 

(%) 

(% 

50 

10 

25 

54.0 

10.8 

-7.4 

-7 

300 

10 

25 

324.0 

10.8 

-7.4 

-7 

700 

10 

25 

756.1 

10.8 

-7.4 

-7 

1000 

10 

25 

1080.1 

10.8 

-7.4 

-7 

3000 

10 

25 

3240.3 

10.8 

-7.4 

-7 

50 

20 

25 

64.7 

25.2 

-22.7 

-20 

300 

20 

25 

388.3 

25.2 

-22.7 

-20 

700 

20 

25 

906.0 

25.2 

-22.7 

-20 

1000 

20 

25 

1294.3 

25.2 

-22.7 

-20 

3000 

20 

25 

3883.0 

25.2 

-22.7 

-20 

50 

10 

35 

48.8 

9.8 

2.4 

2 

300 

10 

35 

292.9 

9.8 

2.4 

2 

700 

10 

35 

683.4 

9.8 

2.4 

2 

1000 

10 

35 

976.2 

9.8 

2.4 

2 

3000 

10 

35 

2928.7 

9.8 

2.4 

2 

50 

20 

35 

58.5 

23.1 

-14.5 

-13 

300 

20 

35 

351.0 

23.1 

-14.5 

-13 

700 

20 

35 

818.9 

23.1 

-14.5 

-13 

1000 

20 

35 

1169.9 

23.1 

-14.5 

-13 

3000 

20 

35 

3509.6 

23.1 

-14.5 

-13 

50 

30 

35 

70.9 

39.3 

-29.5 

-23 

300 

30 

35 

425.6 

39.3 

-29.5 

-23 

700 

30 

35 

993.2 

39.3 

-29.5 

-23 

1000 

30 

35 

1418.8 

39.3 

-29.5 

-23 

3000 

30 

35 

4256.4 

39.3 

-29.5 

-23 

50 

10 

45 

42.1 

8.5 

18.7 

18 

300 

10 

45 

252.8 

8.5 

18.7 

18 

700 

10 

45 

589.9 

8.5 

18.7 

18 

1000 

10 

45 

842.7 

8.5 

18.7 

18 

3000 

10 

45 

2528.1 

8.5 

18.7 

18 

50 

20 

45 

50.5 

20.2 

-1.0 

-0 

300 

20 

45 

303.0 

20.2 

-1.0 

-0 

700 

20 

45 

706.9 

20.2 

-1.0 

-0 

1000 

20 

45 

1009.9 

20.2 

-1.0 

-0 

3000 

20 

45 

3029.6 

20.2 

-1.0 

-0 

envelope  will  always  be  more  conservative  than  Lowe  and  Karafiath's 


procedure,  which  uses  the  more  correct  line.  The  differences 

between  the  two  procedures  will  depend  on  the  values  of  and  c  as  well 
as  on  the  effective  principal  stress  ratio  (K  )  in  the  case  of  Lowe  and 
Karafiath's  procedure. 

It  is  fortuitous  that  the  R  envelope  and  the  line  representing  the 
relationship  between  and  d^c  are  as  similar  as  they  are  over  a  broad 
range  of  values.  This  close  similarity  was  recognized  by  Johnson  (1974) 
and  provided  a  portion  of  the  the  rationale  for  adoption  of  the  R 
envelope  in  the  Corps  of  Engineers'  procedure.  The  numerical  comparisons 
presented  in  Table  3.1  suggest  that  the  R  envelope  should  generally  lead 
to  conservative  estimates  of  shear  strength  relative  to  the 
fundamentally  more  correct  vs.  d^c  relationship,  although  exceptions 
to  this  are  also  shown. 

R  envelope  construction .  At  least  two  different  criteria  appear  to 
be  currently  used  in  constructing  the  R  enveloDe  from  a  series  of 
consol i dated-undra i ned  triaxial  shear  tests.  The  first  criterion  is  the 
one  util : zed  above  and  consists  of  drawing  the  failure  envelope  tangent 
to  the  circles  of  stress  as  illustrated  previously  in  Fig.  3.8.  The 
second  criterion  consists  of  constructing  the  envelope  by  drawing  a  line 
through  points  on  the  circle  of  stress  cor  respond i ng  to  stresses  on  tre 
failure  plane  as  shown  in  Fig.  3.10.  The  Corps  of  Engineer's  Stability 
Manual  suggests  that  both  approaches  are  used:  tnve'opes  a-e  fawn 
tangent  to  the  circles  of  stress  •  n  the  case  ot  "  ..ng  ’  ^tft^a"  j 

through  points  CO’ respond  ’  ng  tc  stresses  :,r  *  *  a  •  '  «  ;• 


a 

1 

b 

a 

3c 

◄ - 

_ ^ 

3c 

◄ - 

- fr. 

(a  -  a  )"  (a  -  a  ).b 

'  1  3,f  '  1  3,f 


a 


3.1U  -  R- Envelope  Constructed  through  Points  Corresponding  to 
Stresses  on  the  failure  Plane 


stresses  on  the  failure  plane.  Many  textbooks  show  the  R  envelope 
(consol idated-undrained  "total  stress  envelope")  as  apparently  being 
drawn  tangent  to  the  circles  of  stress  (e.  g.  Sowers,  1979;  Holtz  and 
Kovacs,  1981;  Terzaghi  and  Peck,  1967).  There  does  not  appear  to  be  any 
universally  accepted  procedure  for  constructing  R  envelopes  and  there 
probably  never  will  be  due  to  the  rather  arbitrary  nature  of  the 
envel ope . 

The  relationship  between  the  intercepts  and  slopes  of  R  envelopes 
drawn  tangent  to  the  circles  of  stress  and  drawn  through  points 
corresponding  to  the  stresses  on  the  failures  plane  can  be  derived 
theoretically.  The  relationship  between  the  "cohesion"  intercepts  is 
given  by, 

COS0  COS0D 

r  =  r  _ _  o  q 

Rf  R  1  -  sin<ji  sin0R 


where,  cR^  represents  the  intercept  of  the  envelope  drawn  through  points 


corresponding  to  stresses  on  the  failure  plane,  cR  and  pR  represent  the 


intercept  and  slope,  respectively,  of  the  envelope  drawn  tangent  to  the 
circles  of  stress,  and  $  is  the  effective  stress  friction  angle.  The 
relationship  between  the  "friction  angles"  is  given  by, 

cos$  sin<|i. 


KR 


0n*  =  tan  [■,  .  - 

Rf  L1  -  sin0  sin0 


-] 


3.10 


R 


where,  represents  the  slope  of  the  envelope  drawn  through  points 


cor-espording  to  stresses  on  the  failure  plane.  Values  of  cR.p  and 


are  shown  in  Table  3.2  for  a  range  in  values  of  cR,  <j>R  and  <j.  The 


percentage  differences  between  the  two  sets  of  intercepts  ( cD  and  cc 


he  two  sets  of  slopes  ($R  and  0D4;)  are  also  shown  in  this  table. 


Rf ' 


iterances,  expressed  as  percentages  are  independent  of  the 


TABLE  3.2 


COMPARISON  OF  THE  SHEAR  STRENGTH  PARAMETERS  FOR  R-ENVELOPES  DRAWN 
TANGENT  TO  CIRCLES  OF  STRESS  (cB,  <t>R)  AND  R-ENVELOPES  DRAWN  THROUGH 
POINTS  CORRESPONDING  TO  STRESSES  Oft  THE  FAILURE  PLANE  (cRf,  *Rf) 


CRf  "  CR 

$Rf  -  ■ 

CR 

*R 

CRf 

*Rf 

CR 

*R 

(psf ) 

(degs) 

(degs) 

(psf) 

(degs) 

(%) 

(%) 

50 

10 

25 

48.2 

9.6 

-3.8 

-3.7 

300 

10 

25 

289.0 

9.6 

-3.8 

-3.7 

700 

10 

25 

674.3 

9.6 

-3.8 

-3.7 

1000 

10 

25 

963.2 

9.6 

-3.8 

-3.7 

3000 

10 

25 

2889.7 

9.6 

-3.8 

-3.7 

50 

20 

25 

49.8 

19.9 

-0.4 

-0.4 

300 

20 

25 

298.7 

19.9 

-0.4 

-0.4 

700 

20 

25 

696.9 

19.9 

-0.4 

-0.4 

1000 

20 

25 

995.6 

19.9 

-0.4 

-0.4 

3000 

20 

25 

2986.7 

19.9 

-0.4 

-0.4 

50 

10 

35 

44.8 

9.0 

-11.6 

-11.4 

300 

10 

35 

268.8 

9.0 

-11.6 

-11.4 

700 

10 

35 

627.2 

9.0 

-11.6 

-11.4 

1000 

10 

35 

895.9 

9.0 

-11.6 

-11.4 

3000 

10 

35 

2687.8 

9.0 

-11.6 

-11.4 

50 

20 

35 

47.9 

19.2 

-4.4 

-4.1 

300 

20 

35 

287.3 

19.2 

-4.4 

-4.1 

700 

20 

35 

670.3 

19.2 

-4.4 

-4.1 

1000 

20 

35 

957.6 

19.2 

-4.4 

-4.1 

3000 

20 

35 

2872.8 

19.2 

-4.4 

-4.1 

50 

30 

“*5 

49.7 

29.9 

-0.5 

-0.4 

300 

30 

35 

298.4 

29.9 

-0.5 

-0.4 

700 

30 

35 

696.3 

29.9 

-0.5 

-0.4 

1000 

30 

35 

994.7 

29.9 

-0.5 

-0.4 

3000 

30 

35 

2984.0 

29.9 

-0.5 

-0.4 

50 

10 

45 

39.7 

8.0 

-26.0 

-25.5 

300 

10 

45 

238.2 

8.0 

-26.0 

-25.5 

700 

10 

45 

555.7 

8.0 

-26.0 

-25.5 

1000 

10 

45 

793.8 

8.0 

-26.0 

-25.5 

3000 

10 

45 

2381.5 

8.0 

-26.0 

-25.5 

50 

20 

45 

43.8 

17.7 

-14.1 

-13.0 

300 

20 

45 

262.9 

17.7 

-14.1 

-13.0 

700 

20 

45 

613.5 

17.7 

-14.1 

-13.0 

1000 

20 

45 

876.4 

17.7 

-14.1 

-13.0 

3000 

20 

45 

2629.3 

17.7 

-14.1 

-13.0 

50 

30 

45 

47.4 

28.7 

-5.6 

-4.6 

300 

30 

45 

284.2 

28.7 

-5.6 

-4.6 

700 

30 

45 

663.1 

28.7 

-5.6 

-4.6 

1000 

30 

45 

947.3 

28.7 

-5.6 

-4.6 

3000 

30 

45 

2841.9 

28.7 

-5.6 

-4.6 

cohesion  values  and  increase  as  the  difference  between  the  friction 
angle  from  the  R  envelope,  $R,  and  the  effective  stress  envelope,  <j>, 
increase.  Differences  of  over  25  percent  are  shown  in  the  extreme  case 
(4>r  =  10  degrees,  $  =  45  degrees). 

Use  of  an  R  envelope  which  is  drawn  through  points  on  the  circles 
of  stress  corresponding  to  stresses  on  the  failure  plane  appears  to  be 
unnecessari ly  conservative  in  instances  when  the  envelope  is  to  then  be 
used  to  represent  the  relationship  between  and  .  The  combined 
effect  of  constructing  the  R  envelope  through  points  representing 
stresses  on  the  failure  plane  and  using  the  envelope  to  express  the 
relationship  between  and  will  be  to  underestimate  the  undrained 
shear  strength  in  almost  all  cases.  This  underestimate  coupled  with  the 
additional  conservatism  associated  with  the  use  of  isotropic 
consolidation  and  an  envelope  representing  a  minimum  of  the  R-S 
strengths  appears  to  excessively  compound  conservatism. 


SECTION  4 


ESTIMATION  OF  ACU  STRENGTH  ENVELOPES  FROM  ICU  TESTS 

(Wright) 

Lowe  and  Karafiath's  (1960a)  procedure  employs  a  series  of 
envelopes,  each  envelope  corresponding  to  a  particular  effective 
principal  stress  ratio  at  consolidation,  <c  (=  5^  /o^)-  enve^°Pes 

may  be  obtained  by  performing  several  series  of  consol i dated-undrai ned 
triaxial  shear  tests  employing  anistotropic  stress  for  consolidation  and 
with  each  series  of  tests  being  performed  at  a  selected  effective 

i 

i> 

principal  stress  ratio.  Testing  of  this  type  is  laborious,  expensive  and 
often  the  scatter  among  results  at  different  effective  principal  stress 
ratios  leads  to  considerable  uncertainty  in  the  strength  envelopes  even 
after  extensive  testing  has  been  completed.  As  an  alternative  several 
procedures  have  been  proposed  for  estimating  the  ACU  envelopes  from  the 
results  of  ICU  triaxial  shear  tests.  Wong,  Duncan  and  Seed  (1983) 
examined  the  available  procedures  for  estimating  the  strength  envelopes 
from  the  results  of  tests  employing  isotropic  consolidation.  They  also 
compared  the  differences  among  the  failure  envelopes  using  data  obtained 
from  a  series  of  tests  performed  on  the  soil  from  Hirfanli  Dam  and 
presented  by  Lowe  and  Karafiath  (1960b). 

The  various  procedures  for  estimating  the  undrained  shear  strength 
for  anisotropic  consolidation  from  the  results  of  tests  employing 
isotropic  consolidation  are  briefly  restated  in  this  chapter.  The 
appropriate  equations  used  to  estimate  the  failure  envelopes  for  various 
values  of  the  effective  principal  stress  ratio  for  anisotropic 
consolidation  are  included.  Most  of  these  equations  were  originally 
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presented  by  Wong,  Duncan  and  Seed  (1983).  In  several  cases  they 
considered  the  effective  cohesion,  c,  to  be  zero;  in  other  cases 
non-zero  cohesion  was  accounted  for.  All  of  the  equations  presented 
herein  for  the  various  procedures  have  been  derived  for  the  case  where 
an  effective  cohesion  value  may  exist.  The  equations  presented  by  Wong, 
Duncan  and  Seed  also  were  all  based  on  the  assumption  that  the  R 
envelope  passes  through  points  on  the  circles  of  stress  corresponding  to 
stresses  on  the  failure  plane  rather  than  being  tangent  to  the  circles 
of  stress.  The  convention  used  by  Wong,  Duncan  and  Seed  has  been 
retained  for  this  section  of  the  report  to  permit  the  work  in  this 
report  to  be  compared  to  their  earlier  work. 

Useful  insight  into  the  various  procedures  for  estimating  failure 
envelopes  for  anisotropic  consolidation  was  gained  during  this  study  by 
examining  the  values  of  Skempton's  pore  pressure  coefficient,  R^,  which 
are  implicit  in  each  of  the  interpolation  procedures.  In  general  the 
1 owtr  the  value  of  the  pore  pressure  coefficient,  the  higher  will  be  the 
value  of  the  corresponding  undrained  shear  strength.  The  appropriate 
equations  for  computing  the  pore  pressure  coefficients  were  derived  in 
conjunction  with  the  preparation  of  the  present  paper  and  are  also 
presented  in  this  section. 

Comparisons  were  made  among  both  the  strength  envelopes  (1^-0^  ) 
and  pore  pressure  coefficients  for  the  various  procedures  for  several 
different  soils.  These  comparisons  are  presented  at  the  end  of  this 
section . 

TAYLOR'S  PROCEDURE 

Taylor  assumed  that  the  effective  stress  paths  for  anisotropic 
consolidation  are  coincident  with  the  effective  stress  paths  for 


isotropic  consolidation  for  the  portion  of  the  loading  which  is 
undrained,  i.  e.  after  consolidation.  Accordingly,  the  undrained  shear 
strength  obtained  with  isotropic  consolidation  will  be  identical  to  the 
undrained  shear  strength  obtained  for  any  specimen  which  is  consolidated 
to  an  effective  stress  state  existing  at  any  stage  in  the  ICU  test. 
versus  envelopes . 

Points  corresponding  to  the  line  for  a  given  anisotropic 

consolidation  stress  ratio,  K  ,  are  obtained  from  the  effective  stress 
paths  for  specimens  which  were  isotropically  consolidated  by  identifying 
the  points  on  the  stress  paths  where  the  effective  principal  stress 
ratio  (d^/dj)  is  equal  to  the  desired  stress  ratio  (K  ).  The  effective 
stresses  corresponding  to  the  selected  point  on  the  effective  stress 
path  ( d^ ,  d^c )  are  then  considered  to  represent  the  effective  stress  at 
consolidation  for  the  selected  anisotropic  consolidation  stress  ratio; 
the  shear  strength  corresponding  to  the  selected  stress  path  from  the 
isotropically  consolidated  specimen  is  then  considered  to  also  be  the 
shear  strength  for  the  anisotropical ly  consolidated  specimen. 

Taylor's  procedure  formally  requires  that  the  actual  stress  path  be 
used  to  estimate  the  strength  envelopes.  However,  if  the  effective 
stress  path  is  assumed  to  be  linear  {A  is  constant  during  shear)  the 
envelopes  for  anisotropic  consol idation  can  be  estimated  directly  from 
the  R  envelope  and  effective  stress  shear  strength  parameters,  c  and  <j>. 
To  estimate  the  failure  envelopes  for  a  linear  stress  path,  the 
effective  minor  principal  stress  at  consol i dati on  corresponding  to  the 
equivalent  test  with  isotropic  consolidation  must  first  be  found.  The 
equivalent  isotropic  consolidation  pressure  will  be  termed,  j ,  and 
can  be  computed  by  solving  the  following  equation: 


4.1 


d,  +  “ 
1  0 


O  ^  t  O 

3c  I  3c 


3c  I 


+  d,  to3cI(Kc  " 


(Kr  -  l)d? 
i  MKc  -  l)d1  +- 


j0  .  +  d_ 

3c  1  3 


where,  a,  is  the  effective  consolidation  pressure  for  an 
an i sotropTcal 1 y  consolidated  specimen  and, 

dj  =  j  cos$(cotp^^  -  coteji) 

d^  =  ^(ccot*  -  c^cot<j>^ )  ( cos$cot<ji^  +  sin<j>  -  1) 
d3  =  cRf cotdRf 


4.2 

4.3 

4.4 


The  undrained  shear  strength,  is  then  computed  from  the  following 

equation : 


ff  1 


q3cl  tandRf  +  cRf 
(1  -  sin$)tan<j>n.r/cos<ti 


4.5 


where  o^  .  is  determined  by  solving  Eq.  4.1. 

Pore  Pressure  Coefficients,  . 

Once  the  equivalent  effective  minor  principal  stress,  o^j,  is 
calculated  for  any  values  of  o3c  (or  o^c)  and  Kc  as  described  above, 
the  corresponding  pore  pressure  coefficient  can  be  calculated  from  the 
following  equation,  which  is  obtained  for  isotropic  consolidation: 


Af  =  d.  +  = 
'  1  o 


3d  +  d3 


4.6 


The  equivalent  consolidation  pressure,  °3cj>  is,  again,  obtained  from 
the  solution  to  Eq.  4.1. 

LOWE  AND  KARAFIATH'S  PROCEDURE 


Lowe  and  Karafiath  expressed  the  pore  water  pressures  at  failure  by 
an  equation  of  the  form, 

AUf  =  ( Ao if  "  Ao3f)  +  ^2 
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where  ip^"  and  "ip^"  are  pore  pressure  coefficients  .  The  pore  pressure 
coefficients  are  determined  from  ICU  tests  and  are  related  to  the  R 
strength  parameters  (cR^  and  <pR^)  as  follows, 

1(1.  =  ^cos0(cot<(iR^  -  cot<p)  4.7 

>l>2  -  ~  (cRf  cot*Rf  -  c  cot*)  4.8 

Lowe  and  Karafiath  assumed  that  the  pore  pressure  coefficients  and 

'>2"  are  the  same  f°r  isotropic  and  anisotropic  consolidation. 
t versus  envelopes . 

The  relationship  between  the  undrained  shear  strength  and  the 
effective  stress  on  the  failure  plane  at  consolidation  is  given  by, 

2of  c 

C(K  tlHI  -  l)sin»]  *  5[(1  +  el>  cot?]  •  e2 


<1 +  n) cot* +  tani>  ■  hit 


where , 


N-  -  K 


1  +  sin*  '^1 


1  ’  1  +  (Kc  -  l)i>1 


N_  -  K 


C£C0S<>  1  +  sTn*  +  2  N*>1  +  *2 


1  +  (Kc  -  1)^1 


N  =  ta'i2(45  + 


^Lowe  and  Karafiath  originally  used  the  designation,  A^,  for 
However,  the  designation  ip.  has  been  adopted  for  this  paper  to  avoid 
confusion  with  Skempton's  pore  pressure  coefficient,  A,  A,  A^,  etc.  The 
parameter,  *  is  not  equivalent  to  any  of  Skempton's  pore  pressure 
coeff  i  dents1 


Pore  Pressure  Coefficients,  A 


Pore  pressures  coefficients  are  calculated  by  first  calculating  tne 
pore  pressure  at  failure  from  the  following  equation: 


iu. 


(°3c^  +  '  Kca3c)u’l  1 


1  -  (i  -  N  )<v 
0  1 
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The  corresponding  difference  in  the  principal  stress  changes,  c  - . f  - 
comPuted  from  the  following  equation: 

“if  -  “3f  =  <°3c  •  “f><\  ‘  >>  *  -  »3c|Kc  •  »  4-!4 

The  pore  pressure  coefficient  is  then  obtained  by  dividing  the  cnange  in 

pore  pressure  at  failure  by  the  difference  in  the  principal  stress 

changes  as, 

AUr 


f  talf  -  Ao3f 


4.15 


MODIFIED  KARAFIATH  PROCEDURE 

This  procedure  was  suggested  by  Wong,  Duncan  and  Seed.  They  assumed 
that  the  relationship  between  Skempton’s  pore  pressure  coefficient  A* 


and  the  effective  stress  on  the  failure  plane  at  consolidation,  of„,  was 


independent  of  the  effective  principal  stress  ratio  at  consolidation 


(K  ).  Accordingly,  the  relationship  between  the  pore  pressure 


coefficient  and  was  uniquely  defined  and  could  be  calculated  from 
the  results  of  tests  performed  on  isotropically  consolidated  specimens. 
versus  a^c  envelopes. 

The  following  equation  has  been  derived,  expressing  the 
relationship  between  the  undrained  shear  strength,  t  ,  ^ ,  and  the 


effective  stress  on  the  failure  plane  at  consolidation,  a 


f  c ' 


2o 


f  c 


C(Kc  +  1)  -  (Kc 
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ff 
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(1  +  a3)  cot^  +  tan<j>  -  CQ-^- 


4.18 


1  "  (1  +  si n<j>) [  1  +  A(Kc  -  1)] 

Acv  N  [(1  -  sin<j)K  +1+  si  no] 
,,  -  - $ - - - 

2  ( 1  +  sinip) [1  +  A(  K  -  1 ) ] 


and  N-  is  given  by  Eq.  4.12. 
0 


Pore  Pressure  Coefficients,  Af 


For  isotropically  consolidated  specimens  the  pore  pressure 
coefficient  at  failure  is  expressed  in  terms  of  the  effective 
consolidation  pressure,  ,  the  strength  parameters  from  the  R 
envelope,  c^  and  <j>-p  and  the  effective  stress  friction  angle, 
d„ 


Af  =  d.  +  - 
t  1  o 


3c  +  d3 


by 

4.19 


where,  d^,  d^  and  d^  are  give  by  Eqs.  4.2,  4.3  and  4.4,  respectively. 
Inasmuch  as  Eq.  4.19  is  written  for  the  case  of  isotropic  consol i oat i on . 


the  consolidation  pressure,  ,  is  also  equal  to  the  effective  normal 


stress  on  the  failure  plane  at  consolidation.  In  the  case  of  anisotropic 
consolidation  Eq.  4.19  is  rewritten  by  replacing  the  effective  stress, 

,  by  the  effective  normal  stress  on  the  failure  plane  at 
consolidation,  b^c .  Thus,  the  expression  for  A^  now  becomes. 


Af  =  d,  +  — 
t  1  o 


fc  +  d3 


4.20 


The  effective  normal  stress  on  the  failure  plane  at  consolidation  in  Eq. 
4.20  can  be  expressed  in  terms  of  the  effective  minor  principal  stress 


at  consolidation  (o^  )  and  the  effective  principal  stress  ratio  at 


consol i dati on  by , 


Thus,  Eq.  4.20  can  r-e  '  ■- adulate  A,  for  any  values  of  c,  and  K  . 

f  3c  c 

E^cVi  -  NOORANY  AND  SEED'S  STJDY 

This  approacn  was  digested  by  Wong,  Duncan  and  Seed,  based  on 

earlier  work  of  Noc'd'-,  }• -,  ~ea.  hhe  procedure  is  based  on  the 

assumption  that  the  >*  ■  r  ; o  petween  the  pore  water  pressure 

coefficient  ana  tre  •  •  •  •  \  •  .  e  normal  stress  on  the  failure  plane  at 

consolidation  is  the  sane  *  •  ’sotropic  and  anisotropic  consolidation. 

Accordingly,  the  assumpf  employed  in  this  approach  are  identical  to 

those  used  in  the  Modi*'-.:  -  wat'ath  procedure:  only  the  form  of  the 

algebraic  equations  usee  *  impute  various  quantities  are  different. 

’ versus  enveippes  . 

The  undrained  shear  s*  •  ■•'gin  ;s  computed  from: 

2o^cCOS(isi  n  ;  r  i  -  A  ^  (  K  -  1)] 

l(K  +  :  )  -  -  1  )smc] 

,  =  - S - - - ■  . -  4.22 

tf  1  +  ■  ils'ni 

where,  the  pore  pressure  -ei*;cieot,  A,,  is  obtained  from. 
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2  sine 
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4 . 23 

Equat i ons 

4 . 22  and  -  .  '  •  ■ 

extensions  to  the  equations 

presented 

by  Wong, 

Duncan  and  .ret 

’’  j d e  an 

effective  cohesion,  c. 

The 

undra i ned 

shear  s  t  r  e  n  g  t  h 

•:  o""t:  .  ted  f  rom 

Eq.  4.22  is  identical 

to  the 

undra i ned 

shear  strength 

r ted  from 

Eq .  4.16  al though  the 

two 

equations  appear  to  oe  so^ewnat  different 


Pore  Pressure  Coefficients,  A  f  . 

The  pore  pressure  coefficient  at  failure  1  s  computed  f-'om  Eq.  4  23. 
which  was  given  above.  The  pore  pressure  coefficients  computed  from  Eqs. 
4.20  and  4.23  are  identical,  although  the  equations  are.  again,  somewhat 
different  in  algebraic  form. 

EXTENSION  OF  H.  BENJAMIN'S  STUDY 

Tnis  procedure  was  also  suggested  by  Wong,  Duncan  and  Seed,  and  is 
based  on  earlier  work  by  Benjamin  (1975).  The  procedure  is  based  on  the 
assumption  that  the  ratio  of  the  undrained  shear  strength  (s  )  to 
effective  major  principal  stress  at  consol i dat ' on  (c.  )  is  the  same  for 
isotropic  and  anisotropic  consolidation. 
versus  envelopes. 

The  undrained  shear  strength  is  related  to  the  effective  stress  on 
the  ‘"allure  plane  at  consolidation  by  the  equation. 


cRf  +  °fc 


2  Kc  tan*Rf 

[(Kc  +  1)  -  (Kc  -  l)sin$] 


[1  -  cot ( 45  +  |)tan<j>Rf] 
Pore  Pressure  Coefficients,  . 


The  pore  pressure  coefficient.  A^,  is  computed  by  first  computing 
the  change  in  pore  pressure  at  failure  from  the  following  equation: 


•  U  f  -  O-  1  -  - j-r - —  +  [ - -  +  CvN-] 

f  3cl  o,  -  1)  cosp  <*> 

JC  cp 


The  corresponding  difference  in  the  principal  stress  changes. 
Ao^f.  is  computed  using  the  following  equation: 

2 1 f f /°3c 

Au  f  -  a o  f  [ - : — —  -  (K  -  1)] 

If  3f  3cL  costf  c  J 
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Finally,  the  pore  pressure  coefficient  is  computed  by  dividing  tne 
change  in  pore  pressure  at  failure  by  the  difference  in  the  principal 
stress  changes  as. 


A 


f 


Au 


f 


Ad 


If 


Ao 


3f 


4.27 


LINEAR  INTERPOLATION  PROCEDURE 

The  linear  i nterpol at i on  procedure  consists  of  i nterpol at i ng 
linearly  between  the  lines  for  isotropic  consolidation  (Kc  = 

1.0)  and  for  anisotropic  consolidation  with  =  K^.  The  envelope 
correspond! ng  to  anisotropic  consolidation  with  Kc  =  K^,  is  identical  to 
the  effective  stress  failure  envelope  determined  from 
consol i dated-drai ned  shear  tests  and  can  be  determined  either  from 
consol idated-drained  tests  or  isotropically  consol i dated-undrai ned  shear 
tests  with  pore  pressure  measurements. 
versus  c.p_  envelopes. 

Lines  corresponding  to  intermediate  consolidation  stress  ratios 
between  unity  and  the  value  at  failure  are  obtained  by  linearly 
interpolating  between  the  lines  of  versus  b^c  obtained  for  isotropic 
consolidation  (K^  =  1.0)  and  for  the  maximum  permissible  effective 
principal  stress  ratio,  i.  e.  the  value,  K^,  at  failure  in  a 
consol idated-drai ned  shear  test.  The  line  representing  the  relationship 
between  and  for  isotropic  consolidation  is  represented  by  the 
equat i on . 


_ c  *  nf|_tanoRf 

ff  1  -  (1  -  s  i  n<t )  tan-i^/costf 


4.28 


The  line  representing  the  relationship  between  and  for 
anisotropic  consolidation  is  simply  the  effective  stress  failure 
envel ope ,  i .  e . 


t  ^  ^  =  c  +  c^ctano 


4.2 


The  maximum  effective  principal  stress  ratio  for  the  effective  stress 
failure  envelope  is  related  to  the  effective  stress  shear  strength 
parameters  and  the  effective  normal  stress  on  the  failure  plane  at 
failure,  (or  at  consolidation  since  =  o^)  by  the  expression, 


.  2  /  n  c  ,  ^  \  ,  C 


=  tan  (45  +  + 


[(K,.  +  1)  -  (K  1  )si  np]  c  - 

- _5: - tan  (45  +  f)  4.3 

Or. 


f  c 


Thus,  if  the  effective  stress  cohesion  value,  c,  is  not  zero,  the 
maximum  value  of  the  effective  principal  stress  ratio  (K^)  actually 


varies  with  the  effective  normal  stress  on  the  failure  plane.  This 


complicates  linear  interpolation  in  that  a  value  of  must  be 


calculated  each  time  an  ' iterpol ati on  is  performed  depending  on  the 
value  of  the  effective  normal  stress  on  the  failure  plane.  Duncan  Wong 
and  Seed  avoided  this  problem  in  at  least  one  instance  where  the 
effective  cohesion  was  not  zero  by  ignoring  the  cohesion  and  computing 
the  effective  principal  stress  ratio  at  failure  from  the  expression. 


Kf  =  tan2(45  +  |) 


By  neglecting  the  cohesion  value  in  computing  a  slightly  cc’se-- 
estimate  is  made  of  the  shear  strength  in  the  linear  i  nte  a* 
Pore  Pressure  Coefficients,  . 

Pore  pressure  coeffirients  are  calculated  i_,  ' 


pore  pressure  at  failure  from  the  following: 


where  is  computed  by  linear  interpolation  and  corresponds  to  the 
value  of  a3c  in  this  equation.  (Note:  the  relationship  between  and 
a^c  in  Eq.  4.32  is  given  by  Eq.  4.21).  The  difference  in  the  principal 
stress  changes  is  computed  from, 

2t, 


Ao,£  -  a =  — -=  -  o-.  (K  -  1) 
If  3f  cos<}>  3cv  c 
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Finally,  the  pore  pressure  coefficient,  A^,  is  computed  from, 
au£ 


Ar  = 


f  Aolf  -  Ao3f 


4.34 


COMPARISON  OF  PROCEDURES 

Undrained  strength  envelopes  in  the  form  of  lines  representing  the 
relationship  between  and  were  computed  using  the  equations 
presented  above  for  several  different  soils.  The  corresponding 
relationships  between  the  pore  pressure  coefficient,  and  o^c  were 
also  computed  and  compared. 

Pilarcitos  Dam  Soil 

Wong,  Duncan  and  Seed  presented  shear  strength  envelopes  for  sandy 
clay  material  from  the  upstream  slope  of  Pilarcitos  Dam.  The  shear 
strength  envelopes  (R  and  S)  for  this  soil  were  curved  and  several  sets 
of  strength  parameters  were  reported,  depending  on  the  range  in  stress 
over  which  the  failure  envelopes  were  drawn.  For  low  stresses  (0-10  psi) 
the  following  values  were  reported: 

Effective  stress  envelope: 
c  =  0,  $  =  45  degrees 

R  envelope: 


Cnf  =  60  psf,  0Df  =  23  degrees 


>  *.<  i->  ».« 


>*  LI 


y| 


The  first  set  of  comparative  calculations  was  performed  using  these  data 
from  Wong,  Duncan  and  Seed. 

versus  envelopes.  Envelopes  for  anisotropic  consolidation 
stress  ratios,  K  ,  of  1,  2,  3,  4,  5,  and  5.83  (=  K^)  are  shown  in  Figs. 
4.1  through  4.6,  for  Taylor's  procedure,  Lowe  and  Karafiath's  procedure, 
the  Modified  Karafiath  procedure,  the  Extension  of  Noorany  and  Seed,  the 
Extension  of  Benjamin,  and  the  linear  interpolation  procedure, 
respectively. 

versus  a^.  The  variations  in  the  pore  pressures  coefficients, 
ftf,  with  the  effective  normal  stress  on  the  failure  plane  at 
consolidation  are  shown  in  Figs.  4.7  through  4.12,  for  Taylor's 
procedures,  Lowe  and  Karafiath's  procedure,  the  Modified  Karafiath 
procedure,  the  Extension  of  Noorany  and  Seed,  the  Extension  of  Benjamin, 
and  the  linear  interpolation  procedures,  respectively.  The  relationships 
are  again  shown  for  consolidation  stress  ratios  of  1,  2,  3,  4,  5  and 
5.83. 

Other  soils. 

V'^iations  in  x^  and  with  6^c  were  also  calculated  for  three 
additional  sets  (Set  Nos.  2,  3  and  4)  of  soil  properties.  Data  for  Set 
No.  2  were  obtained  from  Wong,  Duncan  and  Seed  and  were  based  on  data 
for  a  brown  sandy  clay  from  a  dam  site  in  Rio  Blanco,  Colorado.  Data  for 
Set  No.  3  represent  the  values  obtained  from  the  Pilarcitos  Dam  referred 
to  above  for  a  range  in  stress  form  0  to  100  psi.  The  final  set,  Set  No. 
4,  of  data  were  estimated  based  on  data  for  Hirfanli  Dam  presented  by 
Lowe  and  Karafiath  (1960b).  Data  for  the  three  additional  sets  as  well 
as  the  set  used  previously  are  summarized  in  Table  4.1. 
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Lowe  &  Karafiath's  procedure  -  Material  property  set  1 
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Extension  of  Noorany  and  Seed  -  Material  property  set  1 


Extension  of  Benjamin's  S 
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Procedure  -  Material  property  se 


at  Consolidation,  of  ,  for  Various  Effective  Principal  Consolidation  Stress  Ratios  (K 
Estimated  from  Strei^th  Parameters  for  Isotropic  Consolidation  Using  Taylor's  Procedure 
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at  Consolidation,  a,  ,  for  Various  Effective  Principal  Consolidation  Stress  Ratios  ( 
from  Strength  Parameters  for  Isotropic  Consolidation  Using  the  Modified  Form  of  Lowe 
Karafiath’s  Procedure 


Extension  of  Noorany  and  Seed  -  Material  property  set  1 
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Extension  of  Benjamin's  Study  -  Material  property  set  1 


at  Consolidation,  ,  ,  for  Various  Effective  Principal  Consolidation  Stress  Ratios  (K 
Estimated  from  Strength  Parameters  for  Isotropic  Consolidation  Using  the  Extension  of 
Benjamin's  Study 


Linear  Interpolation  -  Material  property  set  1 


TABLE  4.1 

SUMMARY  OF  FOUR  SETS  OF  SOIL  PROPERTIES  USED  TO 
COMPARE  PROCEDURES  FOR  ESTIMATING  SHEAR  STRENGTH 
ENVELOPES  FOR  ANISOTROPIC  CONSOLIDATION 


c„ 


<t> 


Property  Set 


c 


K 


Values  of  undrained  shear  strength  and  pore  pressure  coefficients 
were  computed  for  each  of  the  sets  of  soil  properties  for  a  single 
anisotropic  consolidation  stress  ratio.  The  stress  ratio  was  computed  as 
the  average  of  the  value  at  failure  (K^)  and  unity,  i.  e. 

Kf  +  1 

Kc  =  ~~~2  4'35 

The  stress  ratio  (K  )  is  shown  in  Table  4.1  for  each  set  of  data.  The 
effective  stresses,  o^c,  used  for  computations  with  each  set  of  data 
ranged  from  zero  to  values  approximately  two-to-three  times  the 
effective  stress  at  the  point  where  the  R  and  S  envelopes  intersect. 

This  range  was  judged  to  be  the  range  of  principal  interest  and  was 
selected  accordingly. 

versus  o^c  envelopes.  Envelopes  are  shown  in  Figs.  4.13  through 
4.16,  for  soil  property  Sets  1  through  4,  respectively.  For  each  set  of 
properties  envelopes  are  shown  for  the  six  procedures  discussed: 

Taylor's  procedures,  Lowe  and  Karafiath's  procedure,  the  Modified 
Karafiath  procedure,  the  Extension  of  Noorany  and  Seed,  the  Extension  of 
Benjamin,  and  the  linear  i nterpol ation  procedures. 

versus  6^.  Corresponding  variations  in  the  pore  pressure 
coefficients,  with  the  effective  normal  stress  on  the  failure  plane 
at  consolidation  are  shown  in  Figs.  4.17  through  4.20,  for  the  four  sets 
of  soil  properties. 

Discussion. 

Several  consistent  patterns  exist  in  the  failure  envelopes  shown  in 
Figs.  4.13  through  4.16.  Lowe  and  Karafiath's  procedure  and  the 
extension  of  Benjamin's  work  appear  to  consistently  overestimate  the 
shear  strength  relative  to  what  is  found  by  the  other  procedures. 


Material  Property  Set  1 
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Material  Property  Set  4 
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Material  Property  Set  3 


ane  at  Consolidation,  ofc,  for  an  Effective  Principal  Consolidation  Stress  Ratio  (Kc)  of 
27  Estimated  from  Strength  Parameters  for  Isotropic  Consolidation  Using  Several  Procedures 
terial  Property  Set  No.  3 


Taylor's  procedure  suggests  higher  strength  at  low  effective  stress  as 
compared  to  the  Modified  Karafiath  procedure,  the  Extension  of  Noorany 
and  Seed,  and  the  linear  i nterpol ati on  procedures.  However,  the  nigher 
strength  estimated  using  Taylor's  procedure  in  the  cases  illustrated  may 
result  from  the  assumption  of  a  linear  stress  path.  The  assumption  of  a 
linear  stress  path  was  made  for  the  present  study  to  simplify  the 
calculations  and  because  in  most  cases  actual  effective  stress  path  data 
were  not  readily  available.  Wong,  Duncan  and  Seed  did  not  observe  the 
higher  strengths  at  low  effective  stresses  with  Taylor's  procedure  when 
they  used  actual  effective  stress  paths  for  the  soil  (from  Hirfanli  Dam) 
which  they  examined.  Taylor's  procedure,  the  Modified  Karafiath 
procedure,  the  Extension  of  Noorany  and  Seed,  and  the  linear 
interpolation  procedures  all  produced  essentially  identical  envelopes  at 
the  higher  stresses  for  all  four  sets  of  soil  properties  considered. 

This  is  substantially  in  agreement  with  the  findings  of  Wong,  Duncan  and 
Seed  for  the  soil  from  Hirfanli  Dam. 

The  pore  pressure  coefficients  shown  in  Figs.  4.17  through  4.20 
also  show  patterns  similar  to  those  observed  for  the  failure  envelopes. 
As  expected,  Lowe  and  Karafiath' s  procedure  and  the  Extension  of 
Benjamin  produce  lower  pore  pressure  coefficients,  which  is  consistent 
with  the  higher  strengths  produced  by  these  procedures.  The  pore 
pressure  coefficients  calculated  using  Taylor's  procedure,  the  Modified 
Karafiath  procedure,  the  Extension  of  Noorany  and  Seed,  and  the  linear 
interpolation  procedures  are  all  similar.  The  deviations  observed  in  the 
failure  envelopes  with  Taylor's  procedure  at  low  stresses  are  much  less 
obvious  when  the  pore  pressure  coefficients  are  examined  and  reflect  the 


sensitivity  of  the  shear  strength  to  small  differences  in  pore  pressure 
coefficient  at  low  stress. 

Comparison  of  the  strength  envelopes  and  pore  pressure  coefficients 
with  measured  values  for  tests  where  anisotropic  consol i dati on  is  used 
would  be  very  helpful  in  establishing  the  validity  of  the  various 
procedures.  There  are  probably  sufficient  data  already  in  the  literature 
to  permit  at  least  some  such  comparisons  to  be  made;  however,  additional 
testing  would  also  be  useful.  Laboratory  tests  where  anisotropic 
consolidation  is  employed  typically  show  sufficient  scatter  to  make  it 
difficult  to  differentiate  between  results  for  one  effective  principal 
stress  ratio  and  another.  The  existence  of  such  scatter  probably 
justifies  such  approximate  procedures  as  those  discussed  in  this 
section.  However,  in  order  to  verify  the  procedures  much  more  accurate 
data  are  required  than  are  often  obtained  in  routine  testing.  The 
additional  effort  required  to  develop  such  information  is  needed  to 
verify  the  approximate  procedures. 

Regardless  of  the  validity  of  the  various  procedures  for  estimating 
the  effects  of  anisotropic  consolidation,  the  results  for  even  the  case 
of  isotropic  consolidation  illustrate  the  significant  variation  in  the 
pore  pressure  coefficient,  A^,  with  consolidation  pressure.  Procedures, 
like  Bishop's  and  Morgenstern ' s  effective  stress  procedure,  which  are 
based  on  the  assumption  of  a  constant  value  for  the  pore  pressure 
coefficient,  clearly  must  be  considered  as  approximate  procedures.  The 
data  for  all  four  soils  considered  show  that  the  pore  pressure 
coefficient  is  actually  negative  at  low  stresses  and  increases  to 
positive  values  as  the  confining  pressure  increases.  Although  Bishop's 
and  Morgenstern ' s  procedure  employs  conservative,  positive  values,  and. 


SECTION  5 


AN  ALTERNATE  APPROACH  FOR  SUDDEN  DRAWDOWN 
(Wright) 

The  effective  stress  approach  first  suggested  by  Bishop  (1954)  and 
later  employed  by  Morgenstern  (1963)  is  attractive  for  several 
fundamental  reasons:  First,  the  procedure  is  based  on  effective 
stresses.  Accordingly,  the  shear  strength  parameters  are  based  on 
effective  stresses  and  can  be  obtained  from  test  results  obtained 
employing  either  consol idated-drained  procedures  or 

consol idated-undrained  procedures  with  pore  pressure  measurements.  The 
effective  stress  shear  strength  parameters  are  generally  much  less 
sensitive  to  effects  of  sample  disturbance  and  compaction  conditions 
(dry  density,  water  content,  compaction  method,  etc.),  than  the  strength 
parameters  for  either  the  R-envelope  or  versus  o^c  relationships. 
Secondly,  Bishop's  and  Morgenstern ' s  procedure  is  based  on  the  use  of 
pore  water  pressures  which  have  a  rational  fundamental  basis  and,  at 
least  conceptually,  reflect  the  tendency  of  soil  to  change  volume  and 
develop  changes  in  pore  water  pressures  during  undrained  shear.  Most  of 
the  other  procedures  for  estimating  pore  water  pressures  do  not  reflect 
the  fact  that  pore  water  pressures  may  change  due  to  induced  shear 
stresses  during  drawdown. 

The  principal  limitations  of  Bishop's  and  Morgenstern ' s  procedure 
result  from  the  assumptions  which  are  made  to  make  the  procedure  simple. 
First,  a  single  value  of  unity  is  adopted  for  the  pore  pressure 
parameter,  B.  This  implies  either  that  A  and  the  ratio  of  the  total 
principal  stress  changes  Ao./Ao^  are  both  constant  and  independent  of 


soil  type,  or  these  two  quantities  vary  in  such  a  way  that  B  is  equal  to 
unity.  The  second  limitation  of  Bishop's  and  Morgenstern ' s  procedure  is 
associated  with  the  assumption  that  the  major  principal  stress  change  is 
equal  to  the  change  in  vertical  stress  at  the  surface  of  the  slope.  As 
discussed  earlier  in  Section  2  this  assumption  is  clearly  not  valid  near 
the  face  of  the  slope. 

In  order  to  eliminate  some  of  the  limitations  imposed  by  the 
assumptions  in  Bishop's  and  Morgenstern ‘ s  procedure  (B  =  1.0,  Ao^  =  Aoy) 
and  to  investigate  the  significance  of  these  assumptions  an  alternate 
procedure  has  been  established  and  used  as  part  of  the  current  study. 

The  alternate  procedure  employs  the  finite  element  method  to  calculate 
the  change  in  stress  in  the  slope  due  to  drawdown.  The  changes  in  stress 
are  then  used  with  appropriate  values  for  the  pore  pressure  coefficients 
(A  and  B)  to  compute  the  changes  in  pore  water  pressure  caused  by  the 
drawdown.  Finally,  the  changes  in  pore  water  pressure  are  added  to  the 
pore  water  pressures  before  drawdown  to  arrive  at  final  values  of  pore 
water  pressures  after  drawdown. 

VALIDITY  OF  SKEMPTON'S  PORE  PRESSURE  EQUATION 

A  fundamental  hypothesis  of  Bishop's  and  Morgenstern ' s  procedure, 
as  well  as  the  alternative  procedure  described  herein,  is  that 
Skempton's  (1954)  equation  and  pore  pressure  coefficients  are  valid  for 
computing  changes  in  pore  water  pressures  due  to  sudden  drawdown.  Before 
proceeding  further  with  the  discussion  of  the  alternative  procedure  it 
is  appropriate  to  first  examine  the  validity  of  Skempton's  equation. 
Skempton's  equation  can  be  written  as. 


where,  Au  is  the  change  in  pore  water  pressure,  Ao^  and  A o^,  are  the 
major  and  minor  principal  stress  changes,  respectively,  and  A  and  B  are 
pore  pressure  coefficients.  Equation  5.1  is  based  on  the  assumption  that 
the  change  in  pore  water  pressure  is  independent  of  changes  in  the 
intermediate  principal  stress,  Ao^.  Henkel  (1960)  and  others  (e.  g. 
Perloff  and  Baron,  1976)  have  proposed  that  changes  in  pore  water 
pressure  are  affected  by  changes  in  the  intermediate  principal  stress. 

To  account  for  the  effects  of  the  intermediate  change  in  principal 
stress  the  following  equation  for  changes  in  pore  water  pressure  has 
been  proposed: 

Aoi  +  Ao2  +  Ao3  /  2  "  2  2 

Au  -  - ^ -  +  a\( Ao^  -  Ac^)  +  (Ac>2  -  Ao^)  +  (Ao3  '  Ao^) 

5.2 

where  Ao^,  and  Ao^  represent  the  major,  intermediate  and  minor 
principal  stress  changes,  respectively1. 

Equations  5.1  and  5.2  are  fundamentally  different  equations  and 
different  values  will  be  predicted  for  changes  in  pore  pressure 
depending  on  the  loading  path  .  For  example,  if  the  pore  pressure 
coefficients  (A  and  a)  are  determined  from  tests  with  one  loading  path, 
e.g.  triaxial  compression,  and  then  used  in  the  respective  equations 
(5.1  and  5.2)  to  predict  the  pore  pressure  changes  for  a  different 
loading  path,  e.  g.  triaxial  extension  or  some  form  of  plane  strain 


Equation  5.2  is  written  in  the  form  for  a  saturated  soil;  a 
different  form  is  required  for  partly  saturated  soils. 

2 

In  this  case  "loading"  path  includes  the  orientation  of  the 
principal  stresses  at  consolidation,  the  orientation  of  the  principal 
stress  changes,  the  intermediate  principal  stress,  and  the  deformation 
conditions  (plane  strain,  triaxial,  etc.) 


■  ■ 


vS\V>»  v'v 


122 


deformation,  different  pore  pressures  will  be  predicted,  depending  on 
which  of  the  two  equations  is  used.  The  problem  of  predicting  ore 
pressures  and  establishing  the  validity  of  Eqs.  5.1  or  5.2  is 
complicated  by  the  fact  that  the  pore  pressure  coefficients  themselves 
may  vary  with  the  loading  path  due  to  inherent  material  anisotropy.  For 
example,  Duncan  and  Seed  (1966)  reported  values  of  A^  of  1.12  for  plane 
strain  tests  in  which  the  principal  stress  at  failure  was  vertical  and 
0.70  for  plane  strain  tests  in  which  the  major  principal  stress  at 
failure  was  horizontal.  Such  variations  in  the  pore  pressure 
coefficients  themselves  make  it  very  difficult  to  verify  that  either  Eq. 
5.1  or  5.2  is  fundamentally  more  correct.  However,  in  the  present  case 
where  the  interest  lies  in  predicting  pore  water  pressures  cue  to  sudden 
drawdown  the  problem  is  somewhat  simpler:  The  purpose  of  us  ng  either 
Eq.  5.1  or  5.2  and  the  associated  pore  pressure  coefficients  is  to 
estimate  changes  in  pore  pressure  using  a  single  value  for  the  pore 
pressure  coefficient,  A  or  a,  regardless  of  the  orientation  of  the 
principal  stress  changes.  Accordingly,  the  most  suitable  equation  (5.1 
or  5.2)  is  the  one  which  produces  a  pore  pressure  coefficient  that  is 
relatively  independent  of  the  loading  path.  The  pore  pressure 
coefficient  which  is  the  more  unique  in  value  for  different  loading 
paths  reflects  the  more  valid  pore  pressure  coefficient  for  the  present 
purposes . 

Van  Saun  (1985)  examined  values  of  the  pore  pressure  coefficients 
(A  and  a)  measured  for  different  loading  paths  from  laboratory  test  data 
for  approximately  20  soils.  In  each  case  both  triaxial  compression  and 
extension  tests  were  performed  on  specimens  at  comparable  consolidation 
pressures  and  overconsolidation  ratios.  The  soils  for  which  data  were 
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examined  encompassed  undisturbed  specimens,  compacted  specimens  and 
specimens  reconsti tuted  and  consolidated  from  a  slurry.  Specimens  ranged 
from  normally  consolidated  to  specimens  with  overconsol i dati on  ratios  of 
24.  Effective  consolidation  pressures  (o^  )  ranged  from  5  to  120  psi . 

The  values  of  Skempton's  pore  pressure  coefficient,  A,  based  on  these 
test  data  for  triaxial  extension  tests  are  plotted  versus  the 
corresponding  values  for  the  compression  tests  in  Fig.  5.1.  A  line  is 
drawn  to  represent  the  line  of  equality  for  the  two  sets  (compression, 
extension)  of  test  values.  A  similar  plot  for  the  other  pore  pressure 
coefficient  (a)  is  shown  in  Fig.  5.2.  The  actual  values  shown  in  Fig. 

5.2  have  been  multiplied  by  the  square  root  of  2;  in  the  case  of 
triaxial  shear  where  two  of  the  principal  stress  changes  are  always 
equal  Eq.  5.2  becomes: 

Au  =  Ao^  +  av2( AOj  -  Ao^)  5.3 

where , 

1 

A°a  =  ^A°l  +  &02  +  Ao3^  5.4 

The  quantity  a>/2,  thus  appears  as  a  coefficient  analogous  to  the 

coefficient  A  in  Eq.  5.1. 

Examination  of  the  data  presented  in  Figs.  5.1  and  5.2  reveals  that 
neither  A  nor  a  are  unique  for  compression  and  extension  nor  does  one 
set  of  coefficients  appear  to  be  more  unique  than  the  other.  However,  at 
the  same  time  both  coefficients  are  reasonably  similar  for  the  two 
loading  paths  (triaxial  compression,  triaxial  extension)  represented  by 
the  data  in  Figs.  5.1  and  5.2.  The  average  difference  between  the  pore 
pressure  coefficient  A^  for  compression  and  extension  tests  is  0.16;  the 
higher  values  generally  being  measured  in  extension  tests.  The  average 
difference  between  the  quantity  av 2  for  compression  and  extension  tests 


Skempton's  Pore  Pressure  Coefficients 


Figure  5.1  -  Comparison  of  Skempton's  Pore  Pressure  Coefficient 
and  Triaxial  Extension  Tests  (After  Van  Saun,  1985 


Octahedral  Pore  Pressure  Coefficients 


Henkel's  Pore  PressureCoeff icient,  a/2,  from  Triaxial 
d  Triaxial  Extension  Tests  (After  Van  Saun,  1985) 


is  0.17;  the  higher  values  generally  being  measured  in  compression 
tests.  For  practical  purposes  and  based  on  the  data  summarized  in  Figs. 
5.1  and  5.2,  Eqs.  5.1  and  5.2  and  their  associated  pore  pressure 
coefficients  appear  to  be  equally  suited  for  use  in  predicting  pore 
water  pressure  due  to  sudden  drawdown. 

Although  it  would  have  been  useful  to  examine  similar  types  of  data 
to  those  presented  in  Figs.  5.1  and  5.2  for  the  case  of  plane  strain, 
such  data  are  very  limited  and  usually  the  intermediate  principal  stress 
is  not  known.  The  intermediate  principal  stress  would  have  to  be 
estimated  to  calculate  the  pore  pressure  coefficient,  a,  from  Eq.  5.2. 

Skempton's  equation  (Eq.  5.1)  was  selected  for  computing  changes  in 
pore  pressure  in  the  present  study  because  of  the  apparently  wider  use 
of  and  experience  with  the  equation  and  the  past  history  of  use  of 
Skempton's  pore  pressure  coefficients  for  sudden  drawdown.  The  data 
shown  in  Figs.  5.1  and  5.2  suggest  that  Eq.  5.2  and  the  pore  pressure 
coefficient,  a,  also  could  have  been  used  with  probably  equal  validity. 

FINITE  ELEMENT  COMPUTATIONS 

In  the  alternative  procedure  developed  herein  the  change  in  stress 
due  to  drawdown  is  computed  first  by  applying  loads  (stresses)  at  the 
face  of  the  slope  to  represent  the  changes  in  stress  caused  by  lowering 
the  adjacent  water  level.  Changes  in  stress  are  applied  as  "tensile" 
stresses  acting  away  from  the  face  of  the  slope  to  represent  the 
unloading  (Fig.  5.3).  For  the  present  study  all  computations  were 
performed  assuming  linear  elasticity  and  using  a  constant  value  for  the 
soil  modulus  (Young's  modulus).  The  actual  value  of  modulus  was 
immaterial  inasmuch  as  only  stresses  were  of  interest;  computed 
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displacements  were  ignored.  Poisson's  ratio  was  assumed  to  be  0.49  to 
approximate  undrained  loading  (zero  volume  change). 

A  finite  element  computer  program  employing  the  "QM5"  isoparametric 
quadrilateral  element  developed  by  Doherty,  Wilson  and  Taylor  (1969)  was 
used  to  compute  the  changes  in  stress.  Changes  in  stress  were  calculated 
at  the  centroids  of  each  element  and,  accordingly  the  pore  pressures 
which  were  subsequently  calculated  were  calculated  at  the  centroids  of 
each  element. 

PORE  PRESSURE  COEFFICIENTS 

Pore  pressure  coefficients  were  assumed  to  vary  within  the  slope, 
depending  on  the  effective  stresses  which  existed  just  prior  to 
drawdown.  Pore  pressure  coefficients  were  computed  at  points 
corresponding  to  the  location  of  the  centroid  of  each  of  the  finite 
elements  which  were  used  to  compute  changes  in  stress.  Pore  pressure 
coefficients  were  computed  using  Taylor's  procedure  which  was  described 
in  Section  4  of  this  paper.  A  linear  effective  stress  path  was  assumed. 
Thus,  the  pore  pressure  coefficient.  A,  could  be  related  to  the  strength 
parameters  from  the  R  envelope  (c^  and  i>^),  the  effective  stress 
friction  angle,  0,  and  the  state  of  stress  (e.  g.  and  Kc).  For 
purposes  of  computing  pore  pressure  coefficients  the  strength  parameters 
associated  with  the  R  envelope  were  considered  to  be  those  for  an  R 
envelope  which  is  drawn  tangent  to  the  circles  of  stress,  rather  than 
through  the  points  corresponding  to  stresses  on  the  failure  plane. 
Accordingly,  the  equations  and  procedures  used  to  compute  the  pore 
pressure  coefficient  were  different  from  those  presented  in  Section  4 
for  Taylor's  procedure,  which  were  based  on  an  envelope  drawn  through 
points  correspond i ng  to  stresses  on  the  failure  plane. 
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To  compute  the  pore  pressure  coefficient  for  a  given  effective 


principal  stress  ratio  at  consolidation  (K  )  and  effective  minor 


principal  stress  (o^)  an  equivalent  effective  minor  principal  stress. 


°3cl’  ^or  the  corrssP°nc^n9  stress  path  for  isotropic  consolidation  is 


found  first.  The  equivalent  effective  minor  principal  stress  is 
determined  by  solving  the  following  equation: 

i) 


+  f_3cldl  +  d2 
°3cl  =  d3c  5, 


o,  T  ( K 
3d '  c 


J  3  c 1 d  3  +  d4 


o-  .d,  +  d- 

1  +  (KC  "  1)  - 
c  a. 


5.5 


'3cld3  +  d4 


where , 


dl  = 


d2  =  2<‘  %  •  CR  V 


(N- 


l)(N.  -  1) 

<p 


d4  =  <N,  -  1)  2cR  \ 


5.6 

5.7 
5.3 
5.9 


and . 


N  =  tan2(45  +  f) 
<t  c 


5.10 


<Jr 


?  R 

N  =  tan^  (  45  +  ~) 
<L 
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Once  the  equivalent  effective  minor  principal  stress  is  calculated,  the 
corresponding  pore  pressure  coefficient  is  calculated  from  the  following 
equation,  which  is  obtained  for  isotropic  consolidation: 


°3cldl  +  d2 


°3cld3  +  d4 
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The  minor  principal  stress  )  and  the  effective  principal  stress 


ratio  for  consolidation  (K  ),  which  are  used  to  compute  the  pore 


pressure  coefficients,  were  estimated  for  the  present  study  as  follows: 
The  effective  vertical  stress  was  calculated  first  assuming  no  head 
losses  in  the  upstream  slope:  total  unit  weights  were  used  above  the 


initial  (pre-drawdown)  water  level  and  submerged  unit  weights  were  used 
below  the  water  level.  The  effective  vertical  stress  was  assumed  to  be 


the  major  principal  stress.  The  effective  minor  principal  stress  was 
calculated  assuming  that  the  effective  principal  stress  ratio  at 
consolidation  was  equal  to  the  average  of  the  value  of  the  effective 
principal  stress  ratio  at  failure  (K^)  and  a  value  of  unity,  i.  e. 


3c  K 
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where,  represents  the  effective  vertical  stress  and, 


K  = 
c 


(1  +  Kf) 


5.14 


Au  =  BAo^  +  A(Aa^  -  Ao ^ 1 


PORE  PRESSURE  COMPUTATIONS 

Once  the  changes  in  stress  and  the  pore  pressure  coefficients  were 
determined  for  the  center  of  each  element,  changes  in  pore  water 
pressure  were  calculated  from  Skempton's  equation: 

5.15 

where  Ao^  and  Ao^  are  principal  stress  changes  calculated  by  the  finite 
element  method,  A  represents  the  pore  pressure  coefficient  determined  in 
the  manner  described  above,  and  B  was  assumed  to  be  equal  to  unity 
(assuming  a  saturated  soil).  The  initial  pore  water  pressures,  uq, 
before  drawdown  were  calculated  assuming  a  horizontal  water  surface  in 
the  slope  with  no  head  losses.  Thus, 


u„  =  y  z 
o  w  w 


5.16 


where  z  is  the  distance  from  the  free  water  surface  to  the  centroid  of 
w 

the  element.  Finally,  the  yore  pressures  after  drawdown,  u^,  were 


calculated  from. 


u,  =  u  +  ;,u 
1  o 


5.17 
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STABILITY  COMPUTATIONS 


Stability  computations  were  performed  using  the  pore  pressures  at 
the  centroids  of  each  finite  element  to  interpolate  values  along  the 
shear  surface.  The  computer  program,  UTEXAS2,  was  used  for  the 
computations  and  includes  provisions  for  automatically  performing  the 
required  interpolation  (Wright,  1986a).  Thus,  once  pore  pressures  were 
determined  at  the  centroids  of  the  finite  elements  all  additional 
computations  were  performed  using  the  existing  computer  program  for 
slope  stability  analyses. 

SUMMARY 

The  finite  element  based  procedure  described  in  this  section 
includes  several  potential  improvements  over  previous  procedures.  The 
procedure  accounts  for  the  variation  in  the  pore  pressure  coefficient, 
A,  with  effective  consolidation  pressures  and  employs  *"ational 
procedures  for  estimating  the  changes  in  total  stress  in  the  slope 
caused  by  drawdown.  A  number  of  assumptions  are  also  made  in  the 
procedure : 

1.  The  soil  is  linearly  elastic  and  has  a  constant  modulus. 

2.  Poisson's  ratio  is  constant  and  equal  to  0.49. 

3.  Skempton's  pore  pressure  coefficients  are  valid  for  computing 
changes  in  pore  water  pressure. 

4.  The  major  principal  stress  before  drawdown  is  equal  to  the 
effective  vertical  overburden  pressure  (yz,  y'z). 

5.  The  effective  principal  stress  ratio  for  consolidation  (before 
drawdown)  is  exactly  midway  between  a  value  of  unity  and  the  value 
at  failure  ( Eq .  5.14). 


6.  Taylor's  procedure  is  valid  for  estimating  the  undrained  shear 
strengths  and  corresponding  pore  pressure  coefficients  for 
anisotropic  consolidation  using  results  from  tests  with  isotropic 
consol i dati on . 

7.  There  is  no  head  loss  before  drawdown;  the  pore  pressures  are 
governed  by  a  horizontal  piezometric  surface  at  the  adjacent  free 
water  1 evel . 

A  number  of  parametric  studies  were  performed  using  the  finite 
element-based  procedure  for  comparison  with  the  results  obtained  by 
other  procedures.  The  results  obtained  with  this  procedure  are  compared 
with  the  results  from  other  procedures  in  the  next  section  of  this 


SECTION  6 


COMPARATIVE  EXAMPLES  AND  CALCULATIONS 
(Wright) 

Three  series  of  example  problems  were  selected  for  comparative 
studies  employing  the  various  procedures  described  in  the  preceding 
sections  of  this  paper.  The  first  series  of  examples  was  selected  and 
used  for  comparison  of  the  effective  stress  procedures.  This  series  was 
entirely  hypothetical  and  selected  for  comparative  purposes  only.  The 
second  series  of  examples  was  based  on  the  Pilarcitos  Dam.  Pilarcitos 
Dam  experienced  a  slide  due  to  sudden  drawdown  and,  thus,  the  actual 
factor  of  safety  could  be  assumed  to  be  close  to  unity.  The  third  series 
of  examples  was  selected  and  used  to  examine  both  the  effective  stress 
analysis  procedures  and  the  two-stage  analysis  procedures.  The  third 
series  of  examples  is  also  hypothetical,  although  the  strength 
properties  are  based  on  actual  data  from  embankment  dams. 

EFFECTIVE  STRESS  PROCEDURES 

The  first  series  of  examples  selected  for  effective  stress 
calculations  are  identical  to  a  series  of  examples  examined  by  Wong, 
Duncan  and  Seed  (1983).  A  slope  height  of  100  feet  and  side  slopes  of 
1.5:1  (horizontal : vertical )  and  3.5:1  were  used.  The  water  level  before 
drawdown  was  assumed  to  be  at  the  crest  of  the  slope;  drawdown  levels 
(L)  of  35,  70  and  100  percent  of  the  initial  water  level  (H)  were 
considered.  Effective  stress  shear  strength  parameters  were  selected  to 

correspond  to  values  of  the  dimensionless  parameter,  \  (=  yHtan<ti/c)  of 

CO) 

5,  10  and  50.  Pertinent  data  for  these  examples  are  shown  in  Fig.  6.1 
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Figure  6.1  -  Parameters  Used  in  First  Series  of  Comparative  Analyses 
Employing  Effective  Stress  Procedures 
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Two  effective  stress  procedures  were  used.  The  first  was  Bishop's 


and  Morgenstern 1 s  procedure.  In  this  case  pore  water  pressures  were 
represented  by  a  piezometric  line  coincident  with  the  face  of  the  slope 
above  the  reservoir  level,  and  coincident  with  the  reservoir  surface 
below  the  reservoir  level  (Fig.  6.2);  the  reservoir  level  represents  the 
reservoir  level  after  drawdown.  The  second  effective  stress  procedure 
was  based  on  pore  water  pressures  obtained  by  solving  the  governing 
differential  equation  for  steady-state  seepage  (Eq.  2.27).  A  finite 
element  procedure  was  used  to  obtain  the  solution.  The  finite  element 
procedure  consisted  of  first  locating  the  phreatic  surface  for  steady 
seepage  before  drawdown  using  Neuman  and  Witherspoon's  (1970)  procedure 
where  the  finite  element  mesh  is  successively  adjusted  until  the 
boundary  conditions  along  the  free  surface  (zero  pressure;  no  flow)  are 
satisfied  (Fig.  6.3a).  Once  the  free  surface  was  located,  the  boundary 
conditions  along  the  free  surface  were  set  to  zero  pressure  and  the 
heads  along  the  upstream  slope  were  set  equal  to  the  values  after 
drawdown  (Fig.  6.3b).  New  heads  and  pore  water  pressures  were  then 
computed.  The  new  heads  were  used  to  calculate  the  pore  water  pressures 
after  drawdown.  This  approach  is  essentially  identical  to  the  approaches 
Cedergren  (1948,  1967)  used  with  flow  nets  and  Desai  (1972,  1977) 
employed  with  the  finite  element  method.  However,  in  the  present  case 
drawdown  was  assumed  to  be  instantaneous  and,  thus,  partial  drainage  was 
not  considered. 

Factors  of  safety  were  calculated  using  the  computer  program, 
UTEXAS2  (Wright,  1986a),  and  assuming  circular  shear  surfaces.  The 
factors  of  safety  are  summarized  in  Table  6.1.  The  ratio  of  the  factor 
of  safety  computed  using  pore  pressures  derived  from  the  steady-state 
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seepage  solution  to  the  factor  of  safety  computed  using  pore  pressures 

obtained  by  Bishop’s  and  Morgenstern ' s  procedure  is  also  shown  in  this 

table.  In  all  cases  the  factors  of  safety  computed  using  the  pore  water 

pressures  from  the  steady-state  flow  solution  are  higher  than  those 

computed  using  Bishop's  and  Morgenstern 1 s  approach.  The  differences 

between  the  factors  of  safety  computed  by  the  two  approaches  generally 

decreases  as  the  slope  angle  decreases  (slope  ratio  increases)  and  as 

the  value  of  A  decreases.  The  largest  difference  is  approximately  29 

percent  (1.5:1  slope,  A  =  50),  while  the  smallest  difference  is  only 

approximately  4  percent  (3.5:1  slope,  A  =  5).  These  differences  are 

c<p 

remarkably  small  considering  the  very  different  bases  of  the  two 
approaches . 

The  differences  between  factors  of  safety  computed  using  pore 
pressures  obtained  by  Bishop's  and  Morgenstern ' s  procedure  and  the 
steady-state  flow  solution  may  in  some  cases  be  much  larger  than  those 
summarized  in  Table  6.1.  In  the  preceding  calculations  the  slope  was 
assumed  to  be  homogeneous  and  to  rest  on  an  impervious  foundation  as 
shown  in  Fig.  6.4a.  If  instead  the  slope  rests  on  a  pervious  foundation, 
a  significant  downward  flow  would  have  been  predicted  with  an  attendant 
significant  reduction  in  pore  water  pressures  (Fig.  6.4b).  Adams  (1982) 
showed  a  similar  effect  may  even  occur  for  an  impervious  foundation  in 
the  case  of  sloping  core  dams.  He  showed  that  for  the  Lokvarka  Dam  (Fig. 
6.5)  the  flow  is  predominately  vertically  downward  except  at  the  base  of 
the  core.  This  resulted  in  nearly  horizontal  equ i potent i al  lines 
throughout  most  of  the  core  (Fig.  6.6).  Such  horizontal  equ i potent i al 
lines  correspond  to  zero  pore  water  pressure.  Consequently,  the  factors 
of  safety  computed  in  this  case  (Lokvarka  Dam)  using  pone  pressures 


Paths  for  Drainage  of  Upstream  Slope  for 
Embankment  on  Impervious  Foundation,  and 
Embankment  on  Relatively  Pervious  Foundat 


TABLE  6.1 


SUMMARY  OF  FACTORS  OF  SAFETY  CALCULATED  BY  EFFECTIVE  STRESS 
METHODS  FOR  SELECTED  EXAMPLES 

F(Steady  Seepage) 


Slope 

Ratio 

"c<}> 

Drawdown 
( percent) 

F( Bi shop/ 
Morgenstern) 

F( Steady 
Seepage) 

F ( Bi shop/Morgens tern ) 

1.5:1 

5 

35 

2.014 

2.232 

1.108 

1.5:1 

10 

35 

1.408 

1.584 

1.125 

1.5:1 

50 

35 

0.715 

0.895 

1.252 

1.5:1 

5 

70 

1.522 

1.732 

1.138 

1.5:1 

10 

70 

1.037 

1.213 

1.170 

1.5:1 

50 

70 

0.544 

0.685 

1.259 

1.5:1 

5 

100 

1.397 

1.610 

1.152 

1.5:1 

10 

100 

0.933 

1.116 

1.196 

1.5:1 

50 

100 

0.484 

0.620 

1.281 

3.5:1 

5 

35 

3.395 

3.525 

1.038 

3.5:1 

10 

35 

2.508 

2.630 

1.049 

3.5:1 

50 

35 

1.457 

- 

- 

3.5:1 

5 

70 

2.617 

2.720 

1.039 

3.5:1 

10 

70 

1.918 

2.003 

1.044 

3.5:1 

50 

70 

1.261 

1.327 

1.052 

3.5:1 

5 

100 

2.551 

2.659 

1.042 

3.5:1 

10 

100 

1.842 

1.928 

1.047 

3.5:1 

50 

100 

1.200 

1.253 

1.044 
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based  on  the  steady-state  flow  equation  would  be  expected  to  be  much 
higher  than  those  based  on  Bishop's  and  Morgenstern ' s  procedure. 

Bishop's  and  Morgenstern ‘ s  procedure  is  equivalent  to  assuming  that  the 
flow  is  horizontal  toward  the  face  of  the  slope  and,  thus,  the 

i 

|  equipotential  lines  are  vertical. 

|  ANALYSES  FOR  PILARCITOS  DAM 

Pilarcitos  Dam  was  examined  previously  by  Wong,  Duncan  and  Seed 
(1983).  They  provide  the  following  description: 

[  "Pilarcitos  Dam  is  a  homogeneous  rolled  earthfill  embankment. 

The  crest  of  the  dam  is  about  78  feet  above  the  upstream  toe.  The 
upstream  toe  has  a  gradient  of  2-1/2  to  1  for  the  lower  58  ft.  (i. 
i  e.  to  El.  678),  and  a  3  to  1  slope  from  this  point  to  the  crest 

(El.  698).  The  water  level  was  lowered  from  El.  692  to  El.  657 
between  October  7  and  November  19,  1969,  at  which  time  a  drawdown 
|  slide  occurred." 

A  cross-section  of  the  embankment  is  shown  in  Fig.  6.7.  The  embankment 
was  considered  to  have  a  total  unit  weight  of  135  pcf.  Effective  stress 
I  shear  strength  parameters  for  the  fill  were:  c  =  0,  <j>  =  45  degrees. 

i 

:  Strength  parameters  for  the  R  envelope  were  given  as:  c^  =  60  psf,  ^  = 

.  23  degrees.  Although  the  information  was  not  given,  the  R  envelope  was 

assumed  to  have  been  constructed  to  be  tangent  to  the  circles  of  stress, 

i.e.  c  =  Cp,  rather  than  c^,  and  <p  =  rather  than  0^. 

f 

j  Factors  of  safety  were  computed  for  the  Pilarcitos  Dam  using  the 

|  following  procedures: 

s 

!  1.  Bishop's  and  Morgenstern ' s  procedure. 


2.  The  effective  stress  procedure  described  in  Section  5:  Stresses 
caused  by  drawdown  were  calculated  by  the  finite  element  method  and 


pore  water  pressure  coefficients  were  used  to  compute  the 
corresponding  changes  in  pore  pressure.  This  procedure  is 
designated  as  the  "Finite  Element  -  A"  procedure. 

3.  The  Corps  of  Engineer's  two-stage  analysis  procedure  employing 
the  bilinear  (composite  R-S)  strength  envelope. 

4.  Lowe  and  Karafiath's  procedure;  linear  interpolation  was  used, 
as  described  in  Section  4,  to  estimate  the  strength  envelopes  for 
values  of  between  unity  and  the  value  at  failure  (K^). 

5.  A  two-stage  analysis  procedure  employing  a  single  versus  o^c 
relationship  based  on  an  effective  principal  stress  ratio  of  2  for 
consolidation.  The  single  envelope  was  estimated  using  linear 
interpolation. 

The  fifth  procedure  employing  a  single  envelope  to  describe  the 
relationship  between  and  o^c  was  suggested  by  Wong,  Duncan  and  Seed 

(1983)  as  an  approximate  way  of  accounting  for  the  effects  of 
anisotropic  consolidation  and  as  an  alternative  to  Lowe  and  Karafiath's 
more  rigorous  approach.  Wong,  Duncan  and  Seed  found  the  approach 
produced  reasonable  results  and,  accordingly,  the  procedure  was  examined 
as  part  of  the  current  studies. 

Values  of  the  factors  of  safety  are  summarized  in  Table  6.2.  The 
critical  shear  surface  (circle)  for  Bishop's  and  Morgenstern ' s  procedure 
was  very  shallow,  corresponding  essentially  to  a  plane  coincident  with 
the  face  of  the  slope.  The  corresponding  critical  circles  for  the 
remaining  four  procedures  listed  above  are  shown  in  respective  order  in 
Figs.  6.8  through  6.11.  The  factors  of  safety  computed  by  all  five 
procedures  were  within  approximately  15  percent  of  unity,  which  would  in 
most  cases  be  considered  sufficiently  close  to  unity  to  indicate 


Pressures  Calculated  Using  Stresses  from  Finite  Element  Computations  and  Pore  Press 
Coefficients  Estimated  from  the  R-Envelope  Using  Taylor's  Procedure 


failure.  The  second,  forth  and  fifth  of  the  five  procedures  listed  above 
produced  factors  of  safety  that  were  in  agreement  within  2  percent  of 
one  another.  Such  close  agreement  (2  percent)  is  encouraging  in  that  the 
second  procedure  is  an  effective  stress  procedure,  while  the  fourth  and 
fifth  procedures  are  two-stage  procedures  where  the  factor  of  safety 
after  drawdown  is  actually  computed  using  total  stresses  and  undrained 
shear  strengths.  However,  all  three  approaches  should  reflect  the 
effects  of  the  pore  pressures  generated  by  undrained  shear  during 
drawdown  and  the  fact  that  the  pore  pressure  coefficient  (A)  varies  with 
the  effective  consolidation  pressure'''.  Accordingly,  such  close  agreement 
is  not  surprising. 

The  factors  of  safety  calculated  by  the  second,  fourth  and  fifth 
procedures  are  in  close  agreement  but  slightly  higher  than  unity,  while 
the  factor  of  safety  calculated  by  the  Corps  of  Engineer's  procedure  is 
somewhat  less  than  unity.  The  relationship  between  the  factors  of  safety 
calculated  by  these  various  procedures  and  a  value  of  unity  is  as 
expected:  The  Corps  of  Engineers  procedure  essentially  takes  a 
lower-bound,  conservative  approach  for  defining  shear  strengths  and 
effectively  ignores  higher  strengths  which  would  develop  at  low  stresses 
due  to  dilatancy.  The  other  approaches  (FEM-A,  Lowe  and  Karafiath,  and 
Single  t^-o^  envelope)  on  the  oiner  hand  are  based  on  the  assumption 

^Although  the  pore  pressure  coefficient  is  not  explicitly 
considered  in  the  two-stage  analysis  procedures,  the  pore  pressure 
coefficients  are  explicitly  expressed.  For  a  given  set  of  effective 
stress  strength  parameters  (c  and  p)  and  undrained  strength  envelope 
( T f f -Of  line  or  and  p^  values)  a  unique  re'ationship  between  the 

pore  prissure  coefficient  and  effective  stress  is  implied  as  discussed 
in  detail  in  Section  3. 
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that  the  soil  is  undrained  and  that  any  effects  of  dilatancy  and  the 
attendant  reduction  in  pore  water  pressures  due  to  dilatancy  can  be 
fully  included.  In  reality  conditions  somewhere  between  these  two 
extremes  may  develop.  The  consistency  and  rationality  of  the  results 
obtained  for  Pilarcitos  Dam  should  encourage  further  studies  for  other 
case  histories. 

PARAMETRIC  STUDIES 

The  final  series  of  examples  was  selected  and  used  to  perform 
analyses  employing  six  different  procedures.  The  procedures  consisted  of 
the  five  procedures  employed  for  Pilarcitos  Dam  plus  the  procedure  based 
on  the  solution  of  the  steady-state  flow  equations  as  described  for  the 
first  series  of  examples.  A  slope  height  of  100  feet  and  side  slopes  of 
1.5:1  and  3.5:1  were  used.  Four  sets  of  soil  properties  were  used  as 
shown  in  Table  6.3.  The  shear  strength  properties  shown  are  identical  to 
those  considered  previously  in  Section  4  where  the  various  procedures 
for  estimating  shear  strength  envelopes  for  anisotropic  consolidation 
were  presented.  The  initial,  pre-drawdown  water  surface  was  assumed  to 
be  at  the  crest  of  the  slope  and  complete,  instantaneous  drawdown  was 
assumed  to  occur  in  all  cases  and  in  all  procedures. 

Factors  of  safety  computed  for  each  set  of  soil  properties  (1 
through  4)  are  summarized  in  Tables  6.4  and  6.5  for  1.5:1  and  3.5:1 
slopes,  respectively.  Two  sets  of  normalized  values  of  the  factor  of 
safety  are  also  shown.  The  first  set  of  normalized  values  was  obtained 
by  dividing  the  computed  factors  of  safety  by  the  corresponding  factor 
of  safety  computed  using  Lowe  and  Karafiath's  procedure.  The  second  set 
was  normalized  by  dividing  the  computed  factor  of  safety  for  each 
procedure  by  the  corresponding  factor  of  safety  computed  by  the  Corps  of 
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TABLE  6.3 


SUMMARY  OF  FOUR  SETS  OF  SOIL  PROPERTIES  USED 
FOR  COMPARATIVE  CALCULATIONS  WITH  VARIOUS  PROCEDURES 


Property  Set  Number  1: 

Total  (saturated)  unit  weight,  y:  135  pcf 
R-Envelope  parameters:  cR  =  60  psf,  0R  =  23  degrees 
Effective  stress  strength  parameters:  c  =  0,  0  =  45  degrees 

Property  Set  Number  2: 

Total  (saturated)  unit  weight,  y:  130  pcf 

R-Envelope  parameters:  cR  =  700  psf,  0R  =  15  degrees 

Effective  stress  strength  parameters:  c  =  200  psf,  0  =  31  degrees 


Property  Set  Number  3: 

Total  (saturated)  unit  weight,  y:  135  pcf 

R-Envelope  parameters:  cR  =  300  psf,  0R  =  15.5  degrees 

Effective  stress  strength  parameters:  c  =  0  psf,  0  =  34  degrees 


Property  Set  Number  4: 

Total  (saturated)  unit  weight,  y:  128  pcf 

R-Envelope  parameters:  cR  =  1,400  psf  ’  0R  "  22’5  degrees 

Effective  stress  strength  parameters:  c  =  0  psf,  0  =  35  degrees 
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Engineers  procedure.  Several  important  observations  can  be  made  from  the 
results  summarized  in  Tables  6.4  and  6.5: 

1.  In  six  of  eight  cases  the  factors  of  safety  computed  by  Lowe  and 
Karafiath's  procedure  and  by  the  Finite  Element-fl  (FEM-A)  procedure 
were  within  15  percent  of  each  other.  In  one  of  the  cases  where  the 
differences  were  larger  (1.5:1  slope,  Property  Set  No.  1) 
significant  errors  may  have  been  introduced  in  extrapolating  and 
interpolating  pore  pressures  at  the  face  of  the  slope  by  the  FEM-fl 
procedure  and  the  results  are  questionable.  In  the  second  case 
(3.5:1  slope  -  Property  Set  No.  4)  the  cohesion  component  for  the  R 
envelope  was  substantial.  In  this  case  the  factors  of  safety 
differed  by  30  percent.  However,  the  generally  small  differences 
between  the  factors  of  safety  calculated  by  Lowe  and  Karafiath's 
approach  and  the  FEM-fi  procedures  are  again  encouraging.  These  two 
approaches  differ  in  that  one  is  an  effective  stress  approach  while 
the  other  is  a  two-stage  approach,  but  both  reflect  the  tendency 
for  the  soil  to  either  dilate  or  compress  depending  on  confining 
stress . 

2.  The  two-stage  analysis  procedure  in  which  a  single  envelope 
(t^^-o^c)  was  used  produced  factors  of  safety  which  were  within  10 
percent  of  the  values  calculated  by  Lowe  and  Karafiath's  procedure 
for  six  of  the  eight  cases  considered.  The  remaining  two  cases 
where  the  differences  were  larger  occurred  for  a  1.5:1  slope:  In 
one  case  (Property  Set  1)  use  of  the  single  envelope  resulted  in  a 
17  percent  lower  factor  of  safety;  in  the  other  case  (Property  Set 
No.  4)  the  single  envelope  resulted  in  a  factor  of  safety  which  was 
37  percent  higher  than  the  value  computed  by  Lowe  and  Karafiath's 
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TABLE  6.4 


summary  of  factors  of  SAFETY 

CALCULATED 

BY  VARIOUS 

PROCEDURES 

METHODS  FOR  FOUR  SETS  OF  MATERIAL  PROPERTIES  -  I 

Factor 

,5:i  Slope 

Property 

of 

F  / 

Set 

Procedure 

Safety 

F ( Lowe ) 

F  (  C  O  r  p  S  1 

1 

Bi shop/Morgenstern 

0.499 

0.75 

1  .13 

1 

Steady  Seepage 

0.5G1 

0.75 

1.14 

1 

FEM-A  (u>=0) 

0.321 

0.48 

C.  73 

1 

FEM-A  (u<-0) 

- 

- 

- 

1 

Corps  of  Engineers 

0.441 

0.66 

1.00 

1 

Lowe  &  Karafiath 

0.664 

1.00 

1.51 

1 

Single  T^f-d^c  envelope 

0.554 

0.83 

1.26 

2 

Bi shop/Morgenstern 

0.564 

0.77 

0.84 

2 

Steady  Seepage 

0.708 

0.97 

1.05 

2 

FEM-A  ( u>-0 ) 

0.771 

1.05 

1 . 15 

2 

FEM-A  (u<0) 

0.802 

1.10 

1 . 19 

2 

Corps  of  Engineers 

0.672 

0.92 

1.00 

2 

Lowe  &  Karafiath 

0.731 

1.00 

1.09 

2 

Single  t^-o^  envelope 

0.787 

1.08 

1.17 

3 

Bi shop/Morgenstern 

0.338 

0.57 

0.71 

3 

Steady  Seepage 

0.338 

0.57 

0.71 

3 

FEM-A  (u>=0) 

0.528 

0.89 

1.12 

3 

FEM-A  (u<0) 

0.549 

0.93 

1.16 

3 

Corps  of  Engineers 

0.473 

0.80 

1.00 

3 

Lowe  &  Karafiath 

0.591 

1.00 

1.25 

3 

Single  t^-c^  envelope 

0.587 

0.99 

1.24 

4 

Bi shop/Morgenstern 

0.313 

0.33 

0.57 

4 

Steady  Seepage 

0.316 

0.33 

0.57 

4 

FEM-A  ( u->  =  0 ) 

0.930 

0.98 

1.68 

4 

FEM-A  ( u<0 ) 

0.991 

1.05 

1.80 

4 

Corps  of  Engineers 

0.552 

0.58 

1.00 

4 

Lowe  &  Karafiath 

0.947 

1.00 

1.72 

4 

Single  envelope 

1.295 

1.37 

2.35 

Note:  All  factors  of  safety  were  calculated  using  Spencer's 
limit  equilibrium  procedure  of  slices. 


i 
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TABLE  5.5 


SUMMARY  OF  FACTORS  OF  SAFETY  CALCULATED  BY  VARIOUS  PROCEDURES 
METHODS  FOR  FOUR  SETS  OF  MATERIAL  PROPERTIES  -  3.5:1  SLOPE 

Factor 


Property 

Set 

Procedure 

of 

Safety 

F  / 

F  ( Lowe) 

F  / 

F(Corps) 

1 

Bi  shop/Morgenstern 

1.750 

1.51 

1.90 

1 

Steady  Seepage 

1.750 

1.51 

1.90 

1 

FEM-A  (lm=0) 

1.145 

0.99 

1.24 

1 

FEM-A  (u'O) 

1.145 

0.99 

1.24 

1 

Corps  of  Engineers 

0.920 

0.79 

1.00 

1 

Lowe  &  Karafiath 

1.160 

1.00 

1.26 

1 

Single  t^-o^  envelope 

1.186 

1.02 

1.29 

1 

C 

Bi shop/Morgenstern 

1.340 

0.97 

1.14 

o 

c 

Steady  Seepage 

1.397 

1.01 

1.19 

2 

FEM-A  ( u>=0) 

1.483 

1.07 

1.26 

2 

FEM-A  (u<0) 

1.483 

1.07 

1.26 

2 

Corps  of  Engineers 

1.174 

0.85 

1.00 

2 

Lowe  &  Karafiath 

1.383 

1.00 

1.18 

2 

Single  envelope 

1.412 

1.02 

1.20 

3 

Bi shop/Morgenstern 

1.181 

1.09 

1.36 

3 

Steady  Seepage 

1.181 

1.09 

1.36 

3 

FEM-A  (u>=0) 

1.237 

1.14 

1.43 

n 

O 

FEM-A  (u<0) 

1.239 

1.14 

1.43 

3 

Corps  of  Engineers 

0.867 

0.80 

1.00 

3 

Lowe  &  Karafiath 

1.088 

1.00 

1.25 

3 

Single  t^-o^  envelope 

1.134 

1.04 

1.31 

4 

Bi shop/Morgenstern 

1.158 

0.49 

0.91 

4 

Steady  Seepage 

1.159 

0.49 

0.91 

4 

FEM-A  ( u->=0) 

1.667 

0.70 

1.31 

4 

FEM-A  (u<0) 

1.669 

0.70 

1.32 

4 

Corps  of  Engineers 

1.268 

0.53 

1.00 

4 

Lowe  &  Karafiath 

2.383 

1.00 

1.88 

4 

Single  t^-o^  envelope 

2.237 

0.94 

1.76 

Note : 

All  factors  of  safety  were 

calculated  using 

Spencer ' s 

limit  equilibrium  procedure  of  slices. 


more  rigorous  procedure.  No  significant  correlations  between  the 
differences  in  the  factors  of  safety  by  the  two  procedures  and  the 
soil  properties  were  evident;  the  differences  would  be  affected  both 
by  the  effective  principal  stress  ratio  and  the  magnitude  of  the 
effective  normal  stresses  along  the  shear  surface  before  drawdown. 

A  single  -1^-o.p  envelope  may  not  be  suitable  for  some  cases  and, 
accordingly,  should  be  used  cautiously. 

3.  The  factors  of  safety  computed  by  the  Corps  of  Engineers 
procedure  ranged  from  approximately  53  percent  to  92  percent  of  the 
value  computed  by  Lowe  and  Karafiath's  procedure.  The  largest 
differences  occur  for  Property  Set  No.  4  which  exhibited  the 
largest  cohesion  intercept  for  the  R-envelope.  The  differences 
between  the  factors  of  safety  computed  by  the  Corps  of  Engineers' 
procedure  and  Lowe  and  Karafiath's  procedure  reflect  the 
conservative  nature  of  the  Corps  of  Engineers'  approach.  The  Corps 
of  Engineers'  approach  deliberately  neglects  the  higher  undrained 
shear  strengths  associated  with  dilatancy  at  low  confining 
pressures . 

4.  The  factors  of  safety  computed  by  effective  stress  procedures 
employing  Morgenstern 1 s  approach  and  the  steady-state  flow 
procedure  for  estimating  pore  water  pressure  were  generally  in 
agreement  within  5  percent  or  better.  The  only  exception  to  this 
agreement  occurred  for  the  1.5:1  slope  with  Property  Set  No.  2 
where  the  pore  pressures  computed  by  Bichop’s  and  Morgenstern 1 s 
approach  produced  a  factor  of  safety  approximately  20  percent  lower 
than  by  the  steady-state  flow  approach  Tn’s  '  a  >’  g  e  >'  difference  (20 
percent)  is  due  to  the  s  ’  ir  ■  f  ’  c  ar  t  effect1  ve  :ohes'0’'  .  a '  ue  fo»' 


Property  Set  No.  2,  which  caused  the  critical  shear  surface  to  pass 
deeper  than  in  other  cases.  As  the  depth  increases  the  differences 
between  the  pore  water  pressures  by  Bishop's  and  Morgenstern ' s 
approach  and  the  steady-state  flow  approach  would  be  expected  to 
become  more  pronounced. 

5.  Factors  of  safety  computed  by  effective  stress  procedures  based 
on  either  Bishop's  and  Morgenstern 1 s  approach  or  the  steady-state 
flow  approach  ranged  from  approximately  33  to  150  percent  of  the 
values  computed  by  Lowe  and  Karafiath's  procedure.  Similarly  the 
values  by  these  effective  stress  procedures  ranged  from  57  to  190 
percent  of  the  values  computed  by  the  Corps  of  Engineers  procedure. 
The  factors  of  safety  computed  by  the  effective  stress  approaches 
generally  tended  to  increase  relative  to  those  computed  by  either 
Lowe  and  Karafiath's  or  the  Corps  of  Engineers  procedures  as  the 
slope  became  flatter  and  as  the  cohesion  intercept  for  the 
R-envelope  decreased,  i.  e.  as  the  tendency  for  dilatancy  to  occur 
diminished.  The  relatively  large  differences  between  the  factors  of 
safety  computed  by  the  two  effective  stress  approaches  and  by  the 
two-stage  approaches  should  be  expected:  Both  effective  stress 
approaches  make  no  distinction  among  soil  types  with  regard  to  the 
tendency  for  pore  pressures  to  vary  depending  on  the  degree  to 
which  the  soil  tends  to  compress  or  expand,  while  the  two-stage 
procedures  are  strongly  affected  by  the  tendency  for  pore  pressures 
to  develop  and  affect  undrained  shear  strength. 

6.  With  only  one  exception  Bishop's  and  Morgenstern ' s  procedure 
produced  factors  of  safety  which  were  either  significantly  lower  or 
within  10  percent  of  the  values  computed  by  Lowe  and  Karafiath's 


procedure.  Accordingly,  the  procedure  seems  to  generally  meet  the 
intended  objective  of  being  at  least  on  the  conservative  (safe) 
side.  The  exception  to  this  conservatism  occurred  in  the  case  of  a 


3.5:1  slope  with  Property  Set  No.  1,  where  c^  was  relatively  small 
(60  psf)  while  <j>  was  relatively  large  (45  degrees).  In  this  case 
Morgenstern ' s  approach  produced  a  factor  of  safety  which  was  50 
percent  higher  than  the  value  computed  by  Lowe  and  Karafiatn's 
procedure  and  90  percent  higher  than  the  value  computed  by  the 
Corps  of  Engineers  procedure. 

7.  The  Corps  of  Engineers  procedure  produced  factors  of  safety 
which  were  either  lower  than  or  within  10  percent  of  the  values 
computed  by  all  other  procedures  with  only  a  few  exceptions.  The 
exceptions  occurrea  with  regard  to  the  effective  stress  procedures 
for  a  1.5:1  slope  and  two  of  the  sets  of  material  properties 
considered.  In  these  cases  the  effective  stress  procedures  produced 
much  lower  factors  of  safety;  however,  there  is  no  basis  to  judge 
that  either  the  effective  stress  procedures  or  Corps  of  Engineer's 


procedure  were  more  correct.  It  appears  likely  that  neither  one  is 
correct,  but  both  may  be  conservative  to  varying  degrees. 


SECTION  7 


SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 
(Duncan  and  Wright) 

SUMMARY 

A  number  of  procedures  exist  for  performing  slope  stability 
analyses  for  sudden  drawdown.  These  various  procedures  differ 
principally  in  the  manner  in  which  the  shear  strength  is  defined  and  can 
be  grouped  into  two  general  categories.  The  first  category  includes 
those  procedures  which  are  based  on  the  use  of  effective  stresses.  The 
effective  stress  procedures  require  that  pore  water  pressures  be 
estimated  at  the  end  of  sudden  drawdown  and  the  differences  among  the 
various  effective  stress  procedures  are  related  directly  to  differences 
in  the  procedures  used  to  estimate  pore  water  pressures.  The  second 
category  of  procedures  for  sudden  drawdown  analysis  is  two-stage 
procedure:  One  set  (stage)  of  calculations  is  performed  for  conditions 
immediately  prior  to  drawdown  and  is  used  to  estimate  effective  stresses 
and  correspondi ng  undrained  strengths;  the  second  set  (stage)  of 
computations  is  performed  for  conditions  immediately  after  drawdown  and 
employs  undrained  shear  strengths  estimated  from  the  first  set  of 
computations. 

Effective  Stress  Analysis  Procedures . 

Most  of  the  effective  stress  analysis  procedures  are  based  on  one 
of  two  approaches  to  estimate  pore  water  pressures.  One  approach  employs 
Skempton's  pore  pressure  coefficients  and  estimated  changes  in  total 
stress  (due  to  drawdown)  to  estimate  changes  in  pore  water  pressure  and 
subsequently  the  final  pore  water  pressures.  Bishop's  and  Morgenstern 1 s 


procedure  and  the  "new"  alternative  procedure  based  on  the  finite 
element  method  to  calculate  stresses,  which  was  described  in  Section  5 
of  this  paper,  represent  the  first  class  of  effective  stress  approaches. 
The  second  class  of  effective  stress  approaches  is  based  on  solutions 
to  some  form  of  the  equations  for  hydraulic  flow.  Most  of  these  are 
based  on  very  simplified  assumptions  and  neglect  the  tendency  for  soil 
to  change  volume  and  develop  pore  water  pressures  due  to  shear  stresses 
induced  by  drawdown.  Although  more  rigorous  solutions  are  possible  and 
should  lead  to  improved  results,  such  rigorous  solutions  have  generally 
not  been  used.  The  most  common  solutions  are  based  on  simply  solving  the 
equation  for  steady-state  flow  and  ignoring  transient  effects. 

An  alternative  procedure  has  been  presented  in  this  paper  as  an 
improvement  over  Bishop's  and  Morgenstern ' s  procedure  whereby  the 
variation  in  the  pore  pressure  coefficient  A^  with  effective  stress  is 
accounted  for  and  the  changes  in  stress  caused  by  removal  of  the  water 
load  at  the  slope  face  are  more  realistic.  However,  the  procedure  has 
not  been  tested  extensively  and  requires  considerably  more  effort  to  use 
than  any  of  the  other  procedures,  effective  stress  or  two-stage, 
described  in  this  paper. 

Two-Stage  Analysis  Procedures 

The  only  difference  among  the  various  two-stage  analysis  procedures 
is  in  the  manner  in  which  the  undrained  shear  strength  is  estimated  from 
the  effective  stresses  before  sudden  drawdown.  Two  principal  approaches 
have  been  used  for  this  purpose:  One  is  the  Corps  of  Engineers 
procedures,  which  adopts  essentially  a  "lower-bound"  approach  for 
estimating  the  relationship  between  undrained  shear  strength  and 
effective  consolidation  pressure.  The  second  approach  is  the  one 


suggested  by  Lowe  and  Karafiath  which  attempts  to  account  for  the 
effects  of  anisotropic  consolidation  on  the  undrained  shear  strength.  A 
third  approach  was  suggested  by  Wong,  Duncan  and  Seed  (1983)  and  was 
used  in  the  comparative  studies  in  Section  6.  The  third  approacti  employs 
a  single  relationship  between  undrained  shear  strength  and  effective 
consolidation  pressure,  which  accounts  in  an  approximate  manner  for  the 
effects  of  anisotropic  consol idation  and  is  somewhat  simpler  than  Lowe 
and  Karafiath's  procedure. 

Several  approaches  have  been  suggested  for  estimating  the  effects 
of  anisotropic  consolidation  on  the  undrained  shear  strength  of  the  soil 
without  performing  extensive  laboratory  tests  with  anisotropic 
consolidation.  Further  study  and  verification  of  the  various  procedures 
would  be  useful;  however,  it  appears  that  at  least  several  may  be  used 
with  reasonable  success  and  confidence  at  the  present  time.  Accordingly, 
it  appears  that  the  effects  of  anisotropic  consolidation  can  be  accounted 
for  at  least  in  an  approximate  manner  in  slope  stability  analyses 
without  the  necessity  for  very  extensive  laboratory  testing. 

CONCLUSIONS 

This  study  provides  a  basis  for  a  number  of  important  conclusions 
regarding  the  accuracy  and  the  degree  of  conservatism  of  procedures  for 
analyzing  rapid  drawdown  slope  stability.  The  writers  believe  that  the 
results  of  the  study  have  brought  the  state-of-the-art  to  a  point  where 
it  is  possible  to  establish  what  is  the  most  suitable  procedure  for 
analysis  of  rapid  drawdown  slope  stability  problems,  without  need  for 
extensive  further  studies. 

One  of  the  facts  brought  out  clearly  by  the  studies  in  this  paper 
is  that  effective  stress  analysis  procedures  that  use  simple  pore 


pressure  assumptions  do  not  properly  reflect  the  actual  behavior  of  the 
soils  in  the  slope.  These  procedures  (Bishop  and  Morgenstern,  Steady 
Seepage)  use  the  same  estimates  of  pore  water  pressure  after  drawdown  no 
matter  what  the  characteristics  of  the  soil  are  with  regard  to 
development  of  pore  pressures  during  undrained  loading.  As  estimated  by 
these  procedures,  the  pore  pressures  after  drawdown  are  the  same  for  all 
soils,  no  matter  whether  they  develop  high  or  low  pore  water  pressures 
during  undrained  loading.  Based  on  the  results  of  the  comparative 
analyses  it  appears  that  these  methods  may  be  excessively  conservative 
in  some  cases,  and  unconservative  in  other  cases.  At  best  they  are 
suitable  only  for  “qui ck-and-di rty"  analyses,  being  easy  to  perform  but 
unreliable  with  regard  to  accuracy. 

The  remaining  effective  stress  method  discussed  (the  FEM-A 
procedure)  resulted  in  factors  of  safety  that  were  consistently  in  good 
agreement  with  those  calculated  using  Lowe  and  Karafiath's  procedure. 
This  finding  is  very  significant.  It  indicates  that  these  two  different 
procedures  for  approaching  the  mechanics  of  the  analyses  lead  to 
essentially  the  same  results  if  they  are  consistent  with  respect  to  soil 
strength  parameters  and  pore  water  pressures.  Pore  water  pressures  are 
considered  explicitly  in  the  effective  stress  analyses  and  implicitly  in 
Lowe  and  Karafiath's  procedure.  Although  the  FEM-A  procedure  is  too 
complex  and  time-consuming  for  routine  use,  it  has  served  a  very  useful 
purpose  by  further  validating  Lowe  and  Karafiath's  procedure,  the  most 
soundly  based  of  the  available  methods  with  respect  to  consistency  of 
mechanics  and  treatment  of  soil  shear  strength. 

The  Corps  of  Engineers'  procedure  is  the  only  one  in  use  that 
excludes  components  of  strength  resulting  from  negative  pore  water 


pressures.  It  accomplishes  this  by  using  the  drained  strength  envelope 
in  the  low  stress  range  where  the  drained  (S)  envelope  is  below  the 
consol i dated-undrai ned  (R)  envelope.  This  procedure  ensures  that  the 
factor  of  safety  will  not  be  overestimated  if  sorre  drainage  occurs 
during  drawdown.  If  some  amount  of  drainage  does  occur  during  drawdown 
negative  pore  pressures  that  develop  in  the  slope  at  shallow  depths 
would  be  reduced  in  magnitude,  and  the  strength  of  the  soil  would  be 
reduced  accordingly. 

Methods  such  as  Lowe  and  Karaf i ath ' s  ana  others  that  include 
components  of  strength  resulting  from  negative  pore  pressures  give 
unconservative  estimates  of  factor  of  safety  except  for  completely 
undrained  conditions.  As  an  example,  consider  the  results  of  the 
analyses  of  Pilarcitos  Dam  shown  in  Table  6.2:  Lowe  and  Karafiath's 
procedure  results  in  a  value  of  F  =  1.14,  while  the  Corps  of  Engineer's 
procedure  results  in  a  value  of  F  =  0.84.  for  the  failure  condition 
where  the  actual  factor  of  safety  must  have  been  less  than  or  equal  to 
unity.  While  it  must  be  recognized  that  these  values  of  factor  of 
safety  might  be  biased  due  to  unknown  inaccuracies  in  evaluation  of  the 
shear  strength  parameters,  it  is  nevertheless  clear  that  unconservative 
estimates  of  safety  factor  may  be  calculated  when  shear  strength 
resulting  from  negative  pore  pressures  are  not  eliminated. 

Based  on  the  results  presented  in  this  paper,  it  is  clear  that  a 
method  of  analysis  of  rapid  drawdown  slope  stability  to  be  used  for 
design  of  embankment  dam  slopes  should  have  the  following 
characteristics: 

(1)  The  procedure  should  employ  a  measure  of  shear  strength  and 

pore  pressure  that  reflects  the  actual  properties  of  the  so^s  in 


the  slope.  Specifically,  the  tendency  of  the  soils  to  develop  low 
or  high  pore  pressures  during  undrained  shear  should  be  reflected 
in  the  shear  strength  or  pore  pressure  parameters  used  in  the 
analyses . 

While  in  principle  this  aspect  of  the  soil  behavior  can  be 
included  in  either  effective  stress  or  total  stress  analyses,  in 
practice  it  is  much  easier  and  more  practical  to  represent  it  by 
using  a  two-stage  analysis  procedure,  with  total  stress  analysis  of 
stability  after  drawdown. 

(2)  The  procedure  should  employ  a  soundly  based  technique  for 
relating  shear  strength  to  effective  consolidation  pressures.  The 
method  of  relating  the  shear  stress  on  the  failure  plane  at  failure 
to  the  effective  stresses  at  the  time  of  consol idation  used  by  Lowe 
and  Karafiath  is,  in  the  writers'  opinion,  the  most  logical  and 
fundamentally  sound  method  available  for  doing  this.  In  contrast, 
the  method  now  used  by  the  Corps  of  Engineers,  involving  use  of  the 
R  envelope,  is  fundamentally  less  sound.  Its  greatest  merit  is 
that  it  can  be  shown  to  be  somewhat  conservative  for  most  (but  not 
al 1 )  condi tions . 

(3)  The  procedure  should  not  include  strength  components  due  to 
negative  pore  water  pressures,  to  avoid  overestimating  factors  of 
safety  for  conditions  where  partial  drainage  occurs  during 
drawdown.  The  Corps  of  Engineers  method  of  using  the  drained 
strength  envelope  in  the  range  of  low  stresses  where  the  undrained 
strength  exceeds  the  drained  strength  is  a  simple  and  practical 
means  of  accomplishing  this  objective. 


A  method  with  these  characteristics  can  be  conceived  as  a 
combination  of  methods  currently  in  use.  It  could  be  described  as  being 
a  modification  of  the  method  currently  used  by  the  Corps  of  Engineers, 
the  modification  being  adoption  of  Lowe  and  Karafiath's  procedure  for 
representing  the  undrained  strength  of  the  soil,  rather  than  using  the 
R-envelope  as  is  now  done.  Lowe  and  Karafiath's  procedure  relates  the 
shear  stress  on  the  failure  plane  at  failure  to  the  effective  stress  on 
the  failure  plane  during  consolidation  and  the  effective  principal 
stress  ratio  during  consolidation.  The  new  method  could  also  be 
described  as  being  a  modification  of  Lowe  and  Karafiath's  procedure,  the 
modification  being  elimination  of  strength  components  due  to  negative 
pore  pressures  through  use  of  the  drained  strength  envelope  in  the  low 
range  of  stresses  where  the  undrained  strength  exceeds  the  drained 
strength . 

RECOMMENDATIONS 

It  is  recommended  that  a  computer  program  be  developed 
i ncoroorati ng  the  features  outlined  above,  and  adopted  for  use  in 
analysis  and  design  by  the  Corps  of  Engineers.  The  program  should  be 
tested  by  analyzing  the  same  benchmark  cases  analyzed  in  the  present 
study  to  insure  that,  as  expected,  it  results  in  factors  of  safety  that 
fall  between  the  values  computed  using  the  current  Corps  of  Engineers 
procedure  and  those  calculated  using  Lowe  and  Karafiath's  procedure. 

The  writer  believes  that  the  method  described  will  incorporate  the  best 
features  of  all  current  procedures,  being  soundly  based  with  regard  to 
treatment  of  shear  strength,  and  reliably  conservative  through  adoption 
of  the  Corps  of  Engineers  procedure  for  eliminating  strength  components 
resulting  from  negative  pore  water  pressures. 


It  is  recommended  that  the  linear  interpolation  method  be  used  to 
establish  the  relationships  among  shear  stress  on  the  failure  plane  at 
failure,  effective  stress  on  the  failure  plane  during  consolidation,  ana 
effective  principal  stress  ratio  during  consolidation.  This  procedure 
requires  only  R-bar  tests  (conventional  isotropically 
consol i dated-undrai ned  triaxial  tests  with  pore  pressure  measurement) 
for  strength  evaluation.  The  required  laboratory  testing  will,  thus.  De 
no  more  difficult  than  required  for  the  analysis  procedure  currently 
used  by  the  Corps  of  Engineers. 

It  is  recommended  that  the  appropriate  minimum  factor  of  safety  for 
use  with  the  new  method  of  analysis  is  F  =  1.00.  The  factors  supporting 
use  of  F  -  1.00  for  rapid  drawdown  conditions  analyzed  by  this  procedure 
are : 

(1)  Like  the  method  currently  used  by  the  Corps  of  Engineers,  the 
method  will  provide  a  conservative  bias  in  the  strength  evaluations 
by  eliminating  components  of  strength  resulting  from  negative  pore 
pressures . 

(2)  There  is  a  further  conservative  bias  in  the  use  of  freshly 
compacted  laboratory  test  specimens  for  strength  evaluations, 
because  compacted  soils  become  harder  and  stronger  with  age  after 
compaction.  Ignoring  this  increase  in  strength  and  stiffness,  as 
is  conventional  procedure,  results  in  inherently  conservative 
estimates  of  shear  strengths  for  compacted  fills. 

(3)  Most  critical  failure  mechanisms  for  slope  instability  during 
rapid  drawdown  are  shallow.  Thus,  most  rapid  drawdown  slides  do 
not  threaten  the  integrity  of  the  dam  and  pose  no  threat  of  loss  of 


the  reservoir.  The  consequence  of  a  rapid  drawdown  slide  in  most 
cases  is  only  the  cost  of  repair,  and  there  is  no  threat  to  nfe. 
The  positions  of  slip  surfaces  should  be  examined  critically  wnen 
the  slope  stability  analyses  are  performed,  and  any  slip  surface 
involving  possible  failure  of  a  large  portion  of  the  embankment,  or 
which  cuts  through  the  crest  and  would  lead  to  loss  of  freeboard,  or  wni 
has  the  potential  for  blocking  an  outlet  works,  should  have  a  fa;  tor  of 
safety  during  rapid  drawdown  greater  than  unity. 

In  summary,  it  is  the  writers'  opinion  that  the  studies  presented 
in  this  report  provide  an  adequate  basis  for  selection  of  a  method  of 
analysis  for  rapid  drawdown  for  use  by  the  Corps  of  Engineers.  This  new 
method  should  incorporate  the  best  features  of  the  method  currently  used 
by  the  Corps,  and  the  method  developed  by  Lowe  and  Karafiath.  It  is 
recommended  that  a  computer  program  be  develoDed  for  performing  analyses 
using  this  procedure,  and  that  the  program  a°n  tne  procedure  be  adopted 
as  standards  for  use  by  the  Corps. 
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